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ABSTRACT 

This volume, Volume I of a three-volume report, describes the 
final results and conclusions of a study of storm water management 
models. The objectives of the study were: 

- to review and modify the U.S. Environmental Protection 
Agency's Storm Water Management Mode) (SWMM) for Canadian 
condi tions; 

- to select and modify a submodel for the simulation of 
quantity and quality of runoff resulting from snowmelt; 

- to develop a data analysis submodel capable of extracting the 
SWMM meteorological data input from existing data banks; 

and , 

- to select and modify a high-speed model for simulation of 
frequency of discharges and overflows, capable of using 
output from the data analysis submodel as input. 

This study has expanded the capabilities of the original SWMM. 
Methods of aggregating the properties of individual subcatchments into a 
single equivalent catchment are described which can lead to a considerable 
reduction in data preparation and computer processing time. A generalized 
quality model was developed, based on SWMM equations and using the 
simplifications obtained in the aggregating techniques. Snowmelt 
quantity and quality models have been integrated with the SWMM. Methods 
of interfacing the SWMM and STORM (U.S. Army Corps of Engineers) models 
have been discussed and demonstrated, and a fully operational package of 
models, capable of simulation all aspects of storm water flow and quality, 
has been prepared for application to Canadian conditions. 

Additional improvements to the SWMM and STORM, and new areas 
for their use in applications such as simulation of inflow Into sanitary 
sewers, continuous simulation and rural hydrology have been suggested. 



RESUME 



Le present volume, le premier d'un rapport qui en compte 
trois, presents les conclusions et las resultats dSfinitifs d'une 
etude pcrtant sur les modules de gestion des eaux pluviales. Cette 
etude avait pour objet 

a. de considfirer le modele de gestion des eaux pluviates 
(SWMM)" de P Env i ronmental Protection Agency des Etats-Unis 
et de le modifier pour 1 'adapter aux conditions du 

mi 1 ieu canad ien, 

b. de choisir et de modifier un sous-modeie pour simuler 
quant i tat ivement et qual i tat i vement le rui ssel iement 
provenant de la fonte des neiges, 

c. de construire un sous-modeie d'analyse des donnfees qui 
puisse extraire des banques existantes les donnfees 
met^orologiques alimentant le SWMM et 

d. de choisir et de modifier un modele uitra-rapide simulant 
la frequence des 6couIements et des d^versements qui 
puisse utiliser a 1' intrant les donnfees de sortie du 
sous-module d'analyse. 

La pr^sente 6tude a permis d'accroTtre le potentiel du SWMM 
original. Le rapport indique comment rassembler les proprifetSs d'un 
certain nombre de sous-col lecteurs en un seul collecteur Equivalent; 
cette m6thode pourrait redu i re consi d^rablement le temps necessaire 
3 la preparation et au traitement des donnfees. Un modele qualitatif 
general a 6t6 mis au point ^ partir des Equations du SWMM et grace aux 
simplifications fournies par les tecliniques de rassemblement . Des 
modules quantitatifs et qualitatifs des eaux de fonte ont pu s'int^grer 
au SWMM. On a examine et dfecrit certaines m^thodes permettant de reller 
le SWMM et le modele STORM (U.S. Army Corps of Engineers). Par ailleurs, 
on a pr6parfi un ensemble de nxjddles entierement fonct lonnels , part i cul ifire- 
ment adaptfes milieu canad ien et qui peuvent simuler tous les aspects de 
l'6couIement et de la quality des eaux pluviales. 



Storm Water Management Model 



ii 



On a propose certaines ameliorations au SWMM et au modele 
STORM ainsi que de nouvelles appii cat ions qu'on peut en fairs comme, 
par exemple, S la simulation des d€versements pluviaux dans les systemes 
d'eaux us€es, S la simulation continue et d I'hydrologie rurale. 
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SUMMARY 

This report documents the findings and conclusions of a study 
for the development of a storm water management model. The terms of 
reference for this study were established by Environment Canada and 
the Ontario Ministry of the Environment under Contract OSSit-00^6 
dated May 30, 197^- The study was conducted jointly by Proctor £ 
Redfern Limited and James F. MacLaren Limited. 

The objectives of the study were: 

(a) to review and modify the U.S. Environmental Protection 
Agency Storm Water Management Model (SWMM) for Canadian 
condi tions; 

(b) to select and modify a submodel for the simulation of 
quantity and quality of runoff resulting from snowmelt; 

(c) to develop a data analysis submodel capable of extracting 
the SWMM meteorological data input from existing data 
banks; and, 

(d) to select and modify a high speed irwdel for simulation 
of frequency of discharges and overflows, with the 
capability to use output from the data analysis submodel 
as Input. 

The report documentation of this study Is contained In three volumes: 
Final Report (Volume I), Technical Background (Volume II), and a User's 
Manual (Volume III). Program listings will be available from the Water 
Resources Branch, Ontario Ministry of the Environment. 

This initial summary presents a concise synopsis of Volume 1 
and serves as an introduction to the study. 

I. BASIC CONCEPTS 

Several factors contribute to the need to review traditional 
methodologies used in municipal drainage design and pollution control: 

(a) New developments are resulting in increased flood 
peaks and extensive downstream flooding. 

(b) Measurements have shown that the concentrations of sus- 
pended solids in surface runoff may exceed those recorded 
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for sanitary sewage, while BOD levels are generally 
similar to those in the effluent from a secondary treat- 
ment facility. Consequently, the effects of storm runoff 
must be considered in the overall cost benefit analysis 
when pollution control measures, such as combined sewer 
separation or the upgrading of treatment plants, are 
considered. in some cases, the storage and treatment 
of urban runoff may be more cost effective than other 
pollution control measures. 
(c) Urban runoff models have been developed with the capability 
of simulating the quantity and quality of storm water flows. 
These models generate overland flow hydrographs for sub- 
watersheds. The accumulation, and washoff, of surface 
pollutants is computed, and the resulting hydrographs 
and pol lutographs are routed through a schematic represen- 
tation of the actual or proposed sewer system. Because 
these models operate with hydrographs, and not merely 
with peak flows, the effects of storage may be considered 
at various surface and subsurface locations. Some models 
are capable of simulating the effect of different treatment 
processes on the storm or combined flows and, subsequently, 
assessing the effects of discharges and overflows in the 
receiving water. 

2. SCOPE OF THE STUDY 

Engineers intending to use modelling techniques face the task 
of selecting the most appropriate model for the particular phase of the 
work under consideration. Once a good model has been selected, it is 
necessary to apply it in a useful and sensible manner, bearing in 
mind its inherent limitations. 

An analysis of different urban runoff models and a review of 
previous studies in Canada, the U.S. and Europe indicated that one of 
the most comprehensive and best documented one-event simulation models 
is the Storm Water Management Model (SWMM) developed for the U.S. Environ- 
mental Protection Agency. The terms of reference of this study required 
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that this model be rendered fully operational and subsequently calibrated 
using the results of measurements recorded in various Canadian municipa- 
lities. The sensitivity of the model to various input parameters was 
to be investigated and methods of aggregating various input data in 
order to simplify data preparation and reduce computer processing time 
were to be proposed. 

The routing or TRANSPORT routine in the SWMM was to be assessed 
in comparison with the more sophisticated computations of the Water 
Resources Engineers (WRE) and Dorsch Hydrograph Volume Method (HVM) 
models, and the appropriate model for particular types of analysis 
ident i f ied. 

It was also required to revise the STORAGE and TREATMENT blocks 
of the SWMM in order to better reflect Canadian conditions. 

The other routines of the model, "Inflow-Infiltration" and 
"Receiving Water" blocks, were to be assessed on the basis of literature 
reviews and general applicability. 

The importance of snowmelt in the simulation of the quantity and 
quality of runoff was recognized as being of particular importance in 
Canada. The SWMM model did not consider snowmelt; therefore, a snowmelt 
routine was selected for integration with the SWMM. 

A data analysis model, capable of processing meteorological 
data directly from the Atmospheric Environment Service data bank into 
a form suitable for input to various urban runoff models was also 
developed. 

The STORM model developed for the U.S. Army Corps of Engineers 
was selected for the simulation of frequencies and magnitudes of 
overflows, and tested on different watersheds. 

3. ASSESSMENT AND SELECTION OF CANADIAN WATERSHEDS 

An initial review of the data requirements of urban runoff 
models and available measurements from urban watersheds in the U.S. 
and Canada revealed that, at the time of the study, the existing data 
from any one watershed was not sufficient for all model testing purposes. 
Since there were very few fully integrated sets of measurements of quantity 
and quality and rainfall data, it was necessary to use data originating 
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from different areas for different aspects of the studies. The adequacy 
of data from Canadian watersheds was found to be comparable with that 
from U.S. watersheds. Canadian drainage areas simulated were: 

(a) the 5't2 acre Bannatyne combined sewer district in 
Winnipeg; 

(b) the 2330 acre West Toronto combined sewer district; 
and, 

(c) the 48 acre Brucewood separate sewer district in North 
Yorit, Toronto. 

There were data voids in all areas, but a certain number of events were 
isolated for which there was reasonable confidence in the measurements. 
The report contains detailed descriptions of these study areas and 
includes a discussion of possible sources of error. More recent data 
acquisition efforts are also briefly reported. 

4. SWHM FLOW SIMULATIONS FOR SELECTED CANADIAN WATERSHEDS 

A review of previous studies in Canada and other countries 
indicated that SWMM flow simulations compared reasonably well with 
measurements and results of other urban runoff models. The only watershed 
previously used for model testing in Canada was about 100 acres in size. 
Studies presented in the report for areas ranging in size from k8 to 2230 
acres indicated that about two-thirds of measured peak flows may be 
simulated by SWMM within +^ 20 percent. 

A sensitivity analysis of the various parameters and default 
values confirmed the suitability of the program defaults in most 
appi I cat ions. 1 

In the case of a system under surcharge, this work indicated 
that the shape of the SWMM computed hydrograph is distorted. Surcharged 
conditions usually occur fairly infrequently and so this weakness In 
the model will not be particularly significant in studies for pollution 
abatement. However, a more comprehensive analysis is required during 
the design of new drainage or relief sewers. 

Simulations of surcharged systems conducted using the HVM 
and WRE models indicated a superior simulation compared to the original 
SWMM simulations. Since inlet hydrographs generated by the three models 
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were very similar, the Improvement resulted from the more sophisticated 
pipe routing techniques used in the HVM and WRE models. The models were 
compared on the basis of the underlying theory, as well as data and computer 
processing requirements. The performance of the WRE model was found 
to be equivalent to that of the HVM. According to information from EPA, 
an updated version of the WRE model will be incorporated in a new 1976 
version of the SWMM, as well as the snowmelt routine developed In this 
report. 

The principal conclusion was that the most recent version of the 
SWMM is an appropriate tool for a broad range of studies, from preliminary 
planning to practical design work (with the exception of design and 
analysis of systems under surcharge). Further field measurements for 
calibration of the model soieiy for flow simulation may be required 
only in special cases. 

5. SWMM QUALITY SIMULATIONS FOR SELECTED CANADIAN WATERSHEDS 

A comprehensive review of quality modelling and available field 
data showed that the simulation of the quality of storm related flows was 
not as accurate as flow simulation. The simulations of storm water 
quality conducted for the Brucewood watershed and for other areas 
emphasized that calibration with meiasurements is required for an 
"order of magnitude" agreement. In the Bannatyne combined sewer district, 
where sewer grades are very flat, the scouring of extensive depositions 
in the trunk lines is a major factor affecting the quality of the 
combined overflow. These phenomena were reasonably well simulated by the 
SWMM. 

The two optional methods included in the SWMM for suspended 
solids computations were compared and, while neither seems to be comprehen- 
sive, some recommendations for the use of the appropriate method in 
different circumstances are given. 

Sensitivity analyses were conducted for variations in the 
input parameters describing pollutant accumulation, such as the length of 
the antecedent dry period, cleaning frequency, specific gravity of 
sediments, washoff and accumulation rates, and sewer grade. The basic 
algorithm used in the SWMM quality simulation appears to account logically 
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for these variations. It seems that the present state-of-the-art in 
this area does not justify sophisticated pollutant routing or the considera- 
tion of pollutant decay while in transit. It appears, therefore, that 
calibration may be facilitated by a more simple approach as described 
in Chapter 10. 

It is concluded that since the model responded in a reasonable 
manner to alternative input conditions, it may be used within the 
documented limitations to compare various control alternatives. 

6. ASSESSMENT OF ADDITIONAL SWMM ROUTINES 

A review of some Canadian data and of the methodology used 
in the SWMM infiltration subroutine highlighted the need for additional 
infiltration measurements. The application of SWMM to the study of 
sanitary sewer inf low- inf i 1 trat ion is briefly discussed. 

The SWMM Receiving Water block was debugged and applied in 
several hypothetical studies. A theoretical comparison with more 
sophisticated models is presented. The model appears to be most useful 
for planning applications, in which different alternatives have to be 
assessed under conditions of limited data. 

7. STORAGE/TREATMENT ROUTINES 

A review of the literature regarding treatment of storm runoff 
and overflows is presented in Volume II of the report. Following an 
analysis of the model and the literature review, several changes in 
the treatment model were made to reflect actual operating practices. 
Modifications have been included In by-passing of flows, modelling of 
solids settling in storage, and the underflow in the swirl concentrator. 
The revised model provides a subroutine designed to estimate the cost of 
treatment units based on the Canadian market. 

Each treatment unit in the model was tested separately and the 
performance of a combination of storage and treatment units for a 
particular storm is illustrated. The limitations and potential uses 
of the model are discussed. 
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8* SNOWMELT QUANTITY AND QUALITY ROUTINES 

Selection of a relatively simple snowmeit model resulted from 
an extensive review of the literature, presented in VolLime It. The nx>del 
was modified and integrated into SWMM and is now available as a user 
option. Required input includes hourly temperature and wind velocity 
information and a description of ground snow distribution and equivalent water 
depth. Preliminary applications of the model using published measurements 
and the limited amount of data available for the Brucewood district 
indicated satisfactory operation. 

The available literature describing the accumulation of contaminants 
in the snow layer and, in particular, their pollution of the snowmeit 
runoff is rather limited. At this stage, an equation similar to that 
used in the model for surface runoff is recommended for the estimation of 
lead and chloride concentrations in the runoff from snowmeit events. 
Comparisons with measurements in Brucewood showed that the model can 
be calibrated if data are available. 

9. A SIMPLIFIED METHOD FOR SWMM MODELLING-AGGREGATED CATCHMENTS 

Simplified procedures leading to reductions in computer processing 
time and data collection requirements were developed. The essential 
problem when using less detail in modelling an area is that more conduits 
are omitted from the analysis and, consequently, less conduit storage is 
available in routing of the storm flows. Compensating surface storage 
may be introduced to malce up for the lost conduit storage, and the 
characteristics of the outflow hydrographs may be preserved even when 
very coarse discretizations are used. 

These essential concepts were tested on four areas: West 
Toronto, Bannatyne, and two hypothetical areas of different sizes. The 
results of these applications for both quantity and quality of runoff 
compared very well with those of the equivalent detailed simulations. 
Also described are a number of simulations based on larger time steps 
for the rainfall input and flow computations. 

The methods recommended for the aggregation of subcatchments 
significantly reduce the costs of data preparation and computer processing 
t ime. 
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10. A GENERALIZED SWMM QUALITY MODEL 

The basic concepts used are similar to those used in SWMM. 
However, this new rradel has been formulated in such a way as to permit 
any input hydrograph to form the basis for quality computations. The 
input flows can originate from measurements or any urban runoff model. 
This generalized quality model does not route pollutants through a pipe 
system and, consequently, requires significantly less computer time than a 
model which includes routing. The independent calibration of all quality 
parameters is permitted, so that the model may be rapidly adjusted to 
local conditions. An overall description of the model is included, as 
well as some verification of the model using the Brucewood data. The 
results are simitar to those obtained using the more complex SWMM. 

11. CONTINUOUS SIMULATION MODELS 

A review of the available continuous long-term simulation models 
led to the selection of STORM for this purpose. The main reasons for 
this choice were the widespread acceptance and simplicity of the model. 
STORM was applied for the West Toronto and Bannatyne areas. The occurrence 
of overflows simulated by the model agreed well with overflow records 
in these districts. The volumes of overflow predicted agreed well with 
measurements and those computed by the SWMM. The model should be used 
in the planning stage of a project. I 

The use of the model for assessing the effects of urbanization 
is briefly discussed and the use of more recent values for rural pollutant 
accumulation is examined in Volume II. The potential for application of 
the SWMM as a quas i -continuous model is also explored. The use of a 
simplified SWMM for the analysis of a considerable number of significant 
events would permit a more complete evaluation of storage and treatment 
than al lowed in STORM. 

The interfacing of the Data Analysis Program (DAP), STORM 
and SWMM in a planning study is illustrated by means of a hypothetical 
test case. STORM is used for a preliminary screening of the events 
and for critical period identification. The detailed simulation of 
individual events is conducted using the SWMM. Further refinements 
in an actual study of an existing system might involve the use of the 
WRE routing routine for the analysis of surcharge. 
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12. DEVELOPMENT OF METEOROLOGICAL DATA ANALYSIS AND PROCESSING PROGRAM 

The large volume of meteorological data required for continuous 
simulation makes computer processing of the input data highly desirable. 
Consequently, a data analysis program was developed to read hourly 
precipitation and daily maximum and minimum temperatures directly 
from the data bank of the Atmospheric Environment Service of Environment 
Canada. The model also screens the precipitation and temperature data 
for errors and constructs a complete record for each station, employing 
data from different sources to fill data voids. The precipitation records 
of up to six stations can be combined, and the computed precipitation 
record tested for consistency. Finally, the processed data is punched 
out in a format suitable for direct input to STORM. Short interval 
rainfall data can also be read and processed. 

13. GENERAL CONCLUSIONS AND RECOMMENDATIONS 

The comparison of the SWMM with measurements for several sewer 
districts, as well as a detailed consideration of the theoretical basis of 
the model, indicated that it simulates flow quantities with a satisfactory 
accuracy over a wide range of conditions. The use of the model for 
flow quantity simulation may, therefore, be recommended for planning and 
design studies related to all aspects of storm water management. 

At the present time, the SWMM simulation of flow quality may 
also be used for evaluation of the pollutional load associated with 
urban runoff and combined sewer overflows. If used within its limitations 
the model is a reasonable tool for comparison of alternatives for 
pollution abatement, including comparison of the costs of different 
combinations of storage and treatment facilities, and the effects of 
discharges in the receiving water. 

Although additional flow measurements are no longer required for 
the substantiation of SWMM flow simulation routines, quality measurements 
will generally be required for the calibration of the quality routines 
in the design phase of a pollution abatement study. Until a comprehensive 
data bank of reliable quality measurements is established, it will not 
be possible to simulate pol lutographs with a consistently good accuracy. 
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The study has also indicated that the number and quantity of 
overflow may be approximately evaluated by the simple STORM model. 
This model should be used in the planning phase of a study for an indication 
of the changes in runoff quality associated with various alternative 
developments. 

This study has expanded the capabilities of the original SWMM 
mode). Methods of aggregating the properties of individual subcatchments 
into a single equivalent catchment can lead to a considerable reduction 
in data preparation and computer processing time. A generalized 
quality model has been developed, based on the SWMM equations and 
employing the simplifications obtained using the aggregating techniques. 
Snowmelt quantity and quality models have been integrated with the 
SWMM. Methods of interfacing the SWMM and STORM models in a planning 
study are discussed and demonstrated, and a fully operational package 
of models capable of simulating all aspects of storm water flow and 
quality is now available for application in Canadian conditions. 

Additional improvements to the SWMM and STORM, and new areas 
for their use in applications such as the simulation of inflow into sanitary 
sewers, continuous simulation and rural hydrology, have been suggested. 
It is probable that both models will continue to be improved because 
of their widespread use and, in particular, through the efforts of the 
EPA and U.S. Army Corps of Engineers. Analysis and discussion of other 
urban runoff models indicate that there are a number of adequate tools 
available and these will probably be used in parallel to the SWMM and 
STORM. Modellers with experience in the application of the SWMM and 
STORM will have become familiar with the use of urban runoff models and 
could rapidly become accustomed to other models. 

Although further research should be carried out in the 
areas of model improvement and data collection, it is considered that 
one of the most productive areas for new research is the investigation, 
by means of the models, of storm water management methods to be 
implemented in practical applications. 
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1. BASIC CONCEPTS OF STORM WATER MANAGEMENT MODELLING 
1.1 General 

The expenditures for drainage works and pollution control 
facilities are among the largest items in the budgets of most municipali- 
ties, and represent a significant percentage of Federal and Provincial 
funding of public works. Design and planning procedures firmly based on 
the fundamental processes governing the quantity and quality of urban 
runoff flows will result in the most cost effective solutions to the 
problems facing planners and decision makers in most urban areas. 

The widespread access to large computers and the recent 
instigation of sampling programs have led to the development of urban runoff 
models, calibrated and validated by comparisons with field data. 
Consequently, there is a definite trend towards the replacement of the 
empirical methods of the last century, such as the Rational Method, 
which are still predominant in Canada, by a more advanced design metho- 
dology. 

The Canadian Urban Drainage Program is administered by the 
Storm and Combined Sewers Subcommittee, working in the frame of the 
Canada-Ontario Agreement on Great Lakes Water Quality, with the following 
principal objectives [I]: 

(a) to define the magnitude of pollution due to storm 
water in the Great Lakes basin; 

(b) to establish priorities and schedules for studies directed 
towards potential solutions; and, 

(c) to develop a strategy for implementing solutions. 

The scope of the present study is to make available to 
Canadian practice a package of models, tested on typical Canadian urban 
watersheds. These must be capable of simulating urban runoff under 
the wide range of meteorological conditions in Canada. The two 
principal models analyzed in this study were the U.S. Environmental 
Protection Agency's SWMM (Storm Water Management Model) and the U.S. 
Army Corps of Engineers' STORM [2, j] . 

This chapter is a review of the basic concepts of modelling. 
The two principal models are also briefly described. 
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1 , 2 The Need for Models - Recent Applications 

The need for a comprehensive approach to the simulation of 
flow quantities and the limitations of the Rational Method have been 
recognized since this method was proposed at the end of the last century. 
For urban areas, however, the gradual replacement of this method by 
methods based on hydrological principles began only in the later 1950's. 
Hydrograph methods had been introduced much earlier in rural hydrology. 
The principal reasons for the time lag appear to be the lack of rainfall 
and flow measurements and the fact that expenditures for the installation 
of storm sewers and culverts were. In the past, less significant than in 
other areas of water resources. 

The first uses of hydrologic models for urban flow simulation 
followed the development of the RRL Model (Road Research Laboratory) 
in the U.K. and the Chicago Model in the U.S. [A, 5]. The Dorsch 
HVM model has since received widespread acceptance in West Germany 
[6]. Many models have been developed in the U.S., such as the EPA's 
SWMM, the WRE model, the University of Cincinnati model. ILLUDAS, 
MIT, HYOROCOMP etc. and are described in recent literature [7]- 
Other models such as the QUURM are being developed by Canadian universities 
[8]. I 

Because of the multiplicity of urban runoff models available, 
several comparative studies have been published, in Australia [11], the 
U.S. [12], and recently in Canada [I3l- The Canadian study indicated 
that for watersheds of less than 100 acres, a number of models, such as 
the RUNOFF block of the SWMM, Dorsch HVM surface routine, RRL, UCUR, 
QUURM, can simulate approximately two-thirds of measured peak flows to 
within +_ 20 percent. Larger discrepancies of +^ 50 percent were found 
with the Rational Method. ' 

Experience gained via practical applications indicates that 
there are more advantages to hydrograph modelling techniques, as compared 
to the traditional approach, than merely increased precision. The new 
methods allow for the simulation of storage, backwater effects and inter- 
connections, and provide the practitioner with a means of developing 
a more imaginative and predictable design at a considerable saving 
[15]. I 



In Canada the first large scale model applications were conducted 
in Toronto for the analysis of relief sewers by the HVM-Dorsch method. 
Subsequently, in Winnipeg, relief sewer design and other drainage studies 
by the SWMM and WRE models resulted in significant savings and the City 
of Winnipeg now recommends the more advanced methods in its Drainage 
Criteria Manual. Simulation models are also being applied for Edmonton, 
Mississauga, Nepean Township, Ottawa, Hamilton, Vaughan, Vancouver and 
at the Toronto International Airport [I't]. 

The Increasing interest in modelling is evidenced by the large 
attendances at modelling seminars and courses organized by the EPA [9]- 
SWMM applications have been conducted in Chicago, Philadelphia, Boston, 
San Francisco and elsewhere. Recently, the City of San Francisco published 
a User's Manual for the San Francisco (WRE) Model [)0]. 

The simulation of storm water quality, together with the 
computation of flow quantities, has begun more recently, largely due to 
the realization that storm sewer discharges constitute a significant 
pollutional load on the receiving waters. Suspended solids flushed 
from streets by storm water may be more concentrated than those in 
sanitary wastes, while BOD levels are often similar to those in the 
effluent from secondary treatment plants. Another problem of concern is 
pollutional load from combined sewers. It is estimated that one-half of 
the urban population of Canada is at present served by combined sewers. 
The runoff flow generated from moderate rainstorms of about 0.1 inches 
per hour can be considerably in excess of the amount of "three times 
dry weather flow" that is normally intercepted and treated at a 
municipal treatment plant. Consequently, combined sewage is frequently 
overflowed to the receiving water. The shock loads imposed on the 
receiving water body can result in short term hazards to health and 
general environmental degradation over a longer period. 

Both quantity and quality simulations are needed for the 
adequate design of drainage in new developments. Traditional drainage 
systems usually result in an increase in peak runoff flows, compared with 
the nonurbanized situation. Often upstream development can cause down- 
stream flooding, and an additional pollutant load to the stream. The 
traditional drainage design approach has been unable to deal with such 



problems, principally because the analysis of storage requires a 
consideration of the flow-time relationship, which is the basis of all 
hydrograph models. 

Today, decision makers in the field of pollution control are 
faced with several recurrent questions, which cannot be adequately 
answered without a comprehensive analysis. These include: 

- How should runoff from new developments be controlled? 

- How is the flow and water quality in the receiving 

water affected by the various alternative drainage schemes 
that might be used in new urban areas? 

- Is sewer separation or construction of road sewers economically 
and/or environmentally justified as compared with other 
options? 

- How can upstream storage control be balanced with in-system 
storage? 

- How can a balance be achieved between cleaning opera- 
tions and downstream treatment practices? 

- What are the optimum levels of treatment and storage? 

* Is advanced treatment of sanitary flows as cost effective 
as the control of urban runoff? 

Some of these and other problems could be clarified as a result 
of extensive monitoring programs. However, these programs are not 
always effective and can be very costly when conducted properly. 

Comprehensive sampling of runoff from urban watersheds in 
Canada has been initiated only recently, and it wilt be years before 
data banks covering all specific conditions can be established. By that 
time, considerable investment could be made on the basis of very 
approximate methods, such as the Rational Method, or rules of thumb such 
as interception and treatment of two to four times the dry weather flows 
and the isolated upgrading of treatment plants. 

Although relatively novel, the use of quality simulation models, 
either for a first hand analysis or in conjunction with a more limited 
measurement program, is considered very helpful in answering the questions 
posed above. Quality simulation urban runoff models have already been 
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used in a number of cities in the U.S., including Washington, San 
Francisco, Boston etc. [l6, 17, 18]. Recently, the U.S. Army Corps 
of Engineers has recommended that STORM be used in all planning studies 
concerning the environmental effects of new urban developments [20]. 
STORM is at present being used in a study, conducted by APWA, for the 
determination of the impact of urban runoff and the potential for storage 
and treatment of runoff in a large number of American and Canadian 
municipalities. STORM was used in the assessment of the environmental 
impact of new developments in Pickering, Ontario [19]. Other practical 
applications of urban runoff quality modelling were carried out in 
Hamilton and Winnipeg. An appreciation of the problems of urban runoff 
pollution has been indicated in tie recent Thames River Basin Report 
[21], Storage ponds are also under consideration as part of a master 
drainage project, involving SWMM flov; simulation, for the Township of 
Vaughan. 

At present, the applications of new storm water management 
techniques in Canada are mainly the result of the isolated efforts of 
individual municipalities, but a large scale implementation program is 
presently under consideration. The present report forms one of the 
background studies in this program. 
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Types of Mode 1 s 



The broadest definition of the function of a model would be 
"the transf (jrma t i on of an input to an output", as indicated by the 
example given in the following diagram, 
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A good model is not only a predictive tool which generates 
additional information; it also gives a better insight into processes 
being studied and generates questions that may lead towards the 
reformulation of the initial problem. 

Different urban runoff simulation models may be used at different 
stages of a project and, consequently, these operate at different levels 
of sophistication and have somewhat different objectives. Because of 
this, it is not sound to use a single model to investigate all facets 
of a problem . 

(a) Planning models are used for an overall assessment of 
the urban runoff problem, as well as for preliminary 
estimates of the cost effectiveness of different abate- 
ment procedures. Planning models may be used for the 
rating of alternative development schemes. Even if the 
results of quality simulations are not exact, they can 
indicate the most effective controls and the best areas 
for their locations. 

(b) Design/Analysis models are oriented towards the detailed 
simulation of storm events. The models may be used for 
accurate predictions of flows and (eventually) concentrations 
at any point in the urban system and, in many cases, 

into the receiving waters. These simulations should 
illustrate the manner in which different abatement 
procedures or design options affect flow and quality at 
the critical locations. Data requirements for these 
more sophisticated models are usually more extensive than 
for the planning type of model. 

(c) Operational models are used to produce actual control 
decisions during a storm event. Rainfall is entered from 
telemetered stations and the model is used to predict 
system responses a short time into the future. 

In this report, only (a) and (b) have been investigated. 
Detailed design/analysis studies usually involve the simulation 
of a limited number of single isolated events, for which the antecendent 



conditions (e.g. soil moisture, pollutant accumulation etc.) have to 
be assumed. Usually, the limitations on the number of events modelled 
are those of computation cost and data requirements. 

Planning studies for pollution abatement purposes and flood 
control require the use of a continuous simulation model. The results 
of such a model are usually presented in a form of a statistical interpretation 
of the system response over a long period. These models are generally 
founded on the same basic fundamental principals as the single-event models. 
Some are high-speed screening models, which are necessarily fairly simple 
in order to minimize computer and data requirements. Screening models may 
eventually be coupled with optimization models in which one or more 
parameters (e.g. capital cost, instream oxygen, etc.) are optimized and 
a decision area is defined. 

Most models used in urban runoff quantity simulation are 
deterministic, (i.e., follow the basis physical phenomena) while others 
are of the so-called "black-box" type, such as unit hydrograph models. 
Some quality models are of a correlation type and these models usually 
only operate successfully when used on the watershed for which they 
were developed. The deterministic models (such as SWMM and STORM) 
have a greater transferability, because of their more fundamental formulation. 
However, all models include parameters which originate from statistical 
studies or reflect some calibration or fitting to measurements. Before 
using a model the potential user should investigate past studies to 
ensure that the model is suitable for the intended application. 

The most efficient model is not always the most sophisticated 
available. More sophisticated models require more voluminous data, 
which may not be available in the initial or planning stages of a study. 
As the framework of the study is more firmly established, and definite 
criteria for pollution and/or hydraulic control are agreed upon, the 
use of the more complex model will be necessary. 

The review of models for potential use in Canadian conditions 
indicated that, in the category of "one-event simulation models", 
the Storm Water Management Model (SWMM), developed by the U.S. Environ- 
mental Protection Agency, was one of the most promising. It is the 
only currently available model integrating all aspects of runoff quantity 



and quality simulation, treatment and effects on the receiving water 
body. The SWMM is the result of many years of effort by several leading 
organizations, including the University of Florida in Gainesville, 
Water Resources Engineers (WRE) and Metcalf and Eddy, Inc. Since the first 
version of the model was published in 1971, it has been continuously 
improved and, at the present time, there are two updated versions of 
the model available. In the near future, a new version of the SWMM will 
incorporate the sophisticated flow routing procedures developed originally 
by WRE for the simulation of surcharged conditions and complex flow 
behaviour. This version is required for more advanced design v/ork 
involving pressurized flow conditions, and may be used for design of 
relief sewers as well as new drainage studies. 

The U.S. EPA has an intensive program of technology transfer 
seminars and publications which ensure that the model has, and will 
continue to receive, widespread attention and frequent testing and 
improvement. SWMM may be used for planning and in some cases for 
des ign/analys i s . 

In the category of planning models a more approximate long 
term simulation model is also required. STORM was developed by Water 
Resources Engineers Inc. for the U.S. Army Corps of Engineers for this 
purpose. This is a simple, high-speed model, which incorporates continuous 
quality modelling routines. Because of the Corps' recommendation for 
the use of STORM in planning studies, the model is certain to receive 
widespread use. 

Brief descriptions of these models are given in the succeeding 
sections. A comparison of different features of SWMM, STORM and many 
other models is summarized In Table 1 [22j . 
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Description of SWMM 



SWMM is a comprehensive mathematical simulation model vihich 
requires a high-speed digital computer for the simulation of storm 
events. The inputs to the model include a rainfall intensity distribution 
in time (hyetograph) and data describing the idealized catchment, 
transport and receiving water systems. The principal objective of the 
model is the complete characterization of the temporal and spatial effects 



TABLE 1. MODEL CHARACTERISTICS [22] 
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of a storm [23]- To achieve this, the flow and associated pollutional 
aspects are represented as continuous curves referred to as hydrographs 
and pel lutographs . 

The urban drainage system may be illustrated as a series of 
interrelated subsystems [2h] , Figure 1. The SWMM allovjs the effect of 
each of these subsystems on the storm waters to be simulated. A simpli- 
fied overview of the model, shown in Figure 2, indicates the major processes 
performed within the model blocks. These are briefly discussed in this 
section; a complete description may be found in Volume I M of this report, 
the User ' s Manual . 

RUNOFF Block - The rainfall hyethograph is distributed evenly 
over the subcatchment within each time interval (typically five 
minutes). The drainage area is characterized by its size, 
degree of imperviousness , slope and several factors describing 
the pollutant accumulation over the area. The available 
depression storage and infiltration potential are satisfied 
before runoff occurs. The overland flow is then considered 
using the kinematic wave formula based on Manning's equation 
and continuity at each time interval. The overland flow may 
be routed through small pipes and gutters in its travel to 
the inlet manhole. The rate of overland flow determines the 
amount of the available surface pollutants incorporated in the 
flow. Thus, at the inlet, a temporal description of the flow 
and the pollutant mass washoff is available. These hydrographs 
and pol lutographs , which are the output from the RUNOFF block, 
form the input to the TRANSPORT block. 

TRANSPORT Block - This is used to represent the physical works 
that convey surface runoff hydrographs and pol lutographs from 
the contributing subcatchments to the point of discharge. This 
may be a storage treatment facility, or a receiving water. 
The inlet flows are modified by combination with flows from other 
drainage areas and by any dry weather flow in the system. The 
flows and pollutant discharge rates are attentuated during 
their passage through the pipe network to an extent depending 
on the capacity of the system and the function of In-system 
storage. 

If 
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The solution procedure follows a modified kinematic wave 
approach in which disturbances are allowed to propagate only 
in the downstream direction. The resultant flow velocity in 
each element controls the deposition or scour of solids. Flows 
may be diverted at various points in the system, either to 
overflow locations or to treatment facilities. 

STORAGE/TREATMENT Block - The output from the TRANSPORT block, 
or some portion of this, such as 3 x DWF {dry weather flow), 
may be routed to the storage units. These are specified by 
simple regular geometries with various alternative inlet and 
outlet controls. Pollutants are considered to settle out in 
the storage facility, while additional improvements in quality 
may be modelled in the variety of optional treatment units. 
The resulting attenuated hydrographs and pol 1 utographs can 
form the input to the RECEIV block. 

RECEIV Block - Any receiving water may be modelled in this 
block. The geometry of the actual system is described using 
a network of nodes and channels, or a simple river system may 
be considered as a series of channels. The output from 
either the TRANSPORT or STORAGE/TREATMENT blocks, or both, 
may be used as the input to different nodes in the receiving 
water and the subsequent responses investigated. The 
equations of continuity, motion and conservation of mass are 
solved in order to establish the variations in water level 
and pollutant concentration at nodes of interest over several 
days. The effect of different treatment policies may be 
compared in a series of runs using only the quality subroutine 
of this block. 

1 .5 Description of STORM 

Figure 3 shows a conceptual view of the urban system considered in 
STORM. The model operates with an hourly record of precipitation and 
temperature, which may extend over a large number of years. The catchment 
is described in terms of land use and impervlousness. Different pollutant 
accumulation rates are associated with the various land uses. The 
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rainfall or snowmelt in excess of the available depression storage is 
transformed directly to runoff at the outlet from the catchment. Flow 
and quality routing is not considered, the poliutograph being directly 
related to the runoff rate in any hour. A treatment rate for this 
runoff may be supplied. Flows in excess of this rate may be stored or 
considered as direct overflow. The water balance between storms is 
very simply considered as the recovery of depression storage by evapo- 
transpi rat ion , while the surface pollutant accumulation is modified 
by street cleaning at a specified frequency. For a given rainfall/ 
snowmelt record, the quantity, quality and number of overflows will vary 
with the treatment rate, storage capacity and land use. 

The model is considered in more detail in Chapter 11. 
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2. SCOPE AND ORGANIZATION OF STUDY 

2. 1 Background 

Early in 197^, the situation with regard to storm water 
management tnodelling in Canada was as follows: 

- There were many urban runoff models available, each with 
a different degree of sophistication, and suitability for 
particular tasks. 

- The urban runoff study carried out under Canada-Ontario 
Agreement Project No. 7't-8-31 [13] concluded that the 
EPA SWMM held most promise, primarily because of its 
comprehensive consideration of the urban system, accuracy, 
excellent documentation and widespread acceptance. 

- Only the RUNOFF block of the SWMM had been tested with 
Canadian data. Comparison of RUNOFF simulations with 
measurements were limited to watersheds less than 100 
acres in size. 

- There was an obvious need to verify the SWMM against 
measurements for larger watersheds and for a broader 
range of meteorologic conditions. 

- A submodel for the simulation of snowmelt quantity and 
quality was not yet available, but would certainly be 
needed by Canadian planners and engineers. 

- Perhaps due to the large number of U.S. researchers working 
on the SWMM and other urban runoff models, there were 
always several versions of each available. Not all of these 
were fully operational. It was, therefore, necessary to 
ensure that all SWMM submodels were operational and 

compat i ble. 

- To ensure that the models responded logically to reasonable 
variations of input parameters, systematic sensitivity 
analyses would be required. 

- The degree of sophistication of model input data had not 
been defined in relation to the level of accuracy required 
in the model output. 
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In general, the more sophisticated the model, the 
greater the cost of data preparation and data processing. 
The quality of results do not always increase with Increasing 
sophistication. Further investigation appeared to be 
needed in the following areas: 

a) the SWMM TRANSPORT routine required comparison with 
the more sophisticated routines used in the Dorsch HVM 
(Hydrograph Volume Method) and WRE (Water Resources 
Engineers) models, particularly with a view to detailed 
design applications. 

b) the Receiving Water Body routine in the SWMM is simpler 
than some comparable water quality models. An assessment 
of its limitations and its applicability to planning 
studies would have been useful. 

c) in SWMM, pollutants are routed through the sewer system 
and pollutant decay is modelled in transit. In most 
practical cases, the available data would not be 
sufficient either in quantity or quality to justify 
such precise simulation. Simplified routines needed 
investigation. 

There was a need for recommendations describing the appropriate 
selection of various optional routines on SWMM in specific 
c i rcumstances . 

Cost routines in the SWMM were based on U.S. cost figures 
for equipment and services. These required modification for 
Canadian conditions. 

There was a requirement for a program of testing and veri- 
fication of the STORM model, similar to that for the SWMM. 
Because of the large amount of meteorological data required, 
especially for long term simulation, a program to extract and 
manage the necessary data was required. Such a program 
should be capable of extracting information from the 
Atmospheric Environment Service data banks and processing 
these in such a format as to be suitable for direct input to the 
long term simulation model STORM, and the SWMM. 
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Terms of Reference 



The following were the terms of reference for this study: 

(1) Review and modification of the U.S. EPA Storm Water Management 
Model (SWMM) for Canadian conditions. 

The Contractor shall review the latest version of the 
EPA SWMM and modify some subsystems of the model as 
discussed in the following paragraphs. The modified 
model has to be fully compatible with the SWMM version 
maintained at the University of Florida, with respect 
to the model structure, interaction between the main 
program and subroutines, and input data requirements. 



{a; 



Executive Block 



The Executive block will be modified only to the extent 
required by addition of new subsystems, e.g., the data 

subsystem, 

(b) Runoff Block 

Runoff block generates surface runoff based on an arbitrary 
rainfall hyetograph, antecedent conditions, land use and 
topography. This block is built of seven subroutines 
which act as follows: 

RUNOFF - Liaison with EXECUTIVE BLOCK 

HYDRO - Coordinates hydrograph computation using 

RHYDRO, WSHED, GUTTER 

RHYDRO - Reads and prepares input data 

WSHED - Computes overland flow 

GUTTER - Computes flow in pipes/gutters 

HCURVE - Prepares computed hydrographs for plotting 

SFQ.UAL - Calculates pollutants washed off the land 

Only minor modifications to this block will be executed, 
namely, modifications resulting from study of the limits of 
a discretization of the watershed elements (equivalent blocks) 
and from changes to the surface water quality subroutine 
(SFQUAL). 



f| Study of Equivalent Catchments 

The level of detail necessary to adequately represent 
watershed elements shall be investigated. The objective 
of this aspect of the work will be to illustrate the 
effect of reducing the number of subcatchments on 
the accuracy of runoff simulation. While this analysis 
will be empirical, confined to testing on one or two 
watersheds, an effort will be made to define the considera- 
tions which would more generally apply to selecting the 
detail necessary to achieve the desired accuracy. 

Necessary modifications will be made to RUNOFF, RHYDRO, 
WSHED and, perhaps HYDRO and GUTTER. 

ii) Modification of the Surface Water Quality Subroutine 
(SFQUAL) 

The surface water quality subroutine was developed using, 
in some cases, limited and localized data. The contractor 
shall review the entire subroutine and shall include 
new estimates of co-efficients and model parameters, 
based on the results of ongoing Canadian as well as other 
studies. In particular, the parameters used to compute 
the following items shall be reviewed: 

Pollutants accumulated on ground surface (prior to the 

storm) 

Availability factor for pollutants on ground (during the storm) 

Estimates of dust-dirt fall 

BOD in catch basins at start of storm and removal rate 

Coliforms in surface runoff 

The surface water quality shall be described by water 
quality parameters such as: BOD, COD, DO, coliforms, 
phosphorus, nitrogen and suspended solids. Parameters may 
be added to, or deleted from the above list pending 
Contractor's review of the current studies and discussion 
with Scientific Authority- 



The Contractor shall identify the basic types of behaviour 
of the above quality parameters and processes such as those 
involving balance or decay with time. The model shall have 
flexibility to enable i t to be operated using any or 
none of these parameters within these groups for specific 
watershed studies. 

(c) Transport Block 

The Transport block routes flows in sewers, routes water 
quality, estimates dry weather flow, infiltration, and 
treats in-system storage. The block is constructed as 
f o 1 1 ows : 

TRANS Liaison with EXECUTIVE BLOCK 

TSTRDT Reads and prepares Input data 

SLOP Sequences sewer elements 

FIRST Computes constants and flow parameters 

in each flow element 
INFIL Estimates and allocates infiltration 
FILTH Computes dry weather flow quality and 

quant i ty 
DWLOAD Computes initial sediment in sewers 
INITIAL Initializes flows, areas, pollutant con 

centrations based on DWF infiltration 
QUAL - Routes pollutants through sewer 
ROUTE - Routes flow through sewer 
PRINT - Prints hydrographs and pol lutographs 
TSTORG - Computes in-sewer flow storage 

The Contractor shall concentrate his effort on the following 
three areas: 

(i) Flow routing model (ROUTE) 

(ii) Quality routing model (QUAL and DWLOAD) 

(iii) Infi Itration model (INFIL) 
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i) Flow Routing Model (ROUTE) 

This subroutine shall be reviewed and its performance 
compared to other routing models. The range of applicability 
of the present EPA subroutine will be defined as well as its 
limitations, possible improvements and increased costs 
associated with more sophisticated procedures. The 
Contractor shall recommend and implement the most 
appropriate {nonpropr ie tory ) routing scheme or schemes and 
will present the reasons for preference. The scope of 
this work shall be further reviewed and discussed with 
the Scientific Authority. 

In selecting the routing scheme, primary importance shall be 
placed on the model capabilities in handling both backwater 
and surcharge conditions, and the following control 
elements: 

Manholes | 

Pump Li ft Stat ions 

Flow Dividers (linear and nonlinear) 

In-line Storage Units 

We i r s 

The routing model shall be capable of executing the 

above calculations for conduits of various shapes frequently 

used in engineering practice. 

if) Quality Routing Model (QUAL and DWLOAD) 

The Contractor shall review the quality routing model and 
apply a sensitivity analysis to the model. The following 
model parameters will be used in the sensitivity analysis: 

Oxygen utilization rate 

Reaeration rate 

Saturation value 

Specific gravity of sewer sediments 

The above analysis may reveal if a simplification of the 
model is feasible and, if so, the appropriate simplification 
will be implemented* 
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i i i) Inf i Itration Model (INFIL) 

The Contractor shall review and test the infiltration 
model by comparing Its predictions to real data. if 
appropriate, the model wilt be simplified. 

{d ) Storage Block 

The Storage block incorporates a number of schemes for 
control and abatement of pollution due to discharges and 
overflows from sewer networks. The main functions of the 
block are as follows: 

(i) Routing of sewage flows through an off-systems 

storage 
(ii) Treatment of sewage 
(iii) Cost estimation for the schemes considered under 

I terns ( i ) and (i i } 

The Contractor shall concentrate his effort on the treatment 
model and cost effectiveness model. 

i) Treatment Model 

The heart of the Treatment Model is subroutine TREAT. It 
computes the movements and removals of water and pollutants 
on a process-by-process basis, every time step. It, in 
turn, makes use of the following subroutines: 

STRAGE - Models storage facilities 

BYPASS & 

TRLING - Link up treatment processes 

SEOIM - Models sedimentation 

KILL - Models coll form reduction 

HIGHRF - Models high rate filter operation 

The treatment process presently handled by the model is: 

Bar racks 
Fine Screens 
Sed imentat ion 
Dissolved air flotation 
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Microstrainers 
High rate f i I ters 
Effluent screens 
Chlorinat ion 

The Contractor shall update each of the process models, 
taking into account variations in treatment efficiency. 
Process models shall be developed for processes not 
included in the present model version but which are 
currently being investigated, such as the disposal of 
effluent. This will require modifications of BYPASS, 
TRUNK, TREAT and perhaps the EXECUTIVE block liaison 
subroutine STORAG which handles the Treatment Model. 

The scope of this development work shall be subject to 
review by the Scientific Authority. 

ii) Cost Effectiveness Models 

Two cost models are used in the SWMM; TSTCST in the 
TRANSPORT BLOCK, and TRCOST in the STORAGE BLOCK. Together 
they make up the Cost Effectiveness Model. 

TRCOST - Computes estimated costs for provision of 
storage and treatment facilities specified 
and operating and maintenance costs 
during storm event modelled. 

TSTCST - Computes capital and operating costs 

of the designated in-line storage units 

Both subroutines employ construction cost indexes, assumed 
land cost, assumed volume/cost relationships, etc., 
arriving at estimates. 

The Contractor shall update and relate to Canadian conditions 
all assumed costs and cost formulas built into the model. 

(e) Receiving Water Block 

The Receiving Water Model consists of two major submodels: 
a flow model (SWFLOW) and a quality model (SWQUAL) . Each 
submodel employs a number of smaller subroutines to effect 
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the specific hydraulic and quality calculations. 

The existing SWMM receiving water block will be included 
as part of the work and the advantages and disadvantages 
of acquiring or developing, and implementing other models, 
will be descr i bed. 

(2) Selection and modification of a submodel for the simulation 
of quantity and quality of runoff resulting from snowmelt. 

Under Canadian conditions, snowmelt events can be accompanied 
by intense rainstorms in spring or occasionally in late 
fall, and may contribute to flooding in urban areas, even 
in the case of small basins. In addition, in the urban 
areas of Canada, water quality during the winter may be 
significantly affected by the use of various chemicals 
on the streets for snowmelting purposes and, possibly 
to some extent, by increased air pollution due to 
heating, inefficient use of cars, etc. 

In view of this, the Contractor shall pay particular 
attention to the problem of modelling snowmelt with 
superimposed rainstorms and to the related problems of 
water qual i ty. 

Snowmelt water quantity modelling Shall consist of modifying 
existing snowmelt models to apply to urban conditions, 
including effects of chemical additives and snow removal. 
The data requirements of the snowmelt model selected shall 
match the capabilities of the Data Analysis submodel 
described below. 

Input data required may consist generally of precipitation, 
temperature, solar radiation, wind, humidity and cloud cover, 
and shall be obtained using the data bank maintained by the 
Meteorological Branch. For those cases in which not all 
the required data are available, the Contractor shall 
Indicate how these data could be estimated. 
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The effect of a rainstorm on the snowmelt conditions shall 
also be modelled, using existing approaches. 

Insofar as practical, the estimated accuracy of the 
snowmelt simulation will be consistent with the capability of 
the SWMM, developed under this contract, to predict runoff 
generating from rainstorms. 

The quality of snowmelt water shall be modelled, accounting 
for street salting and pollutants retention by snow. The 
list of water quality parameters to be considered shall 
be determined by the Contractor on the basis of results of 
current studies and discussions with the Scientific Authority. 

The snowmelt model shall be properly interfaced with the 
rest of the SWMM. 

(3) Development of a Data Analysis Submodel. 

An important aspect of the proposed study is to ensure that 
the model will be compatible with existing data banks. 
From the point of view of real storm analyses or continuous 
simulation, the most important such data bank is the 
CI imatolog ical Bank of the Canadian Atmospheric Environment 
Service. Hourly-sampled cl imatologica I data, describing 
rain events and phenomena influencing snowmelt, are 
available on magnetic tape or punched cards from the 
Atmospheric Environment Service. Subroutines shall be 
written to create and manage a data/base for the detailed 
and high-speed models from such existing data banks. 
These shall include the capability of distributing rainfall 
intensity data to different basin subcatchments . 
Recommendations for future urban rain gauge networks shall 
be made. 

The capabilities of the Data Analysis submode shall not 
be restricted to handling hourly-sampled data since, in 
future surface watershed studies, data may be available for 
shorter sample intervals. 



m 



The Contractor shall provide detailed documentation for 
all the above data analysis subroutines and shall interface 
the data submodel with the rest of the model. Formats 
of data made available through the data submodel shall 
be fully compatible with the data formats required for 
the EPA SWMM. 

The Contractor's effort will be co-ordinated with the EPA 
ongoing R&D project to establish urban hydrology data base. 

^l^) Selection and modification of a high-speed model for simulation 
of frequencies of discharges and overflows. 

The Contractor shall select and modify, where necessary, 
a high-speed model for the long-term simulation of runoff 
events. The model shall have the capability of simulating 
the runoff quantity and quality, overflows, and effects 
of storage and treatment (capacity) on the aforementioned 
phenomena. 

The model developed will be checked against the detailed 
model (SWMM) as well as against real data (e.g., southeast 
Toronto Watershed). 

The Contractor shall provide interfacing between the high- 
speed model and data files and between the high-speed 
model and the detailed model (SWMM). 

(5) Additional considerations. 

(a) Implementation of Computer 

The computer programs developed and modified as part of 

this study shall be written in FORTRAN. The total 

SWMM nxjdel, including integrated submodels selected 

or modified as part of this work, shall be tested, verified, 

and submitted, using a computer system compatible with the 

latest version of the SWMM subniodels. 

The Scientific Atuhority is desirous that the SWMM system 
and extensions incorporated as part of this work be suitable 
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for execution on competitive large-scale computers of 
manufacturers other than the one used in this study. 
To this end, the Contractor shall submit to the Scientific 
Authority, three months before the end of this study, 
suitable terms of reference, including cost estimates 
for the conversion, testing, and verification of the 
SWMM system and extensions as incorporated in this study. 

Insofar as practical, the Contractor will consider the 
idea of program implementation on large-scale computers, 
such as IBM, CDC, UN I VAC, when developing new programs 
during this study. 

The cost of using computers for watershed modelling includes 
both the cost of machine time and the cost of engineering 
time. Interactive time-sharing via low-speed terminal 
will increase computer cost but may provide significant 
savings in engineering cost. The Contractor shall review 
and comment on the relative merits of batch processing with 
normal turn-around in the Toronto area and interactive 
time-sharing, considering the various phases of urban 

watershed modelling. 

I 

Guidelines shall be developed to allow users to estimate 
cpu time on at least one hardware configuration before 
running the SWMM. | 

All documentation of the new submodels shall be written 
using the format of the new User's Manual being compiled 
by the University of Florida. A master tape of the 
source programs, along with an appropriate Operator's 
Manual, shall be included with the final report. 

(b) Ver if i cat ion 

It is expected that each subsystem will be verified by 
comparing its performance to either existing Canadian data 
or to data available elsewhere, pending the approval by 
the Scientific Authority. 
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Where no appropriate verification data are available, 
a field program will be designed around the needs of this 
investigation and will be the subject of a separate terms 
of reference, requiring subsequent proposals. The 
Contractor shall identify such requirements in general 
terms within 30 days and on a continuing basis 
throughout the study. 

The Contractor shall make maximum use of the pertinent field 
data available in Canada or elsewhere. The Contractor's 
efforts shall be coordinated with other Canadian research 
activities (mostly research sponsored by DOE, MOE and NRC) 
and with U.S. Environmental Protection Agency Research 
Program. 

(c) 5WMM Model Parameters - Default Values 

There are numerous default values built into the SWMM to 
allow simulation to continue when the value of a model 
parameter is not known. These default values were often 
based on very limited data during the development of SWMM 
The Contractor shall review the default values and, where 
necessary, shall update them and relate them to Canadian 
conditions and practices. 

2.3 Structure of Report 

The report has been divided into three volumes, as follows: 

Volume I - Final Report 

Volume II - Technical Background 

Volume III - User's Manual 

Volume I gives a description of test applications, sensitivity analyses 
for existing SWMM routines, suggested modifications and simplifications, 
and a description of new programs. A technical reader will no doubt 
wish to read all of the volume. A general reader should be able to obtain 
a good appreciation of the study by reading only the Summary, Basic 
Concepts (Chapter 1), and the Conclusions and Recommendations (Chapter 13). 



29 



Volume M describes previous research and reviews critically 
the state-of-the-art for the benefit of more technically oriented readers. 

Volume Ml contains specific information regarding the computer 
programming aspects of the study for those who will be directly involved 
in simulations. This manual includes the EPA manual for SWMM, modified 
to include some of the recommendations that emanated from this study. 
It also describes how to use STORM and the other programs developed. 

Program listings for SWMM and STORM and for DAP and GSQM are 
presently available from the Water Resources Branch, Ontario Ministry 
of the Environment, ; 
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3. ASSESSMENT AND SELECTION OF STUDY AREAS 
3< 1 General 

Since the publication of a preliminary appraisal of storm 
water overflows by the U.S. Department of Health, Education and Welfare 
in I96A [25], a considerable number of monitoring and sampling programs 
have been instituted in the U.S. and Canada. Initially, these programs 
were oriented towards overall characterization of pollutant discharges, 
annual averages, and regional variations, rather than towards precise 
definition of the complex interrelationships determining urban runoff 
flow and quality from single events. However, the recent advances 
in computer model) ing of the urban runoff process have created the 
need for precise, short interval storm water quantity and quality data 
for single events, for the calibration, validation and refinement of 
these urban runoff models. Much of the existing body of storm water 
data is unsuitable for modelling purposes, and researchers have been 
forced to rely on limited data, or data from several sources. For 
example, Brandstetter [7] attempted to select three catchments with 
reliable rainfall, runoff and quality data on which to base rigorous 
comparisons between various models. Because of the limited availability 
of adequate data, the study was restricted to only two areas: Oakdale, 
a small 13"acre catchment with no quality data; and. Bloody Run, a 
large 2380-acre catchment, the data from which appeared somewhat unreliable. 

During the course of this present study, a number of promising 
collection programs were initiated to supplement this scarcity of reliable 
data and are currently underway. The University of Florida is conducting 
a program for collection, classification and dissemination of all current 
data in an effort to provide model users and developers with a comprehensive 
data base. An interim report [26] identifies '*2 potential data sources. 
However, only a few of these are at present available, and these are 
mainly the older published data extensively used during the SWMM development. 
A review of the more promising of these potential sites is contained 
in Volume II of this report. 

In the present study, efforts were directed first to a review 
of available U.S. data, and to a subsequent Investigation of Canadian 



studies which might be potential sources of modelling data. This 
chapter discusses the data requirements for urban runoff modelling 
purposes, and assesses the Canadian data which were used in this study. 

3.2 Data Requirements for a Study Watershed 

The SWMM requires an extensive amount of data to perform a 
detailed single event simulation. These data may be grouped into three 
categor ies: 

(a) Meteorological data describing the primary input to the 
modelling system, short interval precipitation data, 
and temperature data if snowmelt is being modelled. 

(b) Data describing the physical system, the nature and 
characteristics of the runoff producing area and its 
associated sewer system. These include the imperviousness 
ratio of the watershed, ground slopes, sewer system 
geometry, hydraulic roughness and other specific information. 

(c) Data for verification and calibration, consisting basically 
of reliable measurements of storm water flow and 

quality characteristics at specific points in the sewer 
system. These data may be compared with computed values 
to verify or adjust the urban runoff model for local 
cond i t ions. 

3.2.1 Meteorological data 

In assessing the applicability of meteorological data for urban 
runoff modelling, the critical considerations are timing, location of 
gauges, and network density. Analytical studies [35] have been performed 
using statistical precipitation synthesis techniques to determine optima! 
network design. However, all too often there are basic inadequacies 
in the available records. 

Because of the short response times and rapid flow fluctuations 
typical of urban storm water runoff, short time interval rainfall data 
are required for detailed modelling. In most cases a five minute 
interval is sufficient, but intervals as short as one minute may be 
necessary for modelling very small areas. The automatic tipping bucket 
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recording rain gauge is perhaps best suited for recording such short 
interval rainfall information. Each tip of the bucket (0.0) inches) 
is recorded on a moving chart, which may be read to five minute intervals. 
Chart speeds may be increased on some instruments to provide resolution 
at one minute intervals. Weighing type recording rain gauges are 
usually difficult to read at intervals of 15 minutes or less. In 
Canada, tipping bucket gauges are employed almost exclusively. 

Rainfall usually varies considerably from point to point within 
a watershed. A good location for a single rain gauge is in the centre 
of the study area, although for small areas, closely adjacent locations 
may be suitable. As the distance from the gauge to the study area increases, 
the error between the point rainfall measurement and the watershed 
average increases. This increase can be as great as 100? over a 
distance of only several miles [27]. For larger watersheds, it becomes 
even more important to have at least one rain gauge located in the 
study area, since the effects of distance and increased area combine 
to reduce point measurement accuracy. 

The density of the rain gauge network is also an important 
consideration as the catchment sizes increase. For small catchments 
(up to 1000 acres) a single rain gauge should provide reliable information. 
Some of the largest catchments modelled to date have been between 3000 
and 4000 acres. Significant differences in storm precipitation have 
often been recorded between several gauges located within watersheds of 
this size. The use of an individual rain gauge, as in these instances, 
may result in inaccurate simulations. Studies using dense rain gauge 
networks on small areas indicate that use of two rain gauges within the 
catchment decreases the error between point measured and basin average 
precipitation (as determined by many gauges) by as much as 50^ [27]- 
Other factors affecting the desired density of rain gauges would be the 
predominance of certain types of storm patterns over an area, and topographic 
influences such as major elevation differences. These may lead to 
significant errors between point measured and basin average rainfall. 
The desirability of using two rain gauges in an area cannot be overstressed. 
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Temperature data may be required for the modelling of snowmelt. 
Variations in temperature within an urban area are much smaller than 
variations in precipitation. Reliable information may be obtain«>d from 
the nearest meteorological station. For most snowmelt modelling, 
daily maximum and minimum temperatures are sufficient. 

In Canada, the source of most meteorological data is the 
Atmospheric Environment Service of Environment Canada. Types and 
formats of available data are discussed in Chapter 12. 

3.2.2 Physical system data 

Physical system data encompass all the geometric, hydrologic 
and hydraulic information necessary to describe the subcatchment areas, 
sewer system, and the basic connectivity of these elements, for the 
computation of storm water runoff flows. tn addition, information on 
street sweeping practices, land use, population and other parameters 
is also necessary for computation of runoff quality. Such a data base is 
common to most single-event models, and to a lesser extent to the more 
simplified continuous simulation models. Data are generally available 
from topographical and sewer maps, plans and profiles, aerial photos, 
field recona issance information, etc., through the local Municipal 
Engineering office. The physical data requirements for the SWMM are 
fairly extensive, but are relatively straightforward and may sometimes 
be reduced by adopting a more simplified modelling approach, as 
described in Chapter 10 of this report. A detailed summary of the 
physical data required for SWMM modelling is contained in the User's 
Manual {Volume III). 

The most frequent problems with physical data are missing, 
incomplete, or faulty information on basic documents such as plans and 
profiles. This often necessitates costly field verification of critical 
measurements. 

3.2.2 Verification and calibration 

Verification and calibration data provide the means to check the 
operation of the urban runoff model, or a part of It, either during the 
development of its basic routines, or during the application of the model 
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In a particular study. These data must, therefore, be as accurate as 
possible in order to play a useful role. Sufficient latitude is available 
in selection of certain of the model parameters in SWHM to achieve a 
closer agreement between the simulation and measured data for a particular 
watershed by calibration. Chapter k contains a discussion of the 
sensitivity of the SWMM to variations in model parameters. 

3.2.3-1 Flow measurement . The accuracy of flow measurements is 
largely dependent on the measuring device used. Methods applicable to 
sewer environments may be of three types: critical depth devices, tracer 
or acoustic techniques, or simple stage measurement coupled with an 
empirical flow formula. These techniques have been reviewed elsehwere 
by Wenzel [28], Schontzler [29], Marsalek [30], and Lager and Smith 
[31]. Only a brief discussion will be given here. 

Critical depth devices employ established stage discharge 
relationships coupled with a measurement of the stage during the flow 
event. Measurement of the stage may be accomplished with floats, gas bubblers, 
surface probes, capacitance probes, or ultrasonic devices. Accuracies 
of + 5% are possible with carefully calibrated critical depth sections. 

Tracer or acoustic techniques may be used to determine flow 
velocity directly. A knowledge of the flow area yields the flow rate. 
Accurate results may be obtained {^ S%) , although more sophisticated, 
expensive equipment is generally required. 

Stage measurement coupled with an empirical flow formula (such 
as Manning's equation) is the least accurate of all common methods. However, 
the simplicity and economy of this method has gained it wide acceptance 
for small area studies and for purposes unrelated to computer modelling. 
These techniques are susceptible to error because of the assumption of 
steady, uniform flow, which is rarely the case in most sewer systems. 
Discrepancies of up to 20-25% must be expected with this type of 
measurement system. 

A second basic consideration is the location of the measurement 
site. Aside from maintenance considerations such as ease of access, 
space requirements, etc., it is important to establish precisely which 
flow is being measured. Overflows, diversions, pumping stations and 
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storage facilities located upstream must be identified and their impact 
assessed regarding the measured values. The effects of backwater due 
to lakes, streams, or interception must also be evaluated. 

The length of the historical record for a particular site may 
be an important factor, especially when continuous simulation is contemplated. 
In this case, accuracy may be sacrificed in favour of the number of 
recorded events, which can be most useful in estimating annual discharges. 

3-2. 3-2 Qua! i ty samp I i ng . The collection of accurate storm water 
water quality data is perhaps the most difficult aspect of a monitoring 
program. An excellent assessment of currently available samplers was 
given by Shelley and Kirkpatrick [32]. Because no single sampler can 
carry out all the functions that may be required, the suitability of 
individual samplers should be assessed at the beginning of each program. 

Site selection is also of prime importance in establishing 
a sampling program. The location should be as close as possible to the 
flow measuring device. Factors that can affect the representativeness 
of the samples are: flow velocities at the sampling point; location 
in the cross-section at which sampling occurred; existence of a uniform 
conduit section upstream; contribution of tributary inflows; and, 
finally, the type of samples collected. 

Model calibration and verification usually necessitates frequent 
collection of discrete samples in order to determine the variation in 
const i tutents with progression of the storm flow hydrograph. Intervals 
of 15 minutes have been used in many collection programs. The duration 
of the sampling process will depend on the capacity of the sampler (if 
it is automatic), and the requirements of the modeller. It is often 
important to establish the characteristics of the start, peak and end 
of an event. Composite samples, made up of a series of smaller samples, 
yield time averaged values which may not be useful for direct model 
comparisons. However, comparison of mass discharges of pollutant 
recorded over one or several storm events with computer simulations may 
provide a good means of assessing model performance and overflow problems. 
Composite samples can sometimes be quite useful, and less expensive, 
than discrete samples, depending upon study objectives. 
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The number of pollutant constttuents analyzed affects both 
the cost and usefulness of the data. For calibration/verification with 
the SWMM results, data describing three basic constituents are required: 
BOO, suspended solids, and total conforms. These three pollutants can 
be modelled through in all phases of the program. 

The length of sampling record is an additional factor in data 
assessment, since the longer the record, the more data are available 
for establishing storm water conditions and model response. A longer 
sampling period can be especially useful when continuous simulation is 
planned. 

3- 2. 3- 3 Time synchronization . Of some importance in any collection 
program is the synchronization in time of the three measurements: 
rainfall, flow and sample. In the ideal situation, all three would 
be recorded on the same chart. This is possible with some of the newer 
recorders available, which have been used in recent data collection 
programs. In many small scale programs, the individual data are recorded 
separately. Where the instruments are also physically separated, there 
is a possibility of questionable timing of each recording. Apparent time 
shifts of 10 to 15 minutes between the rainfall and flow record are 
usually not a significant factor, and may be accounted for in the model 
calibration process by adjusting the starting time. It is usually 
advantageous when dealing with separate recording elements to employ 
a control clock and apply a check to the timing of each recorder during 
each si te visit. 

3. 3 Canadian Study Areas 

An Investigation of a number of on-going studies, municipal 
programs, and university research was carried out very early in this 
study to determine the availability of Canadian data suitable for modelling. 
Particular emphasis was placed upon Canadian data for several reasons. 
Firstly, conditions in Canada are well known to potential users and this 
study was primarily Canadian in administration and purpose. Secondly, 
all of the reliable American data available at the start of this study 
had already been used in numerous simulations and were currently being 
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employed In on-going American studies such as those of Brandstetter 
[7], Heaney et al [3'*1 . and others. Again, if suitable Canadian data 
did exist, their use would help focus attention in this country on the 
availability and requirement for further Canadian collection studies 
and the use of computer models in general. Finally, much of the 
American data, such as those recorded at Bloody Run, Cincinnati, and 
Atlanta, Georgia, were either incomplete or unreliable in one or more 
of the basic elements. Some of the more reliable U.S. data, such as 
Oakdale, Illinois, were from very small watersheds and therefore limited 
in scope. 

Data from eight Canadian watersheds were reviewed. In some 
cases, the data represented purposeful efforts to collect data for 
modelling, whereas in other cases the collection program was oriented 
towards overall assessment and monitoring studies conducted by municipal 
departments. The areas reviewed are summarized in Table 2. 

TABLE 2. SUMMARY OF AVAILABLE CANADIAN DATA 



Area 


Size 

(acres) 


Type 


Rainfal 1 
Data 


Flow 
Data 


Qual i ty 
Data 


Period of 
Record 


Brucewood 


48 


Storm 




G 


G 


G 


'74- '75 


East York 


40 


Storm 




G 


F 


P 


'74-'75 


Malvern 


58 


Storm 




G 


G 


G 


'73-'75 


Kingston 


89 


Storm 




G 


G 


- 


'73-*75 


Hal ifax 


168 


Comb i ned 




G 


G 


G 


'71-'74 


Ma 1 ton 
Airport 


495 


Storm 




1 

a 


G 


G 


'74-'75 


Bannatyne 


542 


Combined 




p 


G 


F 


'69-'71 


West 
Toronto 


2330 


Combined 




G 


F 


— 


'69- '76 


Montreal 


2880 


Combined 




P 


F 


f 


'73 


G = Good 
















F = Fair 
















P = Poor 
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3 . 3 • 1 Brucewood catchment, North York [39] 

The Brucewood site provides a good source of modeiJIng data. 
Many low intensity, and several medium intensity storms have been monitored 
for both flow and quality. Rainfall data are provided by an automatic 
tipping bucket rain gauge located acjacent to the study area. Flow 
data are obtained using a weir at the outlet, providing continuous flow 
measurement of all events. Samples are collected automatically at 
four minute intervals from the outlet pipe during storm flows. In addition, 
snowmelt studies and investigations of surface pollutant accumulations 
have been carried out on the test area. 

Rainfall and flow information is remotely recorded on synchronized 
charts, making this a good source of urban runoff data. 

3.3.2 Barrington catchment. East York [h]] 

A limited amount of data have been collected from this area. 
Rainfall data were provided by an automatic gauge within the area. 
A flow measuring weir was used to measure the runoff hydrographs, but 
examination of these suggests that, in many cases, readings were incomplete. 
Most storm events were of very low intensity and, therefore, not particu- 
larly useful for detailed modelling. In some cases, large unexplained 
timing shifts occurred between recorded rainfall and runoff. Both 
manual and automatic sampling were employed in collecting these data. 
The collection times of the bacteriological samples were not shown on 
the lab reports and, consequently, these are of no use for model verifica- 
t ion. 

These data were considered to be of limited use for modelling 
purposes. 

3.3.3 Malvern [46] 

Very good runoff data have been collected on this area. 
Rainfall was monitored within the catchment, and recorded on the same 
chart as the runoff. Flow measurement was accomplished using a calibrated 
weir located at the drainage outfall. For the period of l973-'74, 
about 22 storm were monitored and used for runoff modelling. The runoff 
quality monitoring started in 197^- The project will be continued in 
1976. 
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3.3.'* Calvin Park, Kingston [8] 

Very good runoff data have been collected on this area. 
Rainfall was recorded within the catchment, on the same chart as the 
runoff. Flow measurement was accomplished using a weir. Only a 
limited number of high intensity storms were available. No quality 
data were obtained during this monitoring program, but it is understood 
that recently the study has been expanded to include water quality 
measurement. I 

The quantity data were considered very good for modelling 
purposes. 

3.3.5 Halifax [h2, k}, hh, ^5] 

Good data have been collected in the Halifax area in two programs 
the first, in 1970-71, monitored combined storm flows; and the second, 
in 197^, recorded winter runoff. Rainfall data are excellent, supplied 
by an automatic recording rain gauge located within the area, supported 
by several standard rain gauges. Flow was recorded using a weir at the 
inlet to an overflow retention tank. Water quality was sampled both 
automatically and manually at the outlet. Most of the recorded storms 
are of very low intensity, and only a few high intensity storms are 
ava i lable. 

These data are considered suitable for modelling purposes. 
However, other independent studies were already using the best data 
from this area, and repetition of these efforts was considered 
unnecessary. 

3.3.6 Bannatyne, Winnipeg [^O] 

A large amount of data has been collected from several urban 
catchments in Winnipeg. The most complete data, from the Bannatyne 
combined sewer area, were collected from 1969-1971- Originally, 
recording rain gauges located within the study area were available, 
but these records were misplaced by the collecting agency. Data from 
three meteorological stations located outside the study area are avail- 
able. Flows were measured using a weir, and flows for a large number 
of complete storm events were recorded. However, a limited number of 
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events for which both rainfall and flow data are available have been 
identified. An automatic sampler was used to collect water quality 
samples, and data are available describing quality during part of each 
of the events. These data were recorded separately from the flow 
charts with some unexplained time shifts occurring. 

Some of these data were considered to be suitable for 
modelling. The weakest element o'^ the collected data was the rainfall, 
but reasonable agreement between the nearby rain gauges was evident 
for several of the recorded events. Quality data were incomplete in 
some events due to the short operating time of the sampler used, but 
in view of the state-of-the-art of quality modelling, these data 
were also considered to be of some use. 

3. 3-7 Toronto International Airport [36] 

Storm water monitoring was carried out for several areas on 
the Malton Airport site, during 197A and 1975, as part of an assessment 
of environmental problems. One of the larger catchments was an area 
of ^95 acres consisting of level grassy areas (521) and of flat impervious 
surfaces such as runways and maintenance buildings (^48^). 

Flow levels were measured using a Manning "Dipper" level 
recorder located in a 72-inch outfall sewer. Flow was determined using 
Manning's formula. An automatic recording rain gauge was located 
nearby at the airport meteorological station. Data from several storms 
are available for this catchment end were considered to be suitable 
for modelling purposes {subject to the limitations of the flow recording 
technique used). Quality sampling was conducted using automatic samplers 
at the outfall to this area. 

3.3.8 West Toronto [37] 

The West Toronto watershed has been monitored for overflow 
quantity for the last six years. The data collected for this large area 
were of interest both for detailed single-event modelling and for 
continuous modelling. Very good rainfall data are available for this 
area from a recording rain gauge operated by the Atmospheric Environment 
Service (AES), located approximately in the centre of the area. The area 
is large enough to warrant additional rainfall information; fortunately 
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correlation with another AES rain qauqe located outside the 
boundary of the area was possible. The area drains to a large outfall 
sewer in which water levels have been recorded during all srorm events. 
An examination revealed that no backwater conditions affected the level 
measurements. However, a complication was a series of diversion weirs 
at the outlet designed to intercept dry weather flow in sewers leading 
to the outfall. However, it was possible to reconstruct an approximate 
outflow hydrograph by computing uniform flow from the level measurements, 
and adding an estimated component representing flows diverted both at 
the outlet and at a number of minor points along the upstream boundary. 
Many major overflows have been recorded at the outfall pipe (108" dia.) 
where the capacity of the intercepting sewers forms only a small propor- 
tion of the outflow. While flow measurements were recorded in the major 
intercepting sewer during the period in which overflows were recorded, 
recent measurements are available to more precisely establish the 
intercepted flow. No quality measurements were available from this 
area. 

Overall, the quantity data from West Toronto were Judged suitable 
for modelling purposes. The inherent inaccuracies in the flow determina- 
tion are noted, but there are important aspects to the Toronto data. 
It is a large area (2330 acres), served by combined sewers and prone 
to surcharge. Few such areas, which are typical of many large city 
areas, have been modelled using the SWMM. The length and continuity 
of the record made it valuable for continuous simulation, in which case 
point accuracy of individual peak flows was of less concern than the 
frequency, duration and total volume of overflow. The area had been 
previously modelled in detail with the Dorsch Hydrograph Volume Method 
(HVM) and, consequently, a considerable amount of physical data was 
available in reduced from. This also presented an opportunity to 
perform an approximate comparison between these two models. 

3-3-9 Papineau-Curotte catchment, Montreal [38] 

Limited measurements were also available for a large area in 
Montreal. Rainfall data were available, but none of the recording 
rain gauges were located within the area. Flow was measured at the 
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outlet using a weir, and data were available for ?0 events of relatively 

low intensity rainfall. Quality data were available from automatic 

sampling of only four of these events. It is not clear exactly how the 

flow measurements were obtained since this area also drains to an interceptor/ 

overflow system. The quality sampling was conducted some distance 

downstream of the study area, and so the contributing area is uncertain. 

The data from this area were considered to be unsuitable for 
the present modelling purposes. 

3.3-10 Current collection programs [33] 

Under the sponsorship of Environment Canada, two urban areas 
have been instrumented for collection of complete urban runoff data. 
This program is specifically designed to provide data for runoff 
mode 1 1 i ng . 

The first of these is a 176 acre area, located in Hamilton, 
served by a combined sewer system. This area is mainly single-family 
residential, and is about k^Z impeivicus. Rainfall is being measured 
with two tipping bucket rain gauges, located within the watershed. 
A modified v-notch weir and a bubble type level sensor are used for 
flow measurement. Automatic quality sampling is being carried out. 
The data from these elements are being logged simultaneously on a single 
recorder device. In addition, data on road and catch basin cleaning, 
snow removal and de-icing, and temperature are being collected. 

A second test area, located in Toronto, has also been instru- 
mented for urban runoff data. This area covers 338 acres of predominantly 
single-family residential land use. The area is >i3t impervious. 
Precipitation measurement is by a single tipping bucket rain gauge and 
a special high resolution recorder. Flow measurement apparatus consist 
of a bubbler type level recorder and a partial trapezoidal weir. 
Automatic sampling is done during storm events, initiated by a 
predetermined water level. Continuous temperature measurement is carried 
out at the outfall site, and information on municipal cleaning practice 
is also being collected. 

Both of these measurement programs have been in operation 
since the late summer of 1975. These are potentially excellent sources 
of runoff data, specif icaMy oriented towards runoff modelling. 



3.3.11 Summary of areas reviewed 

Several sampling programs were identified as being potentially 
suitable for modelling purposes within the context of the present study. 
Data were being collected from the Brucewood area which would be available 
for inclusion in this study. The only large areas with acceptable data 
were in Winnipeg and West Toronto. It was considered that the less 
accurate Toronto data would be very useful in portions of the study that 
did not include detailed model flow verification, such as simplified 
simulation and long term simulation. The Winnipeg data represented 
a medium-sized watershed with good records of quantity and quality, 
which would be useful in various parts of this study. Finally, 
limited data from the Mai ton Airport and Barrington watersheds were 
considered to be suitable for some study purposes, and for general 
quantity simulation. 

A review of data available from U.S. watersheds is contained 
in Volume II. Data from the Lancaster area and the Atlanta area 
were obtained for use in this study. Although these data were considered 
the best available in the U.S. at the start of this study, there were 
some inadequacies. Lancaster flow measurements were considered somewhat 
unreliable, and the absence of good rainfall data for the Atlanta area 
had resulted in poor simulations in previous work. 

3. k ConcI usions 

(a) Based upon a review of available data, it was found that 
complete sets of measured storm runoff data were scarce 
both in Canada and the U.S.A. Most of the acceptable 
data sets which have been published are for very small 
catchments, and have been modelled many times in pre- 
vious studies. 

(b) Several sets of data from Canadian urban watersheds were 
identified as being suitable, at least in part, for 
modelling within the terms of the study. 

(c) There is a continuing need for more runoff data for 
Canadian watersheds from programs specifically designed 
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for model linq needs. Existing rnotiitorinq programs should 
be continued, 
(d) The measurement of rainfall is often overlooked in setting 
up monitoring sites. In general, at least one automatic 
recording rain gauge should be located In or immediately 
adjacent to the monitored watershed. For large catch- 
ments, greater than 1000 acres, two rain gauges should be 
used . 



k. SWMM FLOW SIMULATION FOR SELECTED CANADIAN WATERSHEDS 

k. I General 

Accurate flow simulation is a prerequisite for most studies 
of storm drainage and poHution abatement. The first comparisons of 
flow simulations given by different urban runoff models were carried 
out by those involved in the construction of the models and their scope 
was mainly confined to calibration and verification, based on only a 
few events. Sixteen different flow simulation models are discussed 
in a recent study conducted by Battel le Pacific Northwest Laboratories 

[7]. 

A study comparing different flow simulation models was sponso- 
red by Canada Centre for Inland Waters [i3]- This study investigated 
a much larger number of storm events than other previous model comparisons, 
and the models were evaluated from the standpoint of a potential user 
rather than on a strictly theoretical basis. Small urban catchments 
(less than 100 acres) were used for these comparisons. The urban 
runoff models studies in detail were limited to nonproprietary, 
single-event, runoff simulation models. The following runoff models 
were selected on the grounds of availability, technical sophistication, 
and recent successful application: 

(a) Road Research Laboratory Model (RRLM) 

(b) Storm Water Management Model (SWMM) 

(c) University of Cincinnati Urban Runoff Model (UCURM) 

In addition, a limited number of simulations was carried 
out using the Dorsch HVM (Hydrograph Volume Method) model and the 
Unit Pulse model developed at (Jueen's University. 

The main conclusion of the study was that, in general, flows 
simulated by all urban runoff models agreed fairly well with measured 
runoff hydrographs for small urban catchments. A statistical comparison 
of the simulated hydrographs with measured flows resulted in good 
agreement for SWMM, fair to good for RRLM and fair for UCURM. These 
results are summarized in Table 3- Of the nonproprietary runoff models 
investigated, SWMM gave the best results. Since it is, in general. 
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a more versatile and better documented model, the SWMM was recommended 
for subsequent work on larger areas. 

In large urban areas, routing in the trunk sewer system plays 
a significant role in determining the outflow hydroqraph. In the 
framework of the present study, the SWMM quantity model was assessed 
by the simulation of storm flows on large watersheds by both the 
RUNOFF and TRANSPORT blocks, and results were compared with measurements 
and those of other more sophisticated routing models. 

TABLE 3. A COMPARISON OF RUNOFF MODELS [ 1 3] 



{COMPUTED PEAK FLOW)/ 
(MEASURED PEAK FLOW) 



(COMPUTED TIME TO 
PEAK FLOW) /(MEAS- 
URED TIME TO PEAK 



(COMPUTED RUNOFF 
VOLUME) /(MEASURED 
RUNOFF VOLUME) 



MODEL 



No. of 
Obser- 
vat ions 



Aver- 
age 



Stan- 
dard 


No. of 
Obser- 


Aver- 


Stan- 
dard 


No. of 
Obser- 


Devia- 


vat ions 


age 


Devia- 


vations 


t ion 






t ion 





. Stan- 
Aver- , , 
dard 



age 



Devia- 
tion 



RRL MODEL 
{Non-Cal i brat- 
ed) 



57 1.09 0.25 



57 1.03 0.23 20 



1.09 0.27 



SWM MODEL 
(Non-Ca 1 i brat - 
ed) 



58 1.13 0.23 



58 0.95 0.22 20 



0-96 0.22 



UCUR MODEL 
(Non-Cal ibrat- 
ed) 



27 1.15 0.23 



27 0.86 0.\k 10 0.83 O.U 



SWM MODEL 
(Cal ibrated) 



31 1.01 0.21 



31 0.96 0.11 10 



1.03 0.17 
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Basic Concepts 



The urban watershed drainage system (Figure h) may be considered 
as two distinct systems for computer modelling purposes: 

(a) The basic subcatchment system in wbirh the excess rainfall 
is transformed to overland flow. The overland flow hydro- 
graphs at the drainage manholes form the inlet hydrographs 
for the transport system. 

(b) The transport system in which the flows are routed through 
the trunk sewer networks 



it. 2.1 



The basic subcatchment system 



Each subcatchment is conceptualized as a flow plane over which 
overland flow occurs. Representative values of imperv iousness , average 
ground slope, infiltration, retention storage, and land use are assigned 
to each subcatchment. The rainfall depth over each subcatchment is the 
basic input to the overland flow or RUNOFF model. The model sequentially 
accounts for infiltration and retention depth to determine a net excess 
rainfall depth on the overland flow plane. The excess rain depth is used 
to compute an overland flow rate for each time step. The SWMM uses a 
very simple algorithm for this runoff computation, as shown in Figure 5 i. 
Subca tchments can be used to represent large portions of the entire 
watershed, or individual drainage areas, depending upon the desired degree 
of detail. The small local sewer pipes and gutters are usually ignored 
in large subcatchment s . The SWMM also provides the facility to account 
for simplified routing in local pipes in the RUNOFF block. Instead of 
considering the subcatchment overland flow hydrograph to flow directly 
to the inlet manhole, it can be routed initially through a "GUTTER" 
subroutine which is used to account for the many ditches, lateral pipes 
and street gutters, which are not usually included in the TRANSPORT block 
computat ionSi 



k.1.2 



Routing models 



A routing model combines the runoff flows from several subcatchments 
and routes these flows through the main sewer network. Several urban 
runoff models of varying levels of sophistication have been developed for 
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the simulation of this process. Such TRANSPORT models range from those 
using simple time-offset methods to those employing complex numerical 
schemes for the solution of the dynamic wave equations. While the latter 
approach is applied by some one-event simulation models (WRE model, 
Dorsch HVM) , long term simulation models generally use simpler techniques 
and others do not consider sewer routing at all (e.g. STORM). A compre- 
hensive TRANSPORT model should be able to account for: 

=''(a) the unsteady and nonuniform nature of flow through 

sewers; 
=''(b) the storage capacity of the sewer system; 
(c) the dynamic effects of the backwater condition and 

the energy losses along the sewers and at manholes; 
"(d) control elements including manholes, weirs, orifices, 

pump stations, storage units and gates; 
"(e) the distribution of flow through branched and looped 

systems; 
(f) flow reversal and adverse sewer grades; and, 
-■■(g) the variety of conduit shapes conmonly encountered in 

engineering practice. 

In order to assess the TRANSPORT blocl< of the SWMM, the response of the 
model was compared to that of two other models under varying physical 
and hydrological conditions. The routing models considered in this phase 
of the study are described below. The results of comparing these different 
models are described in Section k.k. 

4.2.2.1 SWMM TRANSPORT block . The TRANSPORT block of the SWMM [^7] 
routes the storm water flows through a converging branch sewer network. 
Dry weather flows and infiltration into the sewer system are computed. 
Two types of flow diversion structures, two storage basins and one lift 
pumping station can be modelled. The model is limited to the simulation 
of single storm events and to a maximum of 160 hydraulic elements. Larger 
systems may be handled by modelling the major areas separately. The flow 
routing is based on the quasi-steady dynamic wave approximation, this being 



Accounted for in SWMM 
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a form of the St. Venant equations (the momentum equation for a dynamic 
wave), with the omission of the local acceleration term (equation 1). 
While an implicit scheme is used to solve for continuity (equation 2), 
an explicit scheme has been adopted for the solution of the momentum 
equation. 



S^ (1) 

(2) 



3h V av 
3x g ax 


= s 

o 


3Q A 

ax t 


= 



where: h = depth 



X = longitudinal coordinate along 
channel bottom 

t = time 

V = velocity 

g = gravitational acceleration 

Q = flow rate 

S = channel bottom slope 
o 

S = friction s lope 

9 
A = cross-sectional area of flow 



A Newton-Raphson technique is applied for the solution of the nonlinear 
continuity equation. The friction slope is evaluated from Manning's 
formula and lateral flow is not simulated. The solution procedure 
basically follows a kinematic wave approach, in which disturbances are 
allowed to propagate only in the downstream direction. As a consequence, 
backwater effects are not modelled beyond the domain of a single conduit, 
and downstream conditions will not affect upstream elements in the normal 
TRANSPORT cCMnputat ions . However, backwater can be approximated at a 
maximum of two locations by specifying a storage element and providing 
appropriate geometric input data. These effects are, again, not transmitted 
to the upstream elements. The program does not simulate pressure flow 
conditions. Surcharging is modelled simply, with flows in excess of the full 
flow conduit capacity stored at the upstream manhole until capacity exists 
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to accept the stored volume. Only treeshaped systems are permitted. 
Accordingly, interconnections, looped networks, and flow reversal cannot 
be simulated. A maximum of three inlet and one outlet conduits can meet 
at any manhole. The geometries of 12 commonly used closed conduit 
shapes are programmed internally, and three additional arbitrary shapes 
may be specified by the user. 

^4.2.2.2 Oorsch Hydrograph Volume Method (HVM) . The HVM model [7, 48, 
49, 50] simulates the system as an interdependent network and the effect 
of conditions in each element (backwater effects, surcharging and 
pressure flow conditions) are modelled. The model accounts for combined 
and storm sewer systems, consisting of loops and converging and diverging 
branches and open channel networks. The flow routing is based on an 
implicit finite difference solution of the energy and continuity equations 
(3, 4), which are applied at each node and along each sewer segment. 

r^ + ~ (5^) ^ - It- ^ 4 (V - V cos=.) = s - s, (3) 

3x ax 2g g 9t gA q of 

- dx + — dx = qdx (4) 

where: h = depth 

x = longitudinal coordinate along channel bottom 

t = t i me 

V = velocity 

g = gravitational acceleration 

q = lateral inflow 

• = junction interior angle 

S = channel bottom slope 
o 

S, = friction slope 

A = cross-sectional area of the flow 

Q = flow rate 

An interative solution technique is used, applying a repetitive 
downstream and upstream computational scheme. A maximum of three inflow 
and two outflow conduits can come together at a node. A maximum of 1100 
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sewer segments can be accommodated in one calculation area. The model 
simulates weirs, orifices, gates, diversion structures, retention and 
drop structures. The surface runoff inflow enters the conduit elements 
laterally from the surface runoff model. The model, like SWMM, also 
simulates dry weather flow. This proprietary model has been applied 
extensively. 

^.2.2.3 WRE Model TRANSPORT block . The WRE model [51, 52] is a modified 
version of the SWMM. The TRANSPORT block has been improved, and the 
new model has additional capabilities. The dynamic wave equations are 
solved by a special explicit finite difference formulation. The 
continuity equation (5) is applied at the nodes (sewer junctions and 
connections) and the momentum equation (6) is applied along the links 
(sewer and channel reaches). 



fill = M 
^3t^t As^ 



3i = _„ AS + 2V ^ 

3t ^ f 3t 



w2 3A „ W 



(5) 
(6) 



where: ZQ_ = net inflow to the junction 

Q = flow rate in the conduit 

V = velocity 

■ 
A = cross-sectional area of the flow 

A = surface area of the node (manhole) 

H = hydraul ic head i 

S, = friction slope 

t = t i me 

The sewer system is simulated as an interdependent network and the effect 
of each element on the system is modelled. This is achieved by applying 
a Newton-Raphson iterative technique in solving the whole system as one 
unit at each time step. Surcharging conditions can be modelled at the 
surcharged junctions. An approximate first-order correction based on 
the Hardy Cross method is applied independently at the surcharged 
junctions. When the hydraulic grade line intersects the ground surface, 
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the excess flow volume is lost from the system. A maximijm of eight conduits 
can join at a node and up to 320 pipes are modelled. The model allows 
for flow reversal and, consequently, pipes are not specified as inlets 
or outlets to a manhole. Looping effects and converging and diverging 
branches can also be simulated. The model allows for the geometries of 
five closed conduit shapes and a trapezoidal open channel. Four hydraulic 
control devices are allowed, namely, weirs, gates, orifices and pumping 
stations. As in the SWMM, all flow enters at junctions and no lateral 
flow is considered. 

if,} SWMM Sensitivity Analysis for Selected Input Parameters 

An appreciation of the sensitivity of the SWMM computations 
to changes in input data is important in setting up and conducting detailed 
simulations. This is especially true when some of the data are unavailable 
or unreliable, and a decision must be made to proceed with a simulation 
or obtain more basic data. Because of the large data requirements, 
additional collection may be a costly and time-consuming process (such 
as obtaining additional topographic mapping). A limited sensitivity 
analysis was conducted on the main RUNOFF and TRANSPORT block parameters 
as a part of this study. The purpose was to increase familiarity with 
the model, and to establish the most sensitive parameters before 
conducting detailed simulations of the West Toronto and Winnipeg areas. 

it. 3.1 RUNOFF blocl< sensitivity 

A review of the basic computational theory upon which the 
RUNOFF block is based indicated the following parameters would directly 
affect the surface runoff from a given area: 

- imperviousness; 

- width of overland flow; 

- infiltration; 

- retention depth; 

- manning ' s ' n ' ; and , 

- ground slope. 
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Simulations were performed in which these parameters were varied over 
practical ranges and the results compared with those obtained using the 
default values built into the SWMM. A hypothetical watershed consisting 
of 25 small 11 acre subcatchments was used in these sensitivity simulations. 
A one hour triangular rainfall pattern with a peak intensity of four 
inches per hour and a total volume of two inches was applied. By 
varying parameters on each of the identical subcatchments, numerous 
sensitivity analyses could be performed in one program run. These 
sensitivity simulations were all run for a duration of 200 minutes, which 
is adequate to simulate the runoff hydrograph from a 60 minute rainfall 
event. At the end of simulation, there is usually some amount of surface 
water storage left on the catchment which would eventually run off. 
Changing any parameter which delays runoff from the subcatchment , such 
as ground slope 'n' or "width", results in a greater amount of surface 
storage remaining at the end of simulation and, therefore, affects the 
volume of runoff during the simulation. This left-over storage would 
eventually run off, however, if the simulation were conducted for an 
indefinite period. The following discussions consider the actual period 
of simulation, and for this reason parameters which do not affect water 
losses, such as 'n', will have an effect on runoff volumes. 

^4.3. I. I Sensitivity to Jnf i I trat ion rates . Figure 6 shows the result of 
the variation of infiltration parameters on the peak flow and volume of 
surface runoff. Different rates of maximum and minimum infiltration were 
investigated, all using the default value for the decay rate of inTiltration 
with time. The rates were varied in multiples of the default values 
(3.0"/hr maximum, 0.52"/hr minimum). The percent imperviousness for 
all three cases was ^0%. i 

The results indicate that the peak flow is not very sensitive 
to the infiltration rates, since only an 8Z variation in peak flow was 
observed. A somewhat greater effect was observed on total runoff volume, 
which varied some 25% over the range of rates applied. For areas with 
a greater amount of pervious surface, the peak flow and volume would be 
more sensitive to infiltration, and the converse would be true for highly 
impervious areas. For lower intensity storms, during which all of the 
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I 
rainfall infiltrates into the pervious area, the runoff would be less 
sensitive to infiltration rates, unless these were reduced enough to 
allow an excess rainfall to develop on the pervious areas. 

if. 3. 1.2 Sensitivity to retention depth . Surface retention depths are 
specified for the pervious and impervious areas modelled. The default 
values are 0.062" and els'*", respectively. Four sets of retention 
depths were simulated, ranging from to ^ times the default values for 
imperviousness and pervious areas. 

The results are shown on Figure 6. The sensitivity observed 
was similar to that for the infiltration. Peak flow was relatively 
insensitive with a variation of 9^, and volume of runoff only moderately 
sensitive with a total variation of 20^. The similarity to the infiltration 
results was anticipated, since both parameters act in the same way to 
modify the rainfall excess depth over the subcatchment . 

'♦.3-1.3 Sensitivity to ground slope . Figure 6 shows the results of 
varying the subcatchment ground slope from 0.001 to 0.10. At higher 
slopes (greater than ]%) , there appears to be a levelling off in the 
curves, indicating reduced sensitivity. For flatter slopes, there 
is considerable sensitivity of peak flow to the ground slope. Runoff 
volume also shows greatest sensitivity at lower slopes, but this is 
attributed to a greater buildup of water on the subcatchment surface 
which would drain off slowly If the simulation were continued for a 
longer time. As the ground slope is increased runoff occurs faster, with 
the result that a greater runoff volume has occurred by the end of the 
simulation. If the simulation were run for an indefinite period, the 
volume of runoff would be the same for all slopes. 

it.3-1-^ Sensitivity to Manning's 'n' . Manning's 'n' values for pervious 
and impervious surfaces were varied in multiples of the default values 
as shown on Figure 6. Both volume and peak flow increase as 'n' 
is reduced, and the sensitivity is greatest for low values. This is 
explained by the fact that lower 'n' values permit faster runoff from 
the subcatchment and, consequently, less surface storage at the end 
of the simulation. For much longer simulations, runoff volume 
would be independent of 'n'. 
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A. 3. 1-5 Sens i t i vl ty to _ jmpe ry lousness . The percent Fmperv lousness of 
a subcatchment has an almost linear effect on computed peak outflow and 
volume. The outflow increases as the proportion of impervious area 
increases. For very low values of imperv iousness , the outflow would 
be primarily dependent on pervious area runoff alone. Figure 6 shows 
the variation of outflow rate and volume with imperviousness. Graham 
et al [53] also found this to be the most influential surface runoff 
parameter , 

^.3.1.6 Sensitivity to width of overland flow . The width of overland 
flow was varied from 5^ of its correct value (1^1^ ft for the hypothetical 
catchment), to itOO% of this value. This parameter appears to be less 
sensitive at higher values than at low ones. Both peak flow and volume 
show similar effects. The width is actually a measure of the ease with 
which surface runoff occurs. Because the continuity equation is solved 
on the surface of the subcatchment in terms of rainfall, surface storage 
and outflow, the resulting sensitivity is nonlinear. Thus, if the width 
is reduced by one-half, the outflow rate for a given depth is also reduced, 
but more excess water builds up on the catchment surface which tends to 
increase the outflow. This effect (further discussed in Chapter 9) is 
employed to provide compensating storage when conduits are neglected in 
simplified simulations. The outflow volume is reduced as the width 
is decreased, because it requires a greater depth of excess water on 
the flow plane to discharge the same quantity of runoff. Therefore, 
at the end of simulation, a greater excess water depth is left on the 
surface. 

^4.3.2 TRANSPORT block sensitivity 

A series of conduits and a single II acre subcatchment were 
used for the sensitivity analysis for the TRANSPORT block. Several 
variations of conduit lengths, number of manholes and different total 
conduit lengths were investigated. This test configuration is illustrated 
on Figure 9. 
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A. 3. 2.1 Sensitivity to conduit length . The same inlet hydrograph was 
routed through five conduits whose lengths varied from 500' to 10,000'. 
This was a test of the consistency of the TRANSPORT model, and the 
results are shown on Figure 7- There was a logical increase in routing 
effect, attenuation of flow peak and a delay in the time to peak, as 
the conduit length increased, 

4.3.2.2 Sensitivity to the number of conduits in a given length . 
Simulations were made in which an identical inlet hydrograph was routed 
through a fixed length, made up of varying numbers of conduits and 
manholes. For example, the hydrograph was routed through total lengths 
of 1000', 2000' and '+000' of pipe, alternatively composed of four, two 
and one conduit sections separated by manholes. The results are shown 
on Figure 8. 

Figure 8A shows little difference if the hydrograph is routed 
through four sections of 250' or a single section of 1000'. Both the 
peak flow and volume are almost identical. On Figure 8B, a slight increase 
in peak flow is observed as more conduits were used in the routing. 
The peak flow is approximately l6% less with a single conduit of 4000' 
than with four 1000' pipes. These results tend to support the maximum 
conduit length of AOOO'-SOOO' recommended in the SWMM User's Manual. 
Use of lengths greater than this will result in less accurate routing 
computations due to numerical instability of the finite difference 
technique used in the solution of equations (I) and (2). 

4.3.2.3 Effect of surcharge . In order to investigate the effect of 
surcharge, a simulation was conducted with the central pipe in a five- 
conduit sequence deliberately reduced in size to produce surcharging. 
The results are shown on Figure 9. Conduits 30? and 308 upstream of 
the surcharged section produced identical hydrographs to those produced 
when no surcharging was modelled. The hydrograph in conduit 309, which 
surcharged, has been truncated at the pipe capacity, and the excess 
volume stored at the upstream manhole for later release. The 
hydrographs in conduits 310 and 3II are similar to that in the 
surcharged pipe, but the effect is delayed according to the travel time 
to each one. These results are consistent with the routing procedure 
employed in the SWMM. 
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h.'i-^.h Sensitivity to pipe slope . Simulations were made in which 
the overland flow hydrograph was routed through 1000 ft of conduit with 
slopes varied sequentially from 0.15; to 10?;. Very little sensitivity 
was observed in the peak flow and volume of runoff after the routing. 
The hundredfold change in slope produced only a 3% difference in peal< 
flow, and a 0.2? change in volume, as shown in the following table. 
This was a limited case, however, in which conduit routing effects were 
small in relation to overland flow routing effects because the pipe used 
was short. A greater sensitivity would be expected where longer pipes 
were used. 

TABLE k. SENSITIVITY OF SWMM TRANSPORT BLOCK TO CONDUIT SLOPE 



Inlet Hydrograph 
after 0. ]% condui t 
after 0-5^ conduit 
after ]% condui t 
after \0% conduit 



Peak Flow 
in cfs 

10.269 
9.932 
10.087 
10.166 
10.238 



Runoff Volume 
in ft^ 

23256.9 
23303.7 
23258.1 
23255.1 
23256.3 



't.3.2.5 Sensitivity to Manning's 'n' . Variations in Manning's 'n' 
produced more significant changes in the resulting hydrograph than did 
either pipe slope or conduit length. Roughness values from 0.013 
(the default value) to 0.052 were used in this simulation, and the 
results are shown on Figure 10. As the roughness coefficient is increased, 
the routed hydrograph is retarded and attenuated. The total volumes 
routed through the conduits were virtually identical for all values of 
'n'. The surface hydrograph is also shown in Figure 10 for comparison,, 



it. 3. 3 



Conclusions and sensitivity analysis 



The percentage of imperviousness had the most significant effect 
on surface runoff of all of the parameters tested. Both the peak flow 
and volume of runoff were almost directly proportional to the amount of 
impervious area. The remaining parameters are ranked in order of their 
decreasing influence: 
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- width of overland flow; 

- infiltration rates; 

- retention depths; 

- roughness coefficients; and, 

- ground slope. 

The SWMM default values for these parameters, which are shown on Figure 6, 
lie on the relatively flat portions of the sensitivity curves, and appear 
to be reasonable estimates when no actual site measurements exist. 

The maximum length of a single conduit modelled in TRANSPORT 
should be limited to ^000 to 5000 feet. Use of longer conduits will 
result in inaccuracies in the routing computations. Subject to 
this limitation, the number of conduits used to model a particular uniform 
section of an actual pipe appears to be fairly flexible. This analysis 
indicates that, over a single reach, the outflow hydrograph is not very 
sensitive to changes in conduit slope, as long as the conduit does not 
become surcharged. 

The roughness coefficient used in the TRANSPORT routing can 
have a significant effect on the outflow hydrograph and, within the limits 
of hydraulic experience, should be considered to be a calibration 
parameter. 

k.k SWMM Flow Simulation 



Simulations were performed on several test areas for which 
adequate storm flow measurements were available. All simulations were 
conducted in a detailed fashion, with each catchment being subdivided 
into a large number of subcatchments , and the major sewer elements were 
modelled in the TRANSPORT block. The default values were used for retention, 
infiltration and surface roughness. 

't.'l.l Bannatyne, Winnipeg 

This 5't2-acre combined sewer and area was subdivided into '•I 
subcatchments (Figure 11) for detailed simulation. As previously noted, 
rainfall data were not available for the area itself, but were obtained 
from three surrounding gauges which were a considerable distance from 
the site. A weighted average procedure based upon distance to the area 
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was applied to the records of these gauges for each event considered. 
Six medium to iarge storm events, for which reasonable agreement existed 
between gauge records, were selected for simulation. It was expected 
that the absence of rainfall data measured in the watershed would be 
a source of discrepancy between computed and measured hydrographs. 
However, the results of these simulations, discussed in the following 
paragraphs, show a reasonable agreement. 

Storm of June 19, 1971 

In this simulation (Figure 12), the computed hydrograph 
agreed quite well with the measured one. The three peaks 
were all simulated, although the initial computed peak was 
about 16% lower than measured. The shape of the computed 
hydrograph was very similar to the measured one. The computed 
volume of storm water was about 1% less than measured, represen- 
ting only O.OI" over the entire watershed. The time to peak 
of the simulated hydrograph agreed well with the measured 
hydrograph for the first peak, and was slightly less than 
measured for the second and third peaks. 

Storm of July 3, 1971 

In this simulation (Figure 13), the computed peak flow and 
time to peak were similar to those of the measured hydrograph. 
The simulated peak was approximately \(>"K less than the measured. 
The total volume of runoff was only about SQ% of the measured 
volume, representing 0.07" of rainfall on the watershed. This 
was attributed to error in the rainfall data. For such a 
short rainfall event, differences between the weighted average 
of three gauges and the actual basin rainfall could be 
large enough to account for this discrepancy. The shape of the 
rising and falling limbs of the computed hydrograph were quite 
similar to those of the measured hydrograph. 

Storm of July 15, 1971 

Both peak flow and volume of runoff were well simulated in this 
double event (Figure 14). The rounded, flat shape of the 
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observed hydrograph was closely simulated by the model, even 
during the rain-free period between the major rainfall pulses. 

Storm of July 17, 1971 

The computed peak flow (Figure 15) was within ]$% of the measured 
peak, but the marked double peak in the observed hydrograph 
was not well defined in the simulated hydrograph. This may 
be attributed to the dampening effect of weighted average 
procedure produced on a storm of this magnitude with such sporadic 
bursts of high intensity rainfall. The computed volume was 
35^ less than measured, representing a difference of O.OV 
of rainfall on the watershed. The shape and timing of the 
rising and falling limbs were similar to those of the measured 
hydrograph. 

Storm of July 28, 1971 

in this simulation (Figure 16), the computed peak flow almost 
exactly corresponded to the observed value, and occurred only 10 
minutes later than the observed. The computed volume was 2^% 
less than observed, representing 0.02" of rainfall on the water- 
shed. The sharply peaked computed hydrograph did not closely 
resemble the flatter, broader observed hydrograph. It was felt 
that errors in the rainfall hyetograph were responsible for 
this discrepancy. 

Storm of September 5, 1971 

In this simulation (Figure 17), the peak flow agreed very well 
with the measured value, but was slighly advanced in time. 
Surcharging was computed by the model in some of the upstream 
conduits. The SWMM model simply stores excess pipe flows in 
the upstream manhole of each surcharged conduit, and releases 
them later or when pipe capacity become available. This explains 
the rapid decay of simulated flow after the peak, and the 
levelling off at between 20 and kO cfs as the surcharged flow 
is released to the outlet. This simulation illustrates that 
even with a limited TRANSPORT routine, the SWMM can give a 
reasonably good result under moderate surcharging. 
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h.k.2 West Toronto 

For the detailed simulations, this large 2330 acre CcTtchment 
was subdivided into 45 subcatchment s (Figure 18). Rainfall data were 
available from a rain gauge within the catchment. These were averaged 
with the data from another gauge outside the area, since it was felt 
that a single gauge could not adequately represent so large an area. 
As discussed in Chapter 3, the accuracy of the flow measurements was only 
fair, and it was necessary to adjust the overflow measurements to allow 
for intercepted and diverted flow volumes. Nevertheless, it was felt 
that this estimated flow could be used for comparative purposes. All 
of the comparisons discussed below represent total simulated storm 
flow versus the corrected total storm outflow. Five rainfall events 
were selected for the detailed simulations, representing a range of 
storm intensities for which adequate measurements existed. 

Storm of June 22, 1973 

The simulated flows agreed very well with the observed hydrograph 
(Figure 19). The computed peak was slightly less than observed, 
as was the computed volume of runoff. Default values were 
used, except for the percentage of Impervious area with zero 
detention, which was increased to kS'^- to advance the rising 
limb. The computed volume was about 16'-'' less than the 
observed. 

Storm of September 23, 1973 

The simulated peak flow and volume of runoff agreed reasonably 
well with the observed values, although the computed hydrograph 
lags by about I5"20 minutes (Figure 20). This may be attri- 
buted to rainfall recording inaccuracies. 

Storm of October 2, 1973 

In this simulation, a good agreement between measured and 
observed values was ootained for both peak flow and volume of 
runoff (Figure 21). The simulated values were only 5'^ less 
than those observed. Default values were used in this simulation. 
The shapes of both hydrographs are also quite similar. 
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Storm of August 1 , 1973 

Extensive surcharging in upstream pipes occurred in the 
simulation of this event. This resulted in the truncated shape 
of the computed hydroqraph, as surcharged flows were fir<;t 
stored in manholes and later released (Figure 22). The shape of 
the rising limb before surcharge occurred, however, closely 
resembled that observed. 

Storm of May 10, 1973 

This simulation also produced surcharging resulting in the 

truncated shape of the outflow hydrograph (Figure 23)- The 

rising limb of the computed hydrograph is advanced by about 
10-15 minutes compared to the observed. 

k.k.} Brucewood, North York 

Recent storm flow measurements were available for the '48 acre 
Brucewood area (Figure 2k), and flow simulations were conducted using the 
default values for the various runoff parameters. All the recorded 
data were for relatively small storms. 

Storm of May 1^, IS?'* 1 

Figure 25 shows the results of simulation for this event, 
which produced one of the highest recorded peaks. The 
computed peak flow was some 18' lower than that recorded, 
and a general lag of 5-'0 minutes was observed between the 
measured and computed hydrographs. The shape of the computed 
hydrograph agreed well with the observed, including a slow 
initial buildup and first peak before the major flow event. 
The volume of runoff computed was approximately 16% less 
than that recorded. 

Storm of August 11, 1975 

Figure 26 shows the results of simulation for this small 
rainfall event. The computed and measured peak flows were 
within 81 of one another, as were the runoff volumes. The 
shape of the computed hydrograph was quite similar to that 
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recorded; however, the falling limb Wfis delayed result inq 
in a broader hydroqraph. 

Storm of August 29. 197- 



A very close agreement between measured and computed peak 
flows is evident in Figure 27- fhe initial p^^aks of this 
event were also we i ! represented by the SWMM, although 
a more rapid flovj recession follov/ing peak fiov/s was observed 
in the measurements. The total computed volume v/as 
approximately lO-' greater than that recorded. 

Storm of November 20, 197^ 

In this simulation (Figure 28), the computed peak flow was 
3 cfs or 33'' less than the recorded value. The first peak 
recorded was larger and better defined than the computed peak. 
The computed volume of flow v«s some 25'-^ greater than recorded. 

k . 5 Comparison of Models 

^4.5.1 RUNOFF flow simulat ion 

The two largest storms used in the SWMM flow simulations for 
the Bannatyne district {June 19 and September 5, 1970 were selected 
as the basis for a comparison of the quantity simulation routines of the 
SWMM, WRE and the Dorsch HVM. The RUNOFF models in the WRE and SWMM 
are similar. Consequently, hydrographs from subwatersheds were compared 
only for the SWMM and Dorsch HVM models. Both models were applied using 
the fine discretization shown in Figure II, The peak flows generated 
in each of the ^\ subcatchments by the SWMM and HVM Dorsch models are 
indicated on Figure 29 and Figure 30, for storms of June 19 and September 
5, respectively. For the former, it appears that the peak flows simulated 
by the Dorsch HVM were similar, but slightly higher than their SWMM 
equivalents. However, for the second storm the results are more evenly 
distributed about the ^5 degree line. Flows draining from individual 
subcatchments were not recorded in the Bannatyne sampling program and, 
consequently, an absolute comparison of the models with measurements 
was not possible. A previous comparison of SWMM RUNOFF and Dorsch HVM 
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with measurements carried out for the 10 acre area, Oakdale, Illinois, 
indicated that differences were not significant from a practical 
view-point. The results of the present comparison confirm that, for 
small watersheds, peak flows simulated by the two runoff models are 
very close. Both models use the kinematic wave approximation for overland 
flow, and pipe routing effects for small areas are not considered 
signi f icant . 

li.S.2 TRANSPORT flow simulation 

The following procedure was used in order to obtain an objective 
comparison of the TRANSPORT flow simulation routines of the SWMH, Oorsch 
HVM and URE models: 

(a) The Bannatyne combined sewer district was discretized 
and the watershed and sewer system data were reduced 
according to the Dorsch HVM procedure. 

(b) The inlet hydrographs to the sewer network were computed 
using the Dorsch HVM routine. 

(c) The TRANSPORT routines of the three models were used to 
route the flow through the system, using the HVM 

inlet hydrographs as input. 

(d) The procedure was carried out for two storm events - 
June 19, which was of low intensity and long duration, 
and September 5, which was of medium intensity and short 
durat ion. 

The results of these comparative simulations are presented in 
Figures 31 and 32. For the first storm there was little difference between 
the simulations, which all agreed well with measurements. For the second 
storm, the system was under surcharge. 

The SWMM hydrograph shows the typical truncation effect due to 
the storage and later release of excess runoff. The WRE and Dorsch models 
reproduced the shape of the recorded hydrograph fairly well. 

The models used were uncal ibrated . However, it appears that 
with limited calibration the WRE and HVM models would accurately reproduce 
the measured hydrograph , whereas in the second, surcharged situation the 
SWMM would not. Limitations of the SWMM as compared to WRE and Oorsch HVM 
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became apparent at high flows when backwater effects or surcharge were 

important. The original SWMM was mainly a planning model develed 

for the study of pollution abatement methods. Comparatively rare events, 

such as 7 to 10 year storms, which are important for the design of 

relief sewers are not necessarily the critical events for the design 

of pollution abatement facilities {storage and/or treatment). Therefore, 

the discrepancies discussed above do not necessarily affect the use 

of SWMM as planning tool. 

For analysis of the existing systems for rare storm events 
or for final design applications, consideration of backwater and surcharge 
may result in savings in the capital cost of new works and the SWMM 
TRANSPORT routine should not be used. In this situation, input hydrographs 
generated using the SWMM can be routed by the Dorsch HVM or WRE TRANSPORT 
models, with superior results. The present work indicates that both 
more sophisticated models are operative and useful in practical applications. 
The non-proprietary WRE TRANSPORT model is currently being incorporated 
by the University of Florida as an optional routine in the SWMM. 

i| . 5 . 3 RUNOFF and TRANSPORT flow simulation 

The limitations of the SWMM TRANSPORT routine in simulating 
flow under surcharged conditions have been discussed. In this section 
a global comparison of the Dorsch HVM and the SWMM is presented. The 
storms of May 10 and September 23, 1973 for the West Toronto area were 
used for these simulations. The May 10 storm caused surcharging, whereas 
the storm of September 23rd did not. The results of the simulations 
of each of the models are shown in Figures 33 and 3'^. 

The Dorsch HVM simulations were based upon a very detailed 
discretization of the 2330 acre catchment that used approximately 2000 
conduits. The SWMM simulation, on the other hand, was based upon 
a discretization using ^45 subcatchments and 71 conduit sections. These 
simulations were conducted with uncalibrated models . The default values 
in the SWMM were employed, and the Dorsch HVM model was run with the standard 
values used in their overall sewer system study for the City of Toronto. 
The measured flows, as discussed in Chapter 3, included approximations 
for diverted flows, and should not be used for absolute model comparisons. 
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For the storm of September 23 {Figure 33), the SWMM model produced 
a sfightly higher peak flow and larger volume than did the Dorsch HVM, 
although the shape of the two computed hydrograph were similar and 
resembled that of the measured on©. No s.ignificant surcharging occurred 
in the SWMM computations. 

For the storm of May 10 (Figure 'ik) , however, the SWMM computed 
hydrographs show the truncation effect due to stored surcharged water being 
released later on, while the Dorscii HVM hydrograph shows a much closer 
resemblance to the shape and magnitude of the measured hydrograph. The 
Dorsch simulation was uncal ibrated , and could be made to represent the 
measured hydrograph very closely with minor calibration. 

These results further indicate that the SWMM can produce a 
response under unsurcharged conditions which is not significantly different 
from that of a more sophisticated model using a far greater amount of 
detail in the computations. Under surcharged conditions, the Dorsch model 
gives a superior simulation to that of the SWMM. 

4,$ Simulation Costs 

Simulation costs for various models have been discussed in 
recent studies [7,5^]- Meaningful comparisons between models are clouded 
by differences in computers used and the rate schedules for on-line and 
off-line equipment. Other factors which affect the cost of a simulation 
are the amount of printout, the extent of surcharged flow simulations, 
and the number of time steps employed. Table 5 summarizes typical run 
costs for different jobs, on the basis of the numerous runs done under 
the terms of this study. The SWMM simulations were performed on an 
IBM 370 model 168 computer at a basic rate of $&M0 per hour. The WRE 
model runs were on a CDC 6600 compjter at a rate of $1 ,080 per hour. 
The Dorsch HVM simulations were conducted by Dorsch, and the computer 
costs reported represent amounts for computer time only* 

Most SWMM simulations involve the RUNOFF and TRANSPORT blocks 
only, although in the future the TREATMENT and RECEIV blocks will be 
increasingly applied. On the basis of the simulations done in this 
study, the following has been prepared as an approximate guide to the 
Cuser in the estimation of the costs of running the SWMM: 
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For runs Involving RUNOFF and TRANSPORT only: 

Cost in Dollars = total no. of elements x 

no. of time steps x $0.0025 

For runs involving RUNOFF, TRANSPORT, STORAGE/TREATMENT 
and RECEIV: 

! 

Cost in Dollars = total no. of elements x no. of 

time steps x $0.0050 



These approximations are based upon the computer system used in this 
study, as discussed above. Several runs will give any user a feel for 
costs on his own computer. I 

TABLE 5. TYPICAL SIMULATION COSTS 
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k.7 Conclusions 



(a) A sensitivity analysis has indicated the relative influence 
of the major parameters in both the RUNOFF and TRANSPORT 
blocks. This analysis has further indicated that the 
existing SWMM default values should be retained^ 

(b) The Dorsch HVM and SWMM generate similar subcatchment 
overland flow hydrographs for a given rainfall pattern. 
This result was expected since both models employ the 
Kinematic Wave approximation for overland flow routing. 

(c) Under normal flow conditions, with no surcharging of the 
sewer system, the Dorsch HVM, WRE and SWMM TRANSPORT 
routines produce similar hydrographs for a given set of 
inlet hydrographs. The quasi-steady dynamic wave form- 
ulation of the SWMM is fairly simple. However, good 
results are produced for conduits in which there are no 
excessive backwater effects. 

(d) In the case of drainage systems under surcharge aiid with 
backwater conditions, the SWMM does not correctly reproduce 
the shape of the recorded hydrograph. The Dorsch HVM and 
WRE models use a pressurized flow routine for surcharged 
conditions, thereby reflecting the increased capacity of 
pipes under surcharge. 

In addition, the solution of the dynamic wave equation 
over the entire sewer system at each time step allows 
for a proper simulation of the dynamic effect of the 
backwater condition. 
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5. SWMM dUALITY SIMULATIONS FOR SELECTED CANADIAN WATERSHEDS 
5. I Genera] 

Studies described in the literature concerning urban storm water 
quality have identified various factors contributing to storm water 
pollution [25, ^2, hi, SS-B^t, 67]- GeneraHy, pollutants can enter 
urban runoff from three major sources - the land surface (both pervious 
and impervious), catch basins, and from sewers in a combined drainage 
system. 

Water quality models for urban runoff can be classified into 
three main categories: empirical, statistical and analytical. A 
discussion and comparison of several existing urban runoff quality models 
is given in a literature survey, presented in Volume II. A descriptive 
summary comparison of various models is presented in Table 6. In general, 
statistical and empirical models cannot be applied for watersheds, other 
than those for which they were developed, with a great amount of 
confidence in the results. Analytically based models may be applied in 
a wider variety of situations due to their more comprehensive nature. 
The SWMM quality routines are based principally on analytical formulations, 
but where the state-of-the-art does not justify these, some empirical 
relationships have been built into the model. To date, testing of the 
SWMM quality models has been restricted to a few specific watersheds 
in the U.S. (for examples see ['t?] and Volume II). In most cases, their 
applications were limited to only few events and were directed towards 
preliminary calibration of model parameters and definition of default 
va 1 ues . 

This section of the study examines the applicability of the 
SWMM quality routines using available quality measurements from selected 
Canadian watersheds. Emphasis has been placed on an assessment of model 
sensitivity to variations in quality parameters and default values, and 
on developing procedures for calibration. Some of the recommendations 
concerning modifications to the model subroutines and input data format 
have been discussed with representatives of the University of Florida 
and have since been implemented, constituting options in the latest 
version of SWMM, The requirement for a more flexible approach to calibra- 
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TABLE 6. COMPARISON OF WATER QUALITY FEATURES OF VARIOUS MODELS 



Model 


Type of Model 
Empirical 


Quality 
Constituents 


Dry weather 
Qual ity 

No 


Storm runoff 
Quality 

Empirical 
relationship 


Qual ity 
Routing In 
Channels 


Sedimentation 

and Scour 

in Channels 


Quality 

Reactions 

In Channels 


Stanley 
[63] 


BOD 


No 


No 


No 


Arnett et al 
[70] 


Semi -empirical 


BOD, suspen- 
ded solids, 
phosphates, 
nitrogen 


No 


Linear 

regression 

equations 

1 


No 


No 


No 


AVCO 

Corporation 
[62] 


Statistical 


BOD, COD 
total and 
suspended 
solids, 
nitrogen, 
phosphates 
total and 
fecal coli- 
forms 


No 


Regression 
techniques, 
1 inear and 
non-linear 
relation- 
ships 


No 


No 


No 


Sartor & 
Boyd 
[56] 


Statistical 


Total 
solids 


No 


Regression 
techniques 


No 


No 


No 


Brands tetter 
[74] 


Statistical 


7 arbitrary 

conservative 

constituents 


Hourly, 
dai ly and 
seasonal 
patterns 


Linear 

regression 

equations 


Pure ad- 

vection, 

mixing 

between 

successive 

time steps 


No 


No 
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Quality 


Sedimentation 


Quality 






Quality 


Dry weather 


Storm runoff 


Routing In 


and Scour 


Reactions 


Model 


Type of Model 
Analytical 


Constituents 


Quality 
No 


Quality 
Non- linear 


Channels 


in Channels 


in Channels 

No 


Corps of 


BOD, suspen- 


No 


No, but land 


Engineers 




ded and 




function of 




surface ero- 




(STORM) 




settleable 




catchment 




sion by univer- 




[66] 




solids, 
nitrogen and 
phosphates 




characteris- 
tics, pollu- 
tant accumu- 
lation and 
runoff 




sal soil loss 
equation 




Hydro comp 


Analytical 


17 consti- 


Inflow con- 


Non-1 inear 


Pure advec- 


No 


Various reac- 


[65] 




tuents 


centration 


function of 


tion, 




tions and 






including 


provided as 


catchment 


weighted 




interactions 




- 


water temp- 
erature 


input data 
or non- 


characteris- 
tics, pollu- 


mixing be- 
tween suc- 


















linear 


tion accumu- 


cessive 












function of 


lation and 


time steps 












inflow 


runoff 














hydro graph 










EPA 
















1-Original 


Analytical 


BOD, sus- 


Hourly, 


Non-1 inear 


Pure advec- 


Suspended 


First-order 


Metcalf & 




pended 


daily and 


function of 


tion, mix- 


solids, con- 


decay without 


Eddy 




solids, coli- 


seasonal 


catchment 


ing between 


sidering 


interactions 


Version 




forms. 


patterns 


characteris- 


successive 


particle 


for BOD 


[47] 




nitrogen and 


provided as 


tics, pollu- 


time steps 


size distri- 








phosphates 


input data 
or computed 
from land 
use 


tion accumu- 
lation and 
runoff. 
Different 




bution 




: 






■ 


equation used 
for pollutant 
removal 
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Quality 


Sedimentation 


Quality 






Quality 


Dry weather 


Storm runoff 


Routing In 


and Scour 


Reactions 


Model 
2-Water 


Type of Model 
Analytical 


Constituents 
COD, Settlea- 


Quality 
Same as EPA 


Quality 
Equation 


Channels 
Same as EPA 


in Channels 


In Channels 
Same as EPA 


Same as EPA 


Resources 




ble 




solved in a 








Engineers 




Solids and 




different 








Version 




oil and 




way; no 








[75] 




grease in 
addition to 
above con- 
stituents 




availability 
factors for 
suspended and 
settleable 
solids 








3-WRE 


Analytical 


Same as EPA 


Same as EPA 


Two methods 


Same as EPA 


Same as EPA 


Same as EPA 


Version 








for sus- 




also inclu- 




modified 








pended solids 




des land 




by Uni- 








computation 




surface ero- 




versity of 












sion by the 




Florida 












universal 




[71] 












soil loss 
equation 




Dorsch 


Analytical 


BOD, Settlea- 


Included but 


Empirical 


Continuity 


Not presently 


No pollutant 


Model 




ble sol ids 


not documen- 


washoff 




included 


decay 


[6 ] 




other pollu- 
tants cur- 
rently being 
added 


ted 


equation 








UCURM 


Analytical 


BOD, Sus- 


No 


Non-linear 


Pure ad- 


Suspended 


No pollutant 


[6^] 




pended and 




function of 


vection 


solids. 


decay 






settleable 




catchment 


mixing 


considering 








solids 




characteri- 
stics, pollu- 
tion accumu- 
lation and 
runoff 


between 
successive 
time steps 


partial 
size distri- 
bution 





tion resulted in the formulation of a generalized surface quality model, 
which is described in Chapter 10. 



5-2 



Review of Some Basic Assumptions 



The state-of-the-art is such that there is not sufficient 
information available to model all the complex relationships between 
urban storm water contaminants, their pollutional characteristics and 
the manner in which they are transported during storms [25, 55, 563- 
Therefore, while the mathematical relationships used in the SWMM are 
the most comprehensive currently available, they are, nevertheless, based 
on simplifying analytical assumptions with some parameters determined 
by statistical means. The relationships presently used in the modelling 
of surface runoff in the SWMM are summarized in Figure 35- These basic 
equations, which are also used in other models, are discussed in more 
detail in the SWMM User's Manual [71]- However, a brief review is 
pertinent to the discussion contained in the following sections. 

The initial mass of dust and dirt accumulated on the surface 
is estimated as a function of the number of dry days, and street cleaning 
frequency and efficiency prior to the storm event. The initial amounts 
of various pollutants, Po , in the dust and dirt are assumed to be in 
the proportions summarized in Table 7 for various land use classifications, 
The most important parameter used for surface runoff quality computations 
is the initial estimate of the total amount of each pollutant, Po , 
on the surface. High initial values of Po will give high initial 
washoff rates. Reference to the equations in Figure 35 indicates that 
variations in the number of dry days and street cleaning practices 
directly control the initial values of Po. 

An exponential decay equation for surface washoff of pollutants 
assumes that the amount washed off in each time step is proportional 
to the total amount remaining on the surface. The rate of washoff is 
controlled by the runoff rate for each time step and an exponent, b. 
An average b value of ^4.6 is assumed in the model (this is discussed 
in more detail in Section 5-3). 

The most recent version of the SWMM, published by the University 
of Florida, contains an optional method for computing suspended solids 
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ALGORITHM OF THE SWMM 
SURBkCE RUNOFF QUALITY MODEL 

FIG. 35 
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TABLE 7. DEFAULT VALUES FOR DUST AND DIRT ACCUMULATION USED IN SWMM^ 



o 



Land Use 



1 . Single Family 



2. Multiple 
Family 

3. Commercial 



4. Industrial 



5. Undeveloped ° 
or Park 



Dust and Dirt 
Accumulation 
lb/day/100 ft Curb 



n.7 



2.3 



3.3 



4.6 



1.5 



(2. or 



(2.0)' 



(2.9)^ 



(6.0)* 



(K0)< 



Coliforms^ 
1.3 X 106 

2.7 X 10^ 



1.7 X 10^ 



1.0 X 10^ 



0.0 



Susp. 

Solids 

1 000 . 
(lll.O)d 

1000.0 
(80.0) 

1000.0 
(170.0) 

1000.0 
(67.0) 

1000.0 
(111.0) 



mg Pollutant Per gm of 
Dust and Dirt 



Sett. 
Solids 

100.0 



100.0 



100.0 



100.0 



100.0 



BOD 



5.0 



3.6 



7.7 



3.0 



5.0 



COD 
40.0 

40.0 

39.0 

40.0 



20.0 



_H__ 
0.48 

0.61 

0.41 

0.43 



0.05 



PO4 
0.05 

0.5 

0.7 
.03 
.01 



Grease 



1.00 



1.00 



1.00 



1.00 



1.00 



a - Most values are based on 1969 APWA Report [61], and included in [7!] 

b - Values for undeveloped and park lands are assumed. 

c - Units for Coliforms are MPN/gram. 

d - Values used in the WRE version of SWMM [75]. 

e - Data from [76]. 



doncent rat ions [71]- As indicated in Figure 35, the original method (ISS=0) 

is consistent with computations for other pollutants since it uses the 

exponential decay equation. However, an availability factor, computed 

as function of runoff intensity, is also incorporated to account for 

the fact that at lower runoff intensities, insufficient energy is available 

to allow the solids remaining on the surface to be treated as completely 

available for washoff. For values of runoff intensity greater than 

0.7 inches per hour, the availability factor is assumed to 1.0 (i.e. 

all solids are available for washoff). The availability factor is used 

only for suspended solids computations. 

The second optional method (lSS=l) for computing suspended 
solids is an empirical equation developed from San Francisco data during 
the initial SWMM development [kj] . The physical significance of many 
parameters is difficult to assess. However, three major underlying 
assumptions are: 

(a) The amount of SS removed during each time step is 
essentially a function of the ratio of the amount 
remaining on the ground at the beginning of the time 
step (Po) divided by the initial amount on the ground 
at the beginning of the simulation period, (Po;). 

(h) The amount removed is a function of a removing coefficient, 
CC, which varies from 0.9 to 0.25, decreasing over a period 
of about five hours from the start of the quality 
s imu lat ion. 

(c) As programmed, the equation could be modified to operate 
over three ranges of runoff intensity. (The reader 
is referred to Figure 35 for clarification.) This 
assumption can sometimes result in some discontinuity 
in the computation of the SS pol 1 utograph. 

For DWF computations, either average values of SS, BOD and 
col i form concentrations can be supplied by the user or these can be 
estimated in the model, based on input data describing population and 
land use characteristics ['»7] • The program also allows for the hourly 
variation of these parameters to be modelled, either by using program 



103 



default values, or by supplying hourly rates. Deposition, scour and 
pollutant decay can also be modelled in the TRANSPORT block. 



5.3 



Sensitivity Analysis 



The SWMM storm water quality simulations are dependent upon 
several parameters. Some of these must be supplied by the user while 
others may be supplied as an option and are also built into the model. 
The internal default values provided in the model are generally based 
on average measurements. 

A sensitivity analysis involves the systematic independent 
variation of variables used in the computations in order to ascertain 
the effect of such variations on the results and to decide whether or 
not the model responds in a logical manner to these changes. The 
sensitivity analysis indicates how changes in specific variables will 
affect the results, and so a knowledge of model sensitivity and a feeling 
for the allowable range over which the more influential parameters may 
be varied can lead to a logical adjustment of the model and to better 
simulations. 

The relative sensitivity of the model to changes in the following 
parameters was assessed in this study in order to assess the potential 
for cal ibrat ion: , 

(a) the exponent, b, in the exponential washoff equation 
(see Figure 35) ; | 

(b) the use of the two available options for SS computations - 
empirical (lSS=l) and exponential {ISS=0); 

(c) street cleaning frequency and number of dry days, 
catch basin BOD; i 

(d) pipe slopes in the TRANSPORT quality computations; and, 

(e) the specific gravity of solid particles in TRANSPORT 
quality computations. , 

A hypothetical test area of 15 acres was used in the sensitivity analysis 
for the surface quality routines (Figure 37)- 
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5.3.1 Sensitivity to washoff exponent "b" 

The developers of the model defined the exponent, b, used in 
the washoff equation {see Figure 35) in the following manner: 

"it was assumed that a uniform runoff of 1/2 inch per hour would wash 
away 90 percent of the pollutant in one hour". (This leads to a value of 
b of ^.6.) "Other assumptions for determining b could be made but this 
original assumption has proven satisfactory in all test applications to 
date for urban areas" {[^7], p. 178). 

This assumption was examined by plotting b versus runoff intensity for 
various removal rates as shown on Figure 36- Reference to this figure 
indicates that a value of b = ^+.6 appears to be a "representative" value. 
For runoff intensities in excess of 0.7 inches/hour, b = ^.6 gives 
removal rates of 100^ in one hour. For smaller runoff intensities, 
and short duration storms, it is possible that a significant amount of 
pollutant may remain on the ground at the end of the storm if a value of 
b of h.6 (or less) is used. For a given rainfall intensity, higher 
values of b will result in a greater amount of total pollutant removal 
in the first few time steps and would hence give higher initial pollutant 
concentrations. 

5.3-2 Sensitivity to type of washoff equation used 

The model includes two equations for suspended solids calculations. 
These are referred to as (1SS=0) and (lS5=l) and may be written as: 

POFF = AV X Po (l-e"^"''l ' ^^) (ISS=0) 

and POFF = ^x cc^ x {AxE + Bxe'') (lSS=l) 

i 

where: POFF = pollutant load washed off in time step, 
At 

AV = an availability factor (see Figure 35) 

Po = surface pollutant load at start of time step 

rj = runoff rate 

b = constant (=4.6 in SVWH) 

Poj = initial pollutant load 

cc = a removing coefficient (maximum = 0.9) 

A,E,B,D = empirical coefficients 
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An inspection of the equations suggests that low or moderate values of 
Po (tSS=l) will tend to give a greater value for POFF in the first time 
step, unless the availability factor, Av, equals 1.0, since the initial 
ratio Po:Po. will always be unity. (The value of AV used in SWMM is 
calculated from the runoff rate and is usually much smaller than 1.0.) 
Also, as Po increases, due to higher pollutant loading rates or longer 
periods of pollutant accumulation, the differences between the two options 
might be expected to decrease. For large values of Po the exponential 
equation might be expected to give higher initial washoff rates, since 
a ratio controls the washoff in the empirical equation. 

The simulation results for the hypothetical area presented 
in Figure 38 indicate that, for SS computations when the initial surface 
loading is relatively small, the empirical equation predicts concentrations 
which are significantly higher than those given by the exponential equation. 
However, for a larger initial surface loading, the exponential decay 
equation gives higher concentrations. The use of the exponential equation 
of SS computations is consistent with computations for other pollutants. 
However, as discussed in Section 5-^, for low initial pollutant loadings 
(as a rule of thumb, say lower than 10-15 dry days prior to the storm), 
calibration to measurements may be possible only by using the empirical 
equation. 



5.3.3 



Sensitivity to initial pollutant accumulation 



The influence of the initial pollutant accumulation on the 
simulation of pollutants washed off is also indicated in Figure 38 
for the hypothetical test catchment of Figure 37. In this example, the 
default pollutant loading rates and compositions were assumed. The 
number of dry days preceding the storm was varied in order to simulate 
different initial surface pollutant loads. A cleaning frequency greater 
than the number of dry days was assumed to avoid modelling of street 
cleaning. The results indicate that the total amount of pollutants 
washed off during this idealized storm is approximately a linear function 
of the preceding number of dry days. ' 

Obviously, street cleaning reduces the amount of surface pollu- 
tant available for washoff and will, consequently, reduce concentrations 
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in the runoff. The simulation results shown on Figure 39 indicate the 
effect of variations in the street cleaning frequency. 

S.'i.h Sensi ti vi ty t o catch basin BOD concentration 

The sensitivity of the BOD calculation to the BOD concentration 
in the catch basins was examined by varying the catch basin BOD concen- 
tration for the 15 acre hypothetical area. A typical catch basin density 
of one per acre and a volume of 30 ft3 were assumed. The BOD concentration 
was varied from 25 to 150 mg/1 for an accumulation period of three days. 
As indicated below in Table 8, the pealc BOD concentrations and total 
loads vary by a relatively small amount with the catch basin concentration, 
(see also ['l/l , Table 5-18). However, for relatively low intensity 
storms where the initial surface loading is small, some calibration may 
be possible by varying the catch basin BOD concentration. 

TABLE 8. EFFECT OF VARIATION OF CATCH BASIN BOD 



Catch basin 


Peak BOD 


Total Mass 


Concentration 


Washed off 


of BOD 


BOD (mg/1) 


(mg/l) 


(lb) 


is 


1A.104 


9.06 


100 


17.227 


9.58 


150 


20.879 


10.60 


Sensitivity to TRANSPORT block p 


arameters 





The sensitivity analyses performed for the TRANSPORT block 
quality routines used data from the Bannatyne combined sewer district 
in Winnipeg (5^2 acres). 

The accumulation of a bed of solids in the sewer during periods 
of dry weather flow was found to vary considerably, depending on the pipe 
slope and number of dry days. For example, pipe slopes in the order 
of 0.0001 ft/ft resulted in simulated accumulations up to 100 times 
greater than those using pipes with slopes of the order of 0.001 ft/ft, 
depending on pipe flow, diameter and length. 
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According to the procedure used in the model, the accumulation of 
solids in the sewers is directly related to the number of days during 
which no storm water flow occurred. Thus, it is evident that the number 
of dry days is also an important parameter for combined systems, as 
illustrated in Table 3- Variations in the number of dry days modelled 
for combined systems can help to show the effects of sewer flushing 
as a pollution control method. 

Routing and decay of nonconservat ive pollutants is modelled 
in the TRANSPORT conduits. However, "some doubt exists as to the 
significance of quality change, other than mixing, during rainfall- 
runoff periods. It has been suggested that the time of transit through 
the sewer system is too short to consider bacterial decay or reduction 
of BOD with time", ([4?], Figure 13-2). | 

The sensitivity to the BOD decay rate was investigated by 
reducing the decay coefficient, D, from 0.20 to 0.0. The Bannatyne June 
19 storm was run for these two values (simulation period 'tOO minutes), 
and the total amounts of BOD flushed out were 5,064 and 5,050 lb, respectively, 
a difference of less than 0.5^. The difference between measured and 
computed BOD values is much greater than this and, hence, D has essen- 
tially no effect . I 

The specific gravity of the solid particulate material in the 
dry weather flow in the TRANSPORT quality routine is fixed at 2.7- The 
effect of reducing this to 2.5 was investigated using the storm of June 19, 
1971 for the Bannatyne district. A reduction in the total accumulated 
suspended solids load in the pipes from 6030 to 't538 lb was observed. 
However, specific gravity is not considered to be useful for calibration 
purposes, unless the total solids discharge is consistently overestimated, 
since a value of 2.7 appears to be the maximum realistic value that 
could be used. 



5.3.6 Conclusions of sensitivity analysis 

In general, the SWMM quality model appears to respond in a 
logical manner to variations in input parameters. However, neither option 
for the computation of suspended solids from surface runoff seems generally 
applicable and the model is somewhat difficult to calibrate. This is due 
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TABLE 9. SIMULAHON RESULTS FOR BANNATYNE WINNIPEG 



STORM 




Recorded 
Flow 
cu ft 


Computed 
Flow 
cu ft 


Ratio* 


BOD 

Recorded 

lb 


BOD 

Computed 

lb 


Ratio* 


SS 

Recorded 

lb 


SS 

Computed 

lb 


Ratio* 


June 19/71 
ISS=1, 8 days 


1) 


261,300 


241,716 


0.93 


3051 


2190 


0.72 


16,585 


11,378 


0.69 


IS5=0, 8 " 


2) 










2078 


0.68 




9,148 


0.53 


ISS=0.11 " 


3) 










2138 


0.70 




11,734 


0.70 


July 15/71 






















ISS=1, 5 days 


1) 


111,600 


101,131 


0.91 


1450 


1124 


0.78 


5.382 


4,613 


0.86 


ISS=0, 5 " 


2} 










1114 


0.77 




4,318 


0.80 


ISS=0.12 " 


3) 










1161 


0.80 




9,147 


1.69 


July 17/71 






















ISS=1, 2 days 


1) 


171,960 


143,184 


0.83 


2285 


1085 


0.48 


9,423 


3,823 


0.41 


ISS=0, 2 " 


2) 










1012 


0.45 




2,453 


0.26 


ISS=0,10 " 


3) 










1120 


0.49 




9,196 


0.98 


July 28/71 






















ISS=1, 8 days 


1) 


227,700 


159,136 


0.7 


2224 


1485 


0.67 


24,988 


9,306 


0.37 


ISS=0, 8 " 


2) 










1434 


0.65 




7,886 


0.32 


ISS=0,11 " 


3) 










1479 


0.67 




10,484 


0.42 


Sept. 5/71 






















ISS=1,16 days 


1) 


447,240 


373,879 


0.84 


27,141 


1616 


0.05 


84,857 


26,205 


0.31 


ISS=0,16 " 


2) 










1620 


0.06 




27.121 


0.32 


ISS=0,39 " 


3) 










3116 


0.11 




64,798 


0.77 



* Ratio of computed to measured values 



both to its complexity and the fact that some values that could logically 

be varied are fixed in the proqram. 

I 

S«ij SWMM Surface Runoff Quality Simulations 

During the study period, there were very few Canadian data 
available for the study of surface runoff quality. However, for testing 
purposes, the preliminary data from the Brucewood catchment, located In 
North York, Toronto, were used. The catchment drains an area of ^8 
acres with an average imperviousness of about h7%. The predominant land 
use in the area is the modern single family unit. Figure ^0 illustrates 
the detailed schemat izat ion and pipe network used to model the Brucewood 
catchment. The Brucewood quantity and quality monitoring program and 
associated measurement problems are fully documented elsewhere [72]. 
Measured quality parameters include SS, BOD and coliforms. Various problems 
with the original automatic water quality sampler were encountered 
early in the measurement program. Therefore, quality samples were sometimes 
talcen manually prior to the summer of 197^. Consequently, for several 
events, quality samples are not available at the beginning of some of 
the storms. All available collform measurements were collected manually. 

For the purpose of this study, six storms were selected for 
comparison of measured data with quality simulations using SUMM . The 
number of quality measurements for other storms were generally insufficient 
for any worthwhile comparisons with simulated values. Of the storms 
selected, three had maximum rainfall intensities less than 0.60 inches 
per hour and only the storm of September 11, 1975 had a peal< rainfall 
intensity greater than 1.0 Inch per hour. As indicated in Table 10, 
measured peak flows ranged from 4 to 25 cfs. Also summarized in Table 10 
are the total flow and pollutant amounts recorded during the quality 
sampl ing per iod . 

Initially, simulations were carried out for the six storms 
using the SWMM default values for the dust and dirt loading rates and 
pollutant composition. The procedure outlined in the User's Manual 
was followed to determine the antecedent number of dry days [71]; 
that is, prior days were counted as "dry" if the total volume of the 
antecedent rainfall was less than one inch. 
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BRUCEWOOO WATERSHED 
SUBCATCHMENT 
DISCRETIZATION 
FIG. 40 




TABLE 10. COMPARISON OF POLLUTANT LOADS FOR BRUCEWOOD 
SEPARATE STORM DRAINAGE SYSTEM 



Storm 


Flow 

Peak 
Recor'd 


(cfs) 

ConiD ' d 


Samp] ing 
Period 

(min.) 


Flow (cu ft) 


BOD (lb) 


SS (lb) 


Total Coll form 
(MPN) 










Recor'iJ 


Comp ' d 


Recor'd 


Comp'd 


Recor'd 


Comp ' d 


Recor'd 


Comp'd 


May 14/74 


11.4 


9.2 


30 


8619 


5068 


2.63 


7.89 


49.27 


117.09** 




- 
















(2.45) 




(9.55)*** 

1 




May 16/74 


4.5 


2.7 


70 


7401 


3853 


4.06 


9.72 
(6.83) 


77.66 


77.47 
(3.93) 


0.72x10^° 


0.37x10^° 


Nov 20/74 


8.7 


11.9 


35 


10362 


13402 


12.65 


33.74 
(10.85) 


137.18 


587.38 
(97.47) 


2.4x10^° 


4.03x10^° 


Aug 29/75 


5.3 


5.1 


88 


5653 


6801 


1.31 


9.13 
(5.4) 


21.7 


78.19 

(4.09) 


- 


- 


Sep 11/75 


25.0 


22.0 


80 


49422 


53649 


10.04 


61.06 

(45.86) 


345.2 


875.14 
(605.8) 


- 





All total washoff amounts have been computed for the sampling period only. 

Flow and pollutant amounts have not been included for the storm of Aug. 11/75, since 

the recorded runoff rate is extremely low during the sampling period. 

ISS = 1 



*** ISS = 



N.B. It is not possible to present a similar comparison for the event of 
Aug. 11/75 due to limited data. 



The general trends observed in these initial runs are summarized 
below and form the basis for subsequent model adjustments* 

(a) For events with a low number of antecedent dry days, 
suspended solids concentrations were severely under- 
estimated when using the exponential equation (tSS=0). 

On the other hand, a more reasonable estimate was obtained 
for the storm of September 11, 1975, with 12 dry days. 

(b) BOD concentrations were generally slightly underestimated 
when the exponential equation was used for suspended 
solids computations. (The model considers S% of the 
computed SS to contribute to BOD.) However, simulated 
BOD concentrations were much closer to measured BOD 
values than were the simulated SS to measured values, 

(c) Simulated coliform concentrations were much higher than 
the measured values. (A programming error was detected 
in the coliform equations in the model. This has been 
rectified in cooperation with the University of Florida, 
May 1975.) 

These observations, in conjunction with various calibration 
efforts and the results of the previous sensitivity analysis, indicated 
that improved simulations could be achieved by making the following 
adjustments: 

(a) Using the empirical (!S5=1) option for suspended solids 
calculations when initial surface pollutant loads are low. 

(b) Supplying measured composition of pollutants wherever 
available. This included using the average measured 
value for catch basin BOD concentration of 60 mg/1 
(compared to the default value of 100 mg/l) for all the 
events and the dust and dirt accumulation rates measured 
prior to the events of August 11 and August 29, 1975- 
For these two events, the measured suspended solids 
fraction in the dust and dirt was taken as the percentage 
finer than the 30 mesh from the sieve analysis (['t7J, 
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p. 185). The measured values for these storms are 
indicated on Figures kk and 'tS- (The BOD composition 
measured in the dust and dirt was essentially the same 
as the default value.) 
(c) The MPN of coliforms in the dust and dirt was changed 
to 0.65 MPN per gram. This values was determined by 
trial and error simulations and used for the three 
events where coliform measurements are available for 
comparison. 

The six storms were simulated with the input to the model 
adjusted as outlined. The final results only are presented on Figures 41 
and kS. (Default values are used unless otherwise noted.) For events with 
the antecedent number of dry days fewer than 10, the use of the empirical 
equation (lSS=l) resulted in much better simulation of suspended solids 
concentrations than the exponential equation, although estimated BOD 
concentrations were only marginally improved, if at all. For the storms 
of May 16, 197'*, November 20, 197'* and September II, 1975 (Figures AZ , 
^3 and 46), simulated and measured coliform MPN's were the same order of 
magnitude. These results indicate that coliform estimates for separate 
sewer systems, using the corrected exponential decay equation, can give 
reasonable results compared to measured values, provided some calibration 
is undertaken. 

In general, the results of these adjusted simulations indicate 
that pol lutographs computed using the SWMM are of the same order of 
magnitude as the measured values. 

5. 5 SWMM Combined Sewer Quality Simulations 

Comparisons of measured and simulated pol lutographs were 
carried out using data from the Bannatyne sewer district of Winnipeg. 
This district is served by combined sewers draining an area of about 5^2 
acres. The predominant land uses are commercial and residential in 
about equal proportions. Only five storm events were considered suitable 
for simulation. 

BOD and SS were the only parameters recorded in the Bannatyne 
area. Peak suspended solids and BOD concentrations were 154 to 9,294 
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FLOW AND POLLUTOGRAPHS 
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FLOW AND POLLUTOGRAPHS 
FOR BRUCEWOOD 
STORM OF SEPT. II, 1975 
FIG. 46 
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mg/1 and 25 to 1,720 mg/1 , respectively. These concentrations, which 
are relatively high, are attributed to scour of accumulated solids from 
the fiat sloped pipe network. 

Simulations were conducted using '+1 subcatchments ranging in 
size from 5 to 23 acres. For all simulations, the catch basin 
contribution of BOD was assumed to be 100 mg/1 and the catch basin 
volume was assumed to be 20 cubic feet [71]. As a result of discussions 
with the City of Winnipeg Works personnel, no street cleaning prior 
to each of the storms was assumed. Preliminary simulations indicated 
that the DWF generated by the model as a function of population and land 
use statistics significantly underestimated the measured DWF. Therefore, 
in order to reflect the measured data, an average DWF of 5-6 cfs with 
SS and BOD concentrations of 558 mg/1 and 391 mg/1, respectively, 
was used (as determined from available measurements). The importance 
of accounting for the DWF rn the simulation is indicated in Figures ^7 
to 51 which compare simulated pollutant concentrations from surface runoff 
alone, and simulated concentrations accounting for both surface runoff 
and DWF, with measured pol lutographs . The simulations for surface runoff 
alone shown on the figures correspond to the option {lSS=l) for the 
SS computa t i ons . 

In this particular area, the contribution of pollutants resulting 
from the flushing of sanitary deposits in the combined sewers outweighs 
the surface contribution. Consequently, the choice of the SS surface 
washoff equation has little effect on the outflow pol I utographs . In 
the curves shown in Figures h7 to 51, the empirical option was used, as 
most of the antecedent dry periods were less than 10-15 days. Two 
different methods for estimating the number of antecedent dry days were 
used. The first method summed the previous days on which zero rainfall 
had occurred. The second method determined the antecedent condition 
by assuming the "dry" period extended back to the date where the 
cumulative total antecedent rainfaM exceeded one inch. The second 
method always gives the number of dry days greater than or equal to 
the first method. Only the pol lutographs based on the second method 
are given on Figures U7 to 51 • A comparison of the total mass of pollu- 
tants washed off from the start of the storm to the end of the sampling 
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period, with the corresponding simulated loads, is presented in Table 10. 

In general, the results of the simulation of pol lutographs based on 

the second method for determining the dry period were quite good. However, 

simulated BOD concentrations were consistently lower than the measured 

values. This might be due to the fact that there is likely to be an 

additional contribution of BOD resulting from the solids settled 

in the sewers. Inclusion in the model of a factor representing this 

additional BOD would require some analysis of deposits from a variety 

of areas. 

It would obviously not be possible to obtain reasonable 
simulations without allowing, in the modelling process, for substantial 
deposition of solids in the combined sewers. Field observations noted 
by engineering personnel in the City of Winnipeg have confirmed the 
occurrence of significant accumulation of solids deposits during dry 
periods . 

§4$ ConcI usions 

{a) The state-of-the-art of storm water quality model lirig 
has not advanced as far as that of quantity modelling. 
However, it is possible to obtain simulated pol lutographs 
similar to measured pol lutographs , for both separate storm 
water and combined discharges. 

(b) The sensitivity analysis indicates that logical changes 
in input parameters are logically reflected in the 
simulated results. Consequently, the model can be regarded 
as being a useful tool for comparison of storm water 
management alternatives. However, caution should be 
exercised in using the model for estimates of the absolute 
value of storm water pollutant loads. 

(c) Neither of the options provided in the SWMM for suspended 
solids computations are adequate over a wide range of 
antecedent conditions. When the initial pollutant loading 
is very high, the exponential equation will simulate higher 
SS concentrations than the empirical equation. Conversely, 
when the initial pollutant loading is fairly low, the 
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empirical equation will simulate higher SS concentrations. 
In general, the empirical equation (lSS=l) may facilitate 
calibration when the dry period preceding the event is 
less than 10-15 days. ' 

(d) The quality model in the SWMM is rather difficult to 
calibrate, partly because of its complexity and partly 
because the values of certain "variables" are fixed in 
the program. When the primary concern is surface runoff 
a simplified approach may be justified. The application 
of a generalized SWMM quality model is discussed in 
Chapter 10. I 

(e) Initial conditions have been shown to control the 
quality computations. Consequently, one event quality 
models should be used in conjunction with continuous long 
term simulation models. These would be used to establish 
initial conditions for the more detailed model. (This 

is discussed more fully in Chapter 13.) 
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6. ASSESSMENT OF THE SWMM INFILTRATION AND RECEIVING WATER BODY ROUTINES 

6.1 General 

The terms of reference of this study required an assessment of 
two additional SWMM routines, INFIL and RECEIV. It was not possible to 
assess the INFIL routine because of data limitations. RECEIV was assessed 
only on the basis of a literature review, following a small number of 
simulations. 

The INFIL subroutine, which forms part of the TRANSPORT block, 
allows the user to model flows due to infiltration into the main conduits 
in a sewer system. Existing data indicate that rates of infiltration 
are highly site specific and, therefore, local measurements may often 
be required. If the importance of infiltration during particular 
circumstances, such as a high storm intensity, is minimal then a 
crude overall estimation of infiltration may be sufficient for modelling 
purposes. In other cases, a more rigorous investigation possibly involving 
measuring programs will be necessary. 

The receiving water model allows for the simulation of quantity 
and quality effects in rivers and estuaries, as well as in unstratified 
lakes and reservoirs. Thus, the effects of storm water discharges and 
combined sewer overflows on the receiving water may be broadly assessed. 
A versatile quality subroutine enables the improvements resulting from 
different combinations of storage and treatment to be rapidly investigated. 
Consequently, receiving water body simulations constitute an important 
factor in cost-benefit analyses for various storm water management 
alternatives. 

6.2 Inf i 1 tration 

Flows resulting from infiltration and inflow reduce the capacity 
of sewer systems and treatment facilities for handling domestic and 
industrial wastewaters. Frequently, infiltration/inflow increase water 
pollution, causing health hazards when untreated wastewaters overflow 
directly to the receiving water. Infiltration may occur due to defective 
pipes, through pipe joints, connections or poorly maintained manhole 
walls. Inflow, which is a separate phenomenon, may enter the system 
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via roof leaders, cellar, yard and area drains, foundation drains and 
cross connections. Additional sources may include catch basins, manhole 
covers, swamp drainage and storm and wash waters. Inflow is accounted 
for in the RUNOFF block, while the infiltration is determined separately 
in the INFIL routine, which forms part of the TRANSPORT block. 

Flows resulting from miscellaneous sources, comprising the 
base dry weather, infiltration, antecedent precipitation, frozen residual 
moisture and high groundwater should be accounted for, where appropriate. 
Initial wortc carried out by Metcalf and Eddy [2j indicated that it Is 
not possible to develop a widely applicable general predictive equation 
for infiltration because particular localized conditions exert a predominant 
influence on this process. Thus, the availability of local data describ- 
ing the contribution of the infiltration to measured flows should take 
precedence over generalized estimates. The INFIL subroutine may be used 
in conjunction with local data or, if this is limited, with an overall 
estimate of the infiltration rate together with information on the total 
pipe length, diameter and number of joints. The contribution of antecedent 
precipitation to infiltration is assumed to depend on the magnitude of 
rainfall and the time elapsed since that rainfall. Local data can be 
used to generate coefficients for a linear regression equation describing 
infiltration due to antecedent rainfall: 

RINFIL = ALF + ALFO x RNO + ALFl x RNI + + ALF9 x RN9 (U 

where ALFN = coefficients determined from regression analysis for 
N = I, 2 ... 9 days antecedent precipitation 

RNN = precipitation on N days prior to estimate (inches) 

After nine days, infiltration resulting from previous storms is considered 
insignificant. Residual moisture in the form of snow or Ice is 
considered to contribute to infiltration when the degree day index 
exceeds 750. The model adds infiltration fron this source, where appro- 
priate, according to a sinusoidal relationship over the melt period. 

In a situation where the groundwater table is above the pipe 
invert, groundwater infiltration Is assumed to supercede infiltration 
from other sources. Estimation of the magnitude of the groundwater 
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contribution is then based on a regression expression involving the 
groundwater head above the sewer invert and measured infiltration rates. 
In the absence of significant data, this source has to be assessed 
in a less rigorous manner (i.e. by calibration). 

Data from several Canadian studies were collected and recommenda- 
tions from the literature reviewed. There was a very wide variation of 
measured and design infiltration rates (Table 11). These data were not 
suitable for establishing regression equations similar to those outlined 
above because the typically very high values indicated the predominance 
of inflow rather than infiltration. Infiltration in combined systems 
is not normally significant when compared with the greater flow volumes 
associated with storm runoff. Consequently, for combined systems infiltra- 
tion can normally be ignored, or the phenomenon could be included as a 
constant rate. Since sanitary sewers with extraneous flows are operating 
as partially combined sewers, mathematical models such as SWMM, which 
simulate both dry weather flow and storm water runoff, may be important 
tools for assessment of system performance. Flow measurements are 
required to enable accurate calibration of the SWMM model to the 
sanitary system under study. Calibration may be achieved by adjusting 
areas with contributing runoff until good agreement is achieved between 
predicted and recorded flows. Design storms may be simulated with the 
calibrated model to identify specific sanitary sewer back-up and other 
system inadequacies and to test the effectiveness of proposed relief 
measures . 

A modelling approach was used in San Francisco to identify 
sources of infiltration/inflow, to assess performance of existing sanitary 
collection and treatment systems, to determine overflow characteristics, 
and to develop cost-effective alternatives for system control [77]. 

6. 3 Receiving Water Body Model 

In order to assess the need for and the efficacy of various storm 
water management and pollution abatement schemes, it is necessary to study 
the effects of pollutant discharges in the receiving water. 
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TABLE 1 1 



SOME MEASURED AND DESIGN INFILTRATION RATES 



Source 

Gore and Storrie, York- 
Pickering Sewer Project, Ont 

J.F. MacLaren Ltd. 
Halifax, N.S. 

Proctor & Redfern Ltd. 
Sault Ste. Marie, Ontario 

J.F. MacLaren Ltd. 
Atlantic Prov. Studies 

Canadian British Ltd. 
Halifax, N.S. 

Phi 1 1 ips £ Roberts Ltd. 
Burl i ngton , Ontar io 

Proctor & Redfern Ltd. 

J.F. MacLaren Ltd. 
Burlington, Ontario 

J.F. MacLaren Ltd. 
Ottawa, Ontario 

Black & Veatch 
Mission Twp., Kansas 
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Pollutant concentrations resulting from an untreated storm water 
discharge or a combined sewer overflow must be known or predicted, so that 
any deleterious effects can be quantified and compared with standards 
or regulations. tf violations are evident, the water quality improvements 
resulting from various storage and treatment policies should be investigated 
and related to the costs of treatment. Such a procedure enables planners 
and decision makers to determine the cost effectiveness of various 

management schemes. 1 

I 

6.3.1 General discussion of receiving water quality models 

Traditionally, water quality models have been limited to the 
solution of the dissolved oxygen balance in rivers subject to constant 
waste discharges. The earliest models gave only steady state solutions 
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for the dissolved oxygen (DO) profile along the river and served to identify 
the location of the maximum DO sag [78]. More recently, dynamic models have 
been used to define the temporal and spatial variations in the DO along 
the river, such as in the Thames River and Saint John River basin studies 
in Canada [21, 79J • €stuarial systems are the most complex aquatic 
systems and, consequently, the procedures for modelling these are 
generally more complex than for simple nontidal river systems. One of 
the most extensively applied estuarial models is that of the Delaware 
River Estuary [80] . This was proposed as a tool to describe the seasonal 
steady state DO concentrations resulting from alternative pollution 
control policies. The model requires the system to be divided into a 
series of perfectly mixed finite segments within which the flow is uni- 
directional. The original model has been expanded to include the predic- 
tion of the fate of nitrogens in the estuary I8l]. The model has also 
been applied for the Potomac River [82] and for Hillsborough Bay in 
Florida [83]- The M.I.T. one-dimensional dynamic network model was 
applied in a study of the St. Lawrence River [Bk] . 

The SWMM receiving water body nrodel is a particularly versatile 
dynamic simulation model. In the RECEIV subroutine of the SWMM, the 
receiving water is discretized by dividing the system into a series of 
one-dimensional channels. Where the actual channel broadens into an 
embayment or estuary, the system may be represented by a network of 
one-dimensional elements, connecting nodal points, approximating the 
volumetric and flow effects in the natural system. Thus, a continuous 
area of water such as an estuary or lake may be simulated. 

The RECEIV subroutine requires the input of the average channel 
dimensions and roughness, as well as nodal characteristics including 
surface area, depth and surface elevation. The surface area of a node 
in a simple river simulation is the sum of the two adjacent channel 
half areas. In a more complicated triangular network, the nodal areas 
are determined manually by the Thiessen polygon method, or the model 
subroutine TRIAN may be used to generate approximate junction areas 
and channel dimensions based upon a minimum of input data. The hydro- 
dynamic information generated in the QUANTITY subroutine forms the basis 
for dUALITV simulations. Variable water withdrawals and additions, 
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waste loadmg inputs and a variety of tidal and nontidal boundary conditions 
may be modelled and the distribution of conservative and nonconservat i ve 
constituents simulated. An explicit solution to the equations of 
continuity and momentum is applied at each time step. Constituent 
concentrations are established from the conservation of mass equation. 

Earlier versions of the model were applied by WRE for the 
simulation of water quality in the Sacremento-San Joaquin Delta [85] 
and for the investigation of potential problems associated with landfilling 
in the San Francisco Bay. The SWMM and associated models have proved to 
operate well in the simulation of the mass transport of conservative 
substances in situations where advection is the major transport mechanism 
[86, 87]. One of the limitations of the SWMM approach is the absence 
of a diffusion term in the mass transport equations. This was emphasized 
in a recent application of SWMM to a network of Finger Fill Canals in 
Florida [88]. For this particular configuration the model failed to 
correctly predict the mass transport in canal systems; the simulated 
transport of a dye cloud lagged far behind the observed movement. The 
authors also applied the Columbia River Model [89,90] which, like 
SWMM, was developed from the original WRE node and channel model. 
However, this model includes a diffusion term in the transport equations. 
The results indicated that "in situations where the net advective transport 
is small, the SWMM is not capable of simulating mass transport". The 
Columbia River model, with an appropriate diffusion constant, was shown 
to adequately simulate mass transfer in these types of systems. 

A good example of a sophisticated two dimensional receiving 
water model is that developed by the Rand Corporation for the simulation 
of coastal and estuarine systems [91]. This model solves the basic 
differential equations of motion and continuity by a time centred diffe- 
rence scheme and an integrated quality simulation routine, based on the 
mass balance equation that solves the quality equations by an alternative 
direction implicit-explicit technique. This model, which uses a grid 
system to represent the water body, has been applied for Jamaica Bay, 
New York City [92]. The model was modified and applied for the simula- 
tion of water quality in Hamilton Harbour, Ontario [93]- 
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Another two-dimensional simulation model was developed originally 
for the U.S. Atomic Energy Commission for predicting the behaviour of 
thermal plumes [9'*]- This model (HYETA) , which uses a horizontal square 
grid system of discretization similar to the Rand model, has been adapted 
for the simulation of water quality constituents. HYETA accounts for 
advective and diffusive transport. A comparison of the results of SWMM 
and HYETA for a hypothetical receiving water is discussed in the next 
sect ion . 

The computational costs and data requirements of different 
models must be balanced against the degree of accuracy required, which 
is usually determined by the scale and stage of development (i.e. screening, 
planning or design) of the project. 

6.3.2 Applications of the model 

Several test applications were carried out in order to test 
the response of the RECEIV block. Figure 52 shows the computed suspended 
solids profiles along a simplified river system. The movement of the 
concentration peak downstream is simply related to the average flow 
rate (i.e. there does not appear to be significant "numerical dispersion"). 
A typical lakefront sewer outfall is shown in Figure 53 in two alternative 
locations (inside and outside the breakwater). The changes in the 
level of pollution of the receiving water resulting from the relocation 
of the sewer outfall outside the breakwater were investigated. Hypothetical 
pol lutographs based on typical recorded storm water pollutant concentrations 
were used in these simulations. The plots shown on Figure 5^+ indicate 
the coliform concentrations at the nodes inside and outside the breakwater 
for each alternative. A coliform decay rate coefficient of 2 day was 
used. The plots show rather qualitatively that discharge into the deeper 
waters outside the breakwater using the greater dilution available is 
preferable to discharge inside the breakwater (assuming a water quality 
standard of about 10^ MPN/100 rr 1 ) . 

For the future applications of the model in practical situa- 
tions, it will be necessary to assess the performance of the model compared 
with measurements and with the results of more sophisticated models. 
It was possible to compare the SWMM and the HYETA models within the 
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scope of the present study, 'he comparison was based upon a hypothetical 
situation in which a river highly polluted by suspended solids (200 mg/1) 
flows into a semi-enclosed embayment with a small harbour at one side 
of the river mouth (Figure 55). The depth varies from three feet 
at the river mouth to six feet at the open end of the system. The 
distribution of suspended solids In the system predicted by the two models 
after 12 hours of simulation is indicated on Figure 55- While the two 
simulations are reasonably in agreement, two dissimilarities are noteworthy; 
SWMM appears to cause more pollutant to be advected into the area behind 
the breakwater, while HYETA predicts a more widespread dispersion through- 
out the whole area. Since verification of these findings is not possible, 
definite conclusions cannot be drawn. Although the fundamental representa- 
tion of the flow field and the inclusion of the diffusive term in the 
HYETA model is probably more accurate, the results of the SWMM indicate 
that it may be used for comparison of alternatives at a screening of 
planning stage. 

6.3.3 Remarks concerning the KECEIV model 

Several ommissions in the Users' Manual [71] became apparent 

during the initial applications of the model: 

(a) The initial tirre for start of hydrograph input from 
cards (TZERO) is controlled by the day cycle on which 
printed output begins. A discharge hydrograph will not 
be read until there is a one day cycle of printed output. 
This is not stated in the input description. 

(b) It is necessary to have at least one stage vs time plot 
when using the tidal boundary condition. This can 
apply for any node in the network (i.e. NPLT > 1), 

(c) The moael computes the decay of any nonconservat i ve 
constituents as a first order process, analogous to the 
biological oxidation of BOO: 

^ = DECAY X C (2) 

dt 
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where C = the concentration of the constituent 
at time t 

DECAY = the first order rate constant 

Any nonconservat i ve pollutant is assigned an equivalent 
conservative constituent by the model. This is 
analogous to the adding of dissolved oxygen (DO) when 
BOD is model led. Of course, when there is zero DO the 
BOD ceases to decay. DO is replaced by reaeration 
which is controlled by the rate constant REAER in 
the model . 

in order to ensure that nonconservat ive pollutants 
{such as coliforms) other than BOD can decay, it is 
necessary to ensure that there is a sufficient supply of 
the equivalent conservative constituents (analogous to 
do). The most effective method for ensuring this was 
found to be the specification of a value for REAER for the 
nonconservat ive constituent, equal to the decay constant, 
DECAY. 

These findings were discussed with representatives of the 
University of Florida and it was understood that subse- 
quent editions of the User's Manual will be modified 
accord ing 1 y , 

The minimum time step possible for a stable solution in 
the QUANTITY subroutine is given by the Courant condition. 

At = Ax 111: 



Hence the minimum value of Ax/ /gh determines the time 
step for hydraulic computations. However, it is normal 
to reduce this time step by multiplying by a safety factor 
of 0.75 to ensure that the solution remains stable. A 
time step of 25 seconds was used in the lake simulation 
described in the previous section. The tidal boundary 
condition is the most appropriate when modelling only 
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part of a large body of water. In the lake example, a 
tidal fluctuation with zero amplitude proved to be the 
most flexible boundary condition. 

6.h Conclusions 

(a) The receiving water body model of the 5WMM does not 

include a diffusion term and should not be applied in 

circumstances in which advection is not the major transport 

mechanism. 

I 

(b) The testing of the model and limited comparison with more 

sophisticated models indicated that, when carefully 

applied, the SWMM may be used as a tool for the assessment 

of the relative effects of various pollution abatement measures 

at a planning stage. 

(c) Several additional user requirements were identified and 
these are being included in updated versions of the model. 
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7. STORAGE/TREATMENT ROUTINES 
7. 1 General 1 

Equations describing the complex physical, chemical and biological 
processes associated with the various storage and treatment units are 
usuaily highly empirical and based upon limited data; consequently, most 
available models are planning aids rather than design tools. A planning 
model should be capable of predicting the water quality improvements result- 
ing from different levels of storage and treatment of a contaminated dis- 
charge and should relate a cost to different degrees of abatement. The 
decision maker should thus be able to select the optimal area for a detailed 
cost-benefit analysis and subsequent capital investment. 

In this chapter, the basic assumptions of the SWMM storage 
and treatment models are discussed and some inconsistencies are rectified. 
A cost estimation subroutine has been developed for Canadian conditions. 
The passage of a particular storm flow hydrograph and pollutograph is 
traced through a hypothetical treatment plant and a set of the results 
of the cost estimation routine is presented. 

7. 2 Storage/Treatment Options and Basic Assumptions 

7.2.1 Storage treatment options 

Storage of discharges, together with the resulting sedimentation 
occurring during storage, with subsequent release into the treatment system 
appears to be the most cost effective method available for reducing storm 
water related pollution [2]. Storage facilities may be constructed upstream 
of the treatment unit in order to provide flow equalization and dilution 
of the typically high "first flush" pollutant concentrations. Such "in- 
line" storage may be available in the sewer system, particularly in systems 
with flat sewer grades, and may be used through appropriate flow constriction 
upstream of overflow points, or may be provided immediately ahead of the 
treatment facilities. Alternatively, storage may be provided "off-line", 
as a means to retain contaminated flows until spare capacity for treatment 
becomes available. "Off-line" storage facilities may tal<e the form of 
conventional sedimentation tanks, underground and void space storage, or 
underwater bag storage. 
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The treatment processes currently employed in the SWMM for 
reduction of pollution from storm sewer discharges and combined sewer 
overflows are briefly described in the following; 

(1) Bar Racks - The main function of bar racks is to trap and 
renove large objects, which might subsequently damage 
downstream equipment, from the waste stream. Bar 

racks are generally installed for all storage, treatment 
and pumping facilities. A typical bar layout provides 
one inch clear openings spaced two inches between centres. 

(2) Swirl Concentrator - The swirl concentrator is primarily 
used for the rapid removal of settleable solids. Pre- 
liminary investigations indicate that this equipment is 
particularly appropriate for the treatment of storm flows, 
responding well to fluctuating flows [95]. The removal 
efficiency improves as the ratio of storm flow to dry 
weather flow decreases. The swirl concentrator may also be 
used at combined sewer outfalls to allow a large portion 

of the solids to be intercepted and delivered to the 
interceptor system. 

(3) Fine Screens - Fine screens were often considered as an 
alternative to sedimentation and have often been used in 
industrial applications. The screen wire may vary from 

6-mesh to 60-mesh, depending on the characteristics of the waste; 
a 50-mesh screen is most common. Problems occur in 
operation when the screens, which may be rotary or drum 
types, become clogged with oil or grease and cannot be 
cleared by normal backwashing. Fine screens generally 
renov; about 38 and 1 6 percent suspended solids and 
BOD, respectively, although the efficiency of a particular 
'jCreoTi i -J dependent on the size of wire mesh and the 
buildup of smaller particlf, on the screen [96]- 
Ih) Sr;r] int-n ta t i on - Sfrd i mc^n t a t i rjn is the mf)s t established and 
w i de'.prfod mf-thod of ri^ductirin tyf sol iris and associated 
BOLi frrjin v;'"js I fWrj t o r s . Cnn vcn t i rjn.i I dcsiqn requires detention 
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times of one to two hours. The removal efficienty is 
dependent upon the residence time, which may vary consider- 
ably during a storm, since in most cases the sedimentation 
process would take place in the storage facility. Removal 
efficiencies in conventional units are approximately 60 
and 30 percent for suspended solids and BOD, respectively, 
for a residence time of two hours. 

(5) Dissolved Air Flotation - Dissolved air flotation has 
several advantages over conventional sedimentation. High 
overflow rates may be used because of the greater effective 
difference in specific gravity between the water and 
combined (air bubble and particle) particle, than that in 
normal sedimentation. High overflow rates allow short 
detention times, and the upward separation and subsequent 
skimming allows for removal of particles less dense than 
water, in addition to surface oil and grease. The character- 
istics of the individual waste stream greatly affect achiev- 
able removal efficiencies and, consequently, pilot plant 
studies will usually be required before this process is 

i nsta I led . 

(6) Micros tra i ners - This unit is similar to the drum screen 
(fine screen) device in operation. The drum rotates about 

a horizontal axis allowing gravity filtration, with continuous 
backwashing at the top of the drum. As with fine screens, 
the removal efficiencies are dependent upon the type of 
media used and the concentration and characteristics of the 
solids in the waste stream. 

(7) High Rate Fi I ters - Filtration removes finer particles than 
those removed in screening processes. In conventional plants 
filtration is usually regarded as a polishing treatment. 
High rate filters generally consist of dual media or mixed 
media. At present, operating costs are high and frequent 
backwashing is required to prevent clogging, but land 
requirements are not extensive and the process is easily 
automated . 
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(8) Biological Treatment - For successful biological treatment 
of storm related wastes, it is necessary to maintain an 
active microbial biomass to use the dissolved organic and 
nonorganic matter in the waste flow. This may be achieved 
either by operating the wet weather unit In parallel with 
the conventional dry weather plant or by using the dry 
weather plant for dry and wet weather flows. The latter 
would be achieved using treatment units able to operate 
efficiently over a wide range of flows and substrate concen- 
trations. Units such as tricl<ling filters tend to operate 
well under conditions of fluctuating waste strength but 
upstream flow equalization is probably essential to prevent 
excessive hydraulic loads. In addition, all biological 
treatment processes necessitate downstream units for 
settling out of cell tissue, and subsequent sludge disposal. 

(9) Effluent Screens - These screens serve only to remove floatable 
and unsightly material that might remain after low levels of 
treatment. Their function is essentially aesthetic 
improvement of the discharge, there being no discernible 
associated improvement in water quality. 

(10) High Rate Disinfection - Disinfection of waste discharges 
is generally achieved by chlorination. This results In a 
reduction in the coliform content of the waste In the order 
of 99.9^- The efficacy of the disinfection depends on the 
free chlorine residual concentration and the time of contact. 
High rate processes usually operate with a 15 minute detention 
period at the design flow. 

7.2.2 Discussion of basic assumptions 

The storage subroutine in SWMM, ST0RA6, allows for in-line and 
off-line storage basins. In-line storage implies all flows from the sewer 
system flow through the storage basin prior to treatment, whereas off-line 
storage Is only used when the treatment capacity is exceeded. Two Idealized 
flow regimes are allowed: plug flow and perfect mixing. Hydraul leal ly , 
the choice of regime In the model is not important, but when pollutant 



concentrations are of interest, the flow regime is of great importance. 
Plug flow is more appropriate for in-line storage and perfect mixing for 
off-1 ine faci 1 i ties. 

The STORAG subroutine recycles flows from the storage basin that 
are in excess of the treatment capacity back to the storage basin. This 
recycling is not always practical; in general, flows in excess of the treat- 
ment capacity are allowed to overflow once the capacity of the storage 
basin is reached. The removal of suspended solids in off-line storage basins 
is considered to occur subject to an exponential decay equation. The 
assumption of the completely mixed basin is not considered appropriate for 
storm water storage facilities and sedimentation would probably be better 
approximated using a mass balance approach. 

The treatment subroutine, TREAT includes all the treatment processes 
discussed in the previous section. The flow chart in Figure 56 indicates 
the possible sequential series of storage and treatment options. Not all 
of the "strings" of treatment processes are technically feasible or econo- 
mically realistic. However, a safety mechanism is available in this 
surboutine whereby the program will terminate upon location of an unfeasible 
"string" of options and will identify economically unsound combinations 
with a warning message. The possible combinations presently available 
are indicated in Figure 57. 

The cost of a specific combination of storage and treatment 
is a function of design flow rate. Unless this is specified by the user, 
the model assesses the peak flow rate from the inlet hydrograph and sizes 
the treatment units accordingly. An evaluation of the range of treatment 
efficiencies possible for a particular storm discharge is obtained by 
specification of a series of design flow rates. 

The cost estimate for a particular combination of storage and 
treatment facilities is determined in the cost subroutine, TRCOST. The 
costs are a function of the design flow rate because, in the model, the 
unit size determines the capital outlay required. 

7.3 Modifications 

In the original storage subroutine, flows released from the 
storage basin, controlled by weirs, constant speed pumps or orifices, could 
either flow directly to the treatment facilities or to a by-pass. The 
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model considered flows discharged from storage, in excess of the treatment 
capability, to be recirculated to the storage basin. This process was 
modified to cause flows in excess of the treatment capacity to be routed 
to the bypass stream. 

The original model removed suspended solids In the storage 
basin according to an exponential decay equation based upon the assumption 
of perfect mixing. The removal process was altered to reflect a more 
empirical basis [97J and this was incorporated in an overall mass balance 
equation: 

(Inflow SS) = {Outflow SS) + (SS removed in storage basin)^'^ (1) 
+ (SS in suspension in storage basin) 

This approach represents actual flow conditions better than the original 
idealized formulations. 

In the original treatment subroutine there was no consideration 
of underflow (the outflow stream containing high solids concentrations due 
to separation in rotary flow field) in the swirl concentrator, leading 
to a discrepancy in the mass balance across this unit. This is not realistic 
and, consequently, an underflow rate equal to l/70th of the inflow rate 
(obtained from empirical data) has been incorporated for this unit. 
Also in the treatment section, the suspended solids removal efficiency 
was raised from 22 to 35 percent to reflect recent reported improvements 
in performance [98]. This unit will, therefore, appear more cost effective 
in the revised model. ' 

The cost estimates in the model for the individual treatment 
units are based on the general relationship 



S = AQ^F 



I 



where A is a base cost factor, Q is the design flow rate for the unit (in 
mgd), B is an indicator of the economy of scale and F is a correction factor 
for the specific time and locality. Cost data for storm water facilities 
in North America are rather limited and a reliable statistically signi- 
ficant data base is not yet available. At present, much of the activity 



where settling efficiency is given by ^-=0.82e 
and Q is the overflow rate in US gpd/ft 
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In storm water storage and treatment is initiated by governmental agencies 
and, consequently, the available figures may not reflect all the costs 
involved in construction. Cost curves and equations were developed 
using the publications of Environment Canada, Provincial Ministries of 
Environment and the U.S. EPA, as well as bid construction costs and 
equipment manufacturers' estimates. The costs incorporated in the revised 
model, where relevant, include estimates for electrical piping, roads and 
other related capital expenditures. The treatment plants assessed in this 
study ranged in capacities from 2.2 to k^.k cfs. The model can accommodate 
design flows up to 770 cfs. Obviously, there will be a considerable degree 
of uncertainty in cost projections for the higher flow range, but generally 
the relative costs of different treatment combinations should remain 
realistic. The present model uses the ENR Index, published in Engineering 
News Record, as a means of standardizing estimates. This is considered 
to be the best available index, although the ENR may be somewhat limited 
in the area of municipal wastewater treatment. The Water Quality Office - 
Sewage Treatment Plant (WQO-STP) Index has proved successful in adjusting 
treatment plant costs to a common time period, but at present insufficient 
data exist to perform a comparison. Table 12 summarizes the capital cost 
functions included In the revised model for the Individual treatment 
options. These were adjusted for Toronto, March 1975, the ENR Index being 
2033. 

7 .^ Testing Storage/Treatment Options 

Each storage and treatment unit was tested individually and in 
series with other units following the modifications to the subroutine. 
This involved running the STORAG/TREAT block with a test hydrograph 
and pollutograph for each process. Some of the untreated and treated 
pel lutographs for the various units are compared in Figures 58 to 61. 

A full scale test of the modified STORAG/TREAT and TRCOST 
routines was conducted using the storm of September 5, 1971 with data for 
the Bannatyne combined sewer system. Figure 62 illustrates the effect 
of the storage basin in attenuating the inflow hydrograph (which is output 
generated by the TRANSPORT block). The treatment plant capacity was 
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TABLE 12. CAPITAL COST FUNCTIONS 



(ENR " 2033. F = 1.0) 



Treatment Options 

1. Bar Racks 

2. Swirl Concentrator 

3. Inlet Pumping 

4. Dissolved Air Flotation 

5. Fine Screens 

6. Sedimentation Tanks 

7. Microstrainer 

8. High Rate Filters 

9. Biological Treatment 

10. Effluent Screens 

11 . Outlet Pumping 

12. Chlorine Contact Tank 

13. High Rate Disinfection 



Cost Functions 
40400 Q*^V 
3330 DIA.F 
58500 Q-^^F 



.43r 



110000 [Q(l + RECIRC)]*^''F, Q^lOmgd* 
35600 [Q(l + RECIRC)]-^^F, Q^lOmgd 

1857QF 

5840000 {^/R)-^^F 

67000 Q'^^F, Q^20 mgd 
27275 Q-'^^F, Q;'20 mgd 

56000 Q-^°F 

1400000Q-^^F 

11200 Q-^^F 

58500 Q-^^F 

28100 Q-^^F 

30700^ '^F 



where 



Q = design flow rate (mgd) 

F = correction factor 

DIA = unit diameter 

R = overflow rate (gpd/ft^) 

RECIRC = recirculated flow 



All volumes in U.S. gallons, 



Cont'd 
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TABLE 12. (CONT'D) 



USEPA - SWMH CAPITAL COST FUNCTIONS 
FOR COMPARISON 
(ENR = 2033. F = 1) 



Treatment Options 



USEPA 

SWMM Cost Function 



1 . Bar Racks 



2. Swirl Concent. 

3. Inlet Pumping 



4. Dissolved Air 

Flotation 

5. Fine Screens 

6. Sedimentation 

Tanks 



7. Microstrainers 

8. Hiqh Rate 

Filters 



22363 + 20330 Q"^^^ F 

Q f TOO mgd 
22363 + 3499 Q F 

Q 2 100 mgd 

3140 Dia. F. 

38680 Q-^^ F 

Q 5 20 mgd 
24755 Q-^^ F 

20 * Q < 100 mgd 
7134 Q F 

q =»■ 100 mgd 

45016 Q F 



18570 q F 

87419 (700q/R) 
q 5 100 mgd 



.91 



874 q (70000/R)-^^ F 
q > 100 mgd 

14699 QF 

163332 Q*^^ F 



Cont'd 
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TABLE 12. (CONT'D) 

USEPA - SWMM CAPITAL COST FUNCTIONS 
FOR COMPARISON 
(ENR = 2033. F = 1) 



Treatment Options 



USEPA 

SWMM Cost Function 



9. Biological 
Treatment 

10. Effluent Screen 



11 . Outlet Pumping 

12. Contact Tank 

13. High Rate 

Disinfection 



79868 Q F 

i 

13763 Q-^^^ F 
Q <■ 100 mgd 

2759 Q F 

Q ? 1 00 mgd 

37306 Q-^" F 

1295 Q F 

3070 Q F 
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specified as S'* cfs. This involves a small overflow over a period of 
about 10 minutes as can be seen from the flat section of the hydrograph 
computed at the inflow the the treatment units. There is also some initial 
removal and dilution of the three pollutants simulated in the storage 
basin. The increases in the concentration of these pollutants towards 
the end of the simulation period are explained by the diminishing 
dilution of dry weather flow as the magnitude of storm flow declines. 

The subsequent performance of a treatment process consisting of 
swirl concentration, microstrai ning and chlorine disinfection Is 
illustrated in Figure 63- The simulated concentration of BOD in the 
discharge stream is always less than 10 mg/1 , compared with peak concentra- 
tions in excess of 200 mg/i in the Inflow to storage. The reduction 
in suspended solids is not so dramatic, with a peak concentration of about 
120 mg/1 in the discharge. However, in general, these results appear 
reasonable. The average coliform concentration in the discharge is 
reported as 3-^ x 10 MPN/100 ml (not plotted). Thus, even after dilution 
in the stream, it Is possible that the effluent from this combination of 
treatment processes would not satisfy primary contact bacteriological 
standards in the receiving water. j 

A summary of the capital, annual operating and storm event 
costs estimated in the COST subroutine for the treatment units specified 
in the example is shown in Table 13- 

:^,.§ Concl us ions 

(a) Modifications have been made in the STORAG and TREAT 
subroutines to better reflect current data on unit 
performance. | 

(b) Each treatment unit has been tested in order to ensure 
the model is fully operational and that variations in the 
input stream are reflected reasonably in the output 
stream. The model may be used for the comparison of 
alternative pollution abatement facilities. 

(c) The TRC05T subroutine has been modified in order to 
better reflect Canadian prices. At present this can give 

the planner only a genera! view of the relative costs involved 
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TABLE 13. SUMMARY OF TREATMENT COSTS 



INPUT 



ASSUMED FUTURE ENGINEERING NEWS RECORD INDICES 
CONSTRUCTION - 20 CITY AVERAGE 



3% 



YEAR 


ENR INDEX 


1970 


1314 


1971 


1346 


1972 


1378 


1973 


1410 


1974 


1442 


1975 


1474 


1976 


1506 


1977 


1538 


1978 


1570 


1979 


1602 


1980 


1634 



COST PARAMETERS 



INTEREST RATE 


= 


7.00 Percent 


AMORTIZATION 1 


'ERIOD = 


25 Years 


CAP. RECOVERY 


FACTOR = 


0.0858 


YEAR OF SIMULATION = 


1974 


SITE LOCATION 


FACTOR = 


1.0000 


UNIT COSTS . . 






LAND 


1000.00 $/acre 




POWER 


0.020 $/kwh 




CHLORINE = 


0.200 $/lb 




POLYMERS = 


1.200 $/lb 




ALUM 


0.03 $/lb 





OUTPUT 



TOTAL LAND REQUIREMENT = 0.90 ACRES 



TREATMENT 

BAR RAK & SW.C. 

NO INLET PUMPING 
BYPASS LEVEL 3 
MICROSTRAINERS 

NO EFFL. SCREENS 

NO OUTLET PUMPS 
CONTACT TANK 



LEVEL 

1 
2 
3 
4 
5 
6 
7 

SUB-TOTAL 

TOTAL 

TOTAL PER 
TRIE ACRE 



CAPITAL COSTS 



INSTAL 

213804. 

0. 

0. 
364931 . 

0. 

0. 
246761 . 



LAND 

72. 

0. 

0. 
774. 

0. 

0. 
141. 



$825496. $988. 
$826484. 

$1524. 



ANNUAL COSTS 



INSTAL 

13764. 

0. 

0. 
31315. 

0. 

0. 
21175. 

$66254. 



LAND MIN MAINT 



5. 

0* 

0. 
54. 

0. 

0. 
10. 

$69. 

$801 60. 

$ 148. 



1604 



7299 



4935 

$13838. 



STORM EVENT COSTS 



CHLORINE 


CHEM 


OTHER 


0. 


0. 


28. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


0. 


34. 


0. 


0. 


0. 


0. 


0. 


0. 


31. 


0. 


15. 



$31. $ 0. $ 76. 

$107. 



$ 0. 



with various treatment alternatives. However, as more 
statistically useful data becomes available following the 
installation of new facilities and pilot plants, the 
treatment and cost estimate models should be updated to 
provide more accurate indications of both efficiency and 
cost of treatment. 
(d) The time step employed in the treatment computations is 
usually the same as that used in the TRANSPORT blocit 
(typically five minutes). Since the equations in the TREAT 
subroutine are based mainly on data obtained from units 
operating at steady state, it is considered that a longer 
time step (say 30 minutes) could be used in treatment 
computations without a significant loss of accuracy. 
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8. SNOWMELT QUANTITY AND QUALITY 
8. I General 

Snowmelt can be accompanied by intense rain storms in spring 
or occasionally in late fall. This may lead to significant flooding in 
urban areas, even in the case of small basins. In urban areas of Canada, 
water quality during winter may be affected by application of chemical 
de-icing agents. In some cases, the snowpack acts to retain various 
pollutants which are then released to the receiving water body during 
the melt period. The prediction of the time of release of melt water 
in some cases is complicated by the fact that chemical additives 
and the mechanical removal of snow affect the natural snowmelt processes. 
The amount and distribution of snow on streets and other impervious and 
pervious surfaces must be known or estimated prior to an assessment of 
melt quanti ties. 

The basic problem is, therefore, to simulate the quantity and 
quality of runoff water resulting from snowmelt or a combination of 
rainfall and snowmelt. Complete reviews of literature concerning the 
quantity and quality of snowmelt are given in Volume II. In the 
following sections the results and conclusions of the literature surveys 
are discussed and the model developed for integration with the SWMM is 
described. 

8. 2 Selection of a Snowmelt Quantity Model 

8.2.1 Summary of the literature review 

The available snowmelt models can be broadly classified as: 

(a) detailed heat budget models requiring extensive 
meteorological input data; 

(b) regression type models of an "intermediate" 

nature accounting for the most important meteorological 

parameters; and, 

(c) empirical models such as the degree-day method. 

The detailed models are based on the following generalized heat budget 
equat ions: 
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H ,^ = H - H 

melt snow cc 



H =H +H, +H +H ^ + H+H C^l 

snow swr Iwr ca cond g p 



where; H = short wave radiation absorbed by the snow pack 
swr 

u 

Iwr = net long wave radiation exchange with the snow pack 

H = convective heat transfer from the air 
ca 

H , = heat of vaporization released by condensate on the snow 
cond c 

surface 

H = heat conducted from the ground 

H = heat content of precipitation 
P 

H = change in heat storage In the snow pack 
snow 

H , = heat available to melt a given amount of snow 
mel t ^ 

H = "cold content" of the snow pack 
cc 

When a snow pack is below the freezing point it has a negative heat content. 

The heat required to raise the temperature of the snow pack to 0°C can be 

computed from equation (I). Any additional heat causes the snow to melt. 

The most important factors governing the heat transfer processes involved 

in equation (I) are air temperature and net radiation. 

Several existing detail led mathematical computer models are based 

on various forms of equation (l). These include Amorocho's model, the 

Stanford Watershed model and the Hydrocomp simulation program [99, 100, 101] 

(These snow accumulation and melt models are described in Volume 11.) The 

various mathematical relationships employed in most of the heat budget 

models are based on the work of the U.S. Army Corps of Engineers [102]. 

Sophisticated models generally require a considerable number of input 

parameters, including: 

- short wave radiation; 

- air temperature; 

- dew point temperature; 

- wind velocity; 

- precipitation; and, 

- cloud cover fraction. 
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The computational time interval for such models is usually limited to one 
hour because the necessary meteorological data are usually not available at 
more frequent intervals. 

Several empirical snowmelt models have been developed based 
on regression relationships. These models incorporate the most important 
meteorological parameters, but are not generally applicable over a range 
of different geographical areas. 

Because of its relative simplicity, the so-called degree-day 
method has received the widest application. A degree-day is defined as 
a deviation of one degree from a given datum temperature over a Zk hour 
period. Snowmelt is given by the following equation: 

MELT = C (T - T ) (2) 

a sn 

where: C = melt coefficient for the particular area, 
T = air temperature, 

T = base temperature for melt, 
sn 

MELT = snowmelt in inches (or centimeters) per day. 

Equation (2) may be adapted for the computation of melt during time periods 
of less than one day. 

The degree-day method gives good results in view of the simplifying 
assumptions and data limitations. For example, Anderson compared a detailed 
heat budget method with the degree-day method and concluded that both methods 
gave quite similar results. Anderson has developed a model which combines 
the advantages of a degree-day computation with the flexibility of the 
heat budget method. This model is described in more detail in Section 

8.2.2 [103, \Qk]. 

One of the major obstacles to an accurate simulation of the 
snowmelt runoff process is the need to satisfactorily account for the 
movement of melt water through the pack and for the resulting runoff. 
This problem is sometimes aggravated by unusual infiltration rates due to 
partially frozen topsoil. The literature survey indicated that, at present, 
there are insufficient data to develop a runoff nradel sophisticated enough 
to account for all these phenomena. Consequently, the SWMM runoff equations 
were modified for use in the present model. 
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8.2.2 Snowmelt quantity mode l 

The selection of an appropriate snowmelt model for integration 
with the SWMM was based on the following considerations: 

(a) A time step of one hour is sufficiently accurate for 
the computation of the total volume of melt water 
during periods of no rainfall. (Intermediate values 
required for short time steps in SWMM calculations are 
interpolated.) This accuracy is also sufficient for 
hydrograph computations during rainfall periods, since the 
rainfall intensity for the short time intervals involved 
is normally significantly greater than the snowmelt 
intensi ty . 

In general, for single event modelling, the snow pack 
distribution and physical parameters are known or may be 
determined in order to define initial conditions prior 
to a storm. The problem is then to describe the physical 
changes in the snow covered areas during the rain storm 
and the resulting effects on the runoff. 

(b) The following assumptions can also be made for single- 
event simulation: 

(i) Frozen ground remains frozen during the given event. 

(ii) The area I distribution of snow cover can be assumed 
available as an input. The areal distribution can 
be accounted for on a subcatchment basis and this 
distribution is assumed constant throughout the 
event simulation. Snow cover in the different sub- 
catchments will be depleted with time during the 
simulation according to the amounts available at the 
beginning of the event, 

(iii) It is assumed that during the simulation, the melt 
volume resulting from the application of de-icing 
salts is negligible. It is assumed that no additional 
salt is applied during the melt period. The effect 
of salt on the free water in the snowpack at the 
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beginning of the storm is accounted for in the model 
input parameters, which can be measured or calibrated 
for a particular event. 

(iv) Incoming solar radiation during a rain on snow event 
is negligible and the incoming long wave radiation 
equals the blackbody radiation at the ambient air 
temperature. 

(v) The dew point temperature and the temperature of the 
rain water are equal to the ambient air temperature. 
The temperature of the melting snow surface is 
32°F. 

(vi) For single-event simulation of shallow (nonmountain 

type) snowpacks, the negative heat component of equation 
(1) is assumed to be negligible. That is, heat applied 
to the snowpack will immediately result in some melt 
which would be either stored in the snow or become 
available for runoff. 

(vii) The air temperature record and average wind speed 
are ava i 1 ab le . 

(c) Paved areas are usually cleared of snow in urban areas. 
Consequently, it is anticipated that the runoff from 
these areas in a subsequent rainfall event would be 
greater than from the snow covered areas. The ability 
to describe bare areas as well as snow covered areas is 
discussed in the user documentation. 

The Anderson snowmelt model was selected on the basis of the above 
assumptions. During rainfall periods melt is given by the following 
equal ion: 

HELT = (T - 32) - (0.007 P AT + 7-5 )i F(u) + 0.001 17AT) 

+ 8.5 F(u) - (e - 0.18) (3) 

a 

where: T = air temperature ( F) 
a 

F(u) = wind function = .006 V&T 

e = vapour pressure of air (inches H ) 
a 9 
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P = rainfall in inches per hour 

X 

X = .000359 PA where PA is the air pressure 

in inches of Hg 

AT = time interval in hours 
V = wind speed in mph. 

The general algorithm for computing snowmelt using the above equations is 
given in Figure 6^. Air pressure and vapour pressure are calculated by 
the model using the relationships summaried in Figure 6'i. (Further 
discussion is contained in Volume II and [lO't].) The melt equation is 
linear with respect to time. Therefore, the melt and associated variables 
may be determined by interpolation at each of the computational time steps 
used in the RUNOFF block. 

Equation (3) should be used for rain on snow events or during 
cloudy periods preceding or following rainfall events. For period of 
no rainfall, the degree-day melt equation (2) is used to compute snowmelt 
by adjusting the melt coefficient to be consistent with the time 
interval for computation. A base temperature between 25 to 32°F can be 
selected, depending on the geographical location and prevailing 
meteorological conditions. A base temperature of 32 F and a degree-hour 
melt coefficient of 0.003 are assumed as default values. Methods for 
estimating the base temperature and melt coefficient can be found in 
the literature (e.g., [105]). It should be noted that use of the degree- 
day option for rain-free periods will give time averaged snowmelt rates. 

8 . 3 Snowmelt Quality Model 

8.3.1 Summary of the literature review 

The pollution of snow and melt waters was not investigated in 
detail until recently and, therefore, the literature is rather limited. 
The main efforts to date have been to measure the actual amount of various 
pollutants in snow in urban areas. The literature found describing the 
modelling of pollutants in snow or in snowmelt waters was very limited. 

The major sources of pollutants in snow are related to de- 
icing chemicals and automobile exhausts. The amount of pollution is related 
to the degree of urbanization in Table \k. Areas of medium population 
density (IOOO-5OOO per sq mile) receive the greatest amount of salting 
11061. 
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TABLE 14. ROAD SALTING RATES IN ONTARIO 



Population Density Rates of Salt Application 

(per sq mile) (lb per app. per lane-mile) 



Less than 1000 75 - 800 

From 1000 to 5000 350 - I800 

More than 5000 '+00 - 1200 



Pollutant concentrations in snow are found to be highest along the most 
heavily travelled roadways. 

Chloride concentrations measured in snow from streets in 
Canadian urban centres have been reported to range from 7 to over 15,000 
mg/I . Recent sampling of runoff during the winter season of 197'i"1975 
at the Brucewood catchment in North York has indicated chloride concentra- 
tions as high as 10,300 mg/1 [107, 108]. 

The concentration of lead in snow samples varied up to 113 mg/1 
along freeways in Ottawa. The lead concentrations in Ottawa storm sewer 
runoff were sometimes as high as 1 mg/1. Recent samples of runoff in the 
Brucewood catchment gave measured lead concentrations up to 2 mg/1 

[36, 109]. 

Generally speaking, chlorides and lead are the most serious 
pollutants associated with snow and snowmelt waters. Table 15 summarizes 
the ranges of various pollutant concentrations which have been measured: 

Detailed records of the amount of various pollutants entering 
the sewer system as a result of snowmelt do not exist. Very few 
snowmelt quality models have been documented. (Comprehensive reference 
material is reported in Volume tl). 

8.3.2 Development of snowmelt quality model 

Since specific information concerning the modelling of 
snowmelt quality was not available, it was decided that the methodology 
used in the SWMM for storm runoff quality computations should be adopted 



TABLE 15. REPORTED SNOW AND SNOWMELT WATER POLLUTANT CONCENTRATIONS 



Pol lutant 



Observed Range {mg/1) 



Chloride 
Lead 

BOD, 

b 

Suspended Solids 
Organic Nitrogen 
Ni trate-N 
Ammon ia-N 
Total Phosphate 
Sod i um 
Phenol 



7-15580 

0.02-1 13 

I. '4-8'* 

l-itOZO 

0.09-^.3 

0.01-1.5 

0.1-0.3 

0.036-3.6 

3.6-9^80 

.002-. 125 



for srtowmelt quality modelling. Chlorides and lead are model led^ tn 
addition to other pollutants normally simulated by the SWMM. The 
basic equation describing pollutant washoff is; 

AP = P (1 - e"''^^'^) I (4) 

where: P = initial surface pollutant load 
o 

b = a constant (usually = h.S) 

At = time interval 

r = runoff rate 

AP = mass of pollutant removed in At 

(These are discussed more fully in Chapter 5.) 

The initial amount of any pollutant (except chloride) on the 
^^round at the start of a storm can be computed from: 

P = P • G, • N- (5) 

o I D I 
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where: P = amount of pollutant at the start of storm 
o 

N = number of dry days since the last storm 
D 

G = total length of gutters in hundreds of feet 

P = accumulation rate of pollutant lb/100 ft of gutter/day. 

If the number of dry days exceeds the cleaning frequency, the 
surface pollutant load is modified according to the number of sweepings 
and the efficienty of sweeping. (The equations are given in Chapter 

5.) 

Since the deposition of chlorides is primarily a function of 
local street salting practices, a different approach is used for the 
estimation of the surface accunul at ion of chlorides. In developing 
the snowmelt quality model it ^as decided to supply as input to the 
model the total amount of salt applied to the catchment prior to the 
event. These data are usually available from municipal records of street 
salting events. Since different salt application rates are applied in 
different streets, the amount of salt applied to each subcatchment 
should be input. Where measurements are not available, the application 
rates given in Table 1^ can be used to obtain initial estimates of 
sal t appl icat ions . 

The total amount of salt applied to the catchment is then 
converted into chlorides by using the following relationship: 

P = p X 600 X ^+53. 6 (6) 

o , o 
cl s 

wherfe: P = total amount of chlorides i n tng at the beginning of 
cl storm, 



P = total amount of salt applied prior to the storm in 
and, 600 X ASB-^ - a unit conversion factor. 



The amount of chlorides washed off during the storm is then computed in 
the same manner as for other pollutants. 

A description of the required input data formats is given in 
the User's Manual, Volume llli 



8. 'i Integration and Verification of the Snowmelt Models 

S.*!.! Integration 

(a) The snowmelt quantity model discussed in Section 8.2.2 
(see Figure 6k) has been integrated with the SWMM 
RUNOFF block. Two equations for snowmelt have been 
included as options in the integrated model. During 
rainfall periods melt is given by equation (3). Snowmelt 
in rain-free periods can be computed from the "degree- 
hour" method, which is simply the degree-day method modified 
for use with shorter intervals, or by using equation (3). 

The principal input parameters are: 

- wind speed (mph) ; 

- ambient temperature ( F) ; 

- distribution of snow; and, 

- snow amount (in inches of water equivalent). 

The SWMM RUNOFF block has been modified to account for 
rainfall plus snowmelt runoff from completely or partially 
snow-covered pervious and impervious areas. The time lag 
to runoff, which is affected by the ability of the snow- 
pack to retain a certain amount of free liquid water, is 
approximated by specifying the free v/ater holding capacity 
of the snowpack as described in the user documentation. 
When the free water holding capacity has been satisfied by 
a combination of rainfall plus snowmelt water, runoff is 
calculated using the existing SWMM RUNOFF methodology. 
The effect of the snowpack on the runoff can also be taken 
into account by considering the Manning's 'n' surface 
roughness coefficient as a calibration parameter. As 
indicated previously, it should be noted that additional 
research is required in order to model more exactly the 
snowmelt plus rainfall runoff phenomena. 

A detailed description of input data and formats is given 
in the User's Manual, Volume III. 
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(b) The runoff quality routines used in SWHM have been 

modified to describe the accumulation and washoff of lead 
and chlorides associated with the snowpack. The 
Snowmelt Quality routines have been integrated with the 
Snowmeit Quantity model and now form a user option in 
the SWMM RUNOFF block. 

8 . i* . 2 Testing and verification of the Anderson Snowme 1 1 Quan t i ty model 

8. '4-2.1 Testing with lysimeter data . The Snowmeit model was initially 

tested using data from the Mer Bleue research site in Ottawa. The 

hourly data available at this site for the winter season of 1973-197'+ 

include snowpack, melt measurements from a lysimeter with a surface 

2 
area of about 13 m , precipitation and air temperature. 

The most critical melt period which occurred during the 1973" 
197't period was the spring snowmeit of March k to 7, inclusive. For 
the purposes of the simulation, it was assumed that equation (3) could 
be applied ± three hours before and after any rainfall occurrence. This 
equation was used because the underlying assumptions are the same when 
the time period between rainfall events is relatively short. For 
other melt periods, the degree-hour equation was used with a melt 
coefficient of 0.005 and a base temperature of 33 F. Equation (3) should 
be used for rain on snow events where there may be rainfall preceded or 
followed by a cloudy period. Otherwise, equation (2) would be used. 

The time history of this event is shown on Figure 65- A total of 
3.9^* inches of melt (S-S't inches of melt plus rainfall) was' ca leu 1 ated 
by the model for the period March '* to March 7. For those hours in which 
rainfall occurred, the total rainfall was over six times greater than 
the calculated amount of melt, even though the rainfall intensities were 
re lat i ve ly sma 1 1 . 

8 . A . 2 . 2 Verification of the SWHh Snowmeit Quantity model for the Brucewood 
area. Some winter runoff sampling was undertaken as part of the Brucewood 
data collection program [108]. The data from this area in North York, 
Toronto, were collected during the unusually mild winter of 197'»-75. A 
limited amount of information describing snow distribution, rainfall and 
snowmeit is available. Only three events were suitable for modelling 
purposes: 
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1) March 12, 1975; 

2) March 19, 1975; and, 

3) February 17, 1975- 

The snow distribution at the start of each event was estimated 
from recorded observations. In all three cases the accumulations of 
snow on the ground were relatively small. Meteorological data were 
obtained from the Downsview Meteorological Station (about 10 miles 
away). Unless otherwise stated, the SWMM default values were used 
with Anderson's equation in the simulations discussed below. Only 
the storm of March 12 vjas plotted, as this had the highest proportion 
of snowmelt in the runoff. 

(a) March 12, 1975 (Figure 69) - the recorded and simulated 
bydrographs agreed fairly well. However, the peak flows 
were only about 1 cf s . The average ambient temperature 
was 3^ F and a snow water equivalent of 0.^3 inches was 
assumed. Twenty-five percent of the impervious area and 
85 percent of the pervious area were covered by snow. 
The proportion of snowmelt contributing to the runoff 
was computed as about two-thirds. 

(b) March 19, 1975 - This storm was of a very low intensity 
(maximum, 0.12 in/hr) and the flow rate was less than 
0.5 cfs for most of the duration of the event. Only 10 
percent of the impervious and 25 percent of the pervious 
areas were snow covered. The water equivalent depth was 
taken as 0.1 incnes for the simulation. The ratio of 
computed snowmelt to total runoff recorded was about 1:2, 
Both the computed and recorded bydrographs oscillated between 
0.2 and 0.5 cfs in response to the occasional rainfall 

i nc rement 5 . 

(c) February 17, 1975 - This storm was of fairly low intensity 
and the proportion of snoxvme 1 1 in the runoff was very 
small. The model simulated the peak flows of about 0.5 
cfs fairly well; again oscillations in the computed and 
recorded f I ov/s were evident. Twenty percent of the 
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impervious and 10 percent of the pervious areas were snow- 
covered and the water equivalent depth of the snow pack 
was taken as 0.2 inches. 

In general, the model appears to simulate the peak flows fairly 
well, even when these are very low, but further testing with more 
significant events Is required for complete validation. 

The snowmelt events monitored at Brucewood were all accompanied 
by rain storms of relatively low intensity. In order to test the 
model for higher intensity rain on snow events, the highest intensity part 
of the rain storm recorded on August 2-3, 1975 over Brucewood was used 
as input to the model. A snowpack of 2.5 inches water equivalent 
was assumed to cover the entire catchment. A temperature of 50 F 
and a wind speed of 10 mph were also assumed. Typical results obtained 
from a set of simulations for various catchment conditions are presented 
on Figure 66. Plot 1 represents the hydrograph generated considering only 
the rainfall. Plot 2 shows the rainfall plus snowmelt hydrograph, and 
plot 3 the hydrograph resulting from snowmelt only. The simulation results 
indicated that the model was working as expected and also demonstrated 
that, even for higher intensity rain storms, the snowmelt component can 
account for a significant portion of the hydrograph. 

8.'».2.3 Verification of the SWHM Snowmelt Quantity model for Toronto 
International Airport . In a 1975 study of environmental problems at 
Toronto International Airport [36], a few snowmelt events were recorded. 
Some of these records were used and interpreted In this study. The 
events of March 15 and 16, 1975 were considered suitable for modelling 
purposes. A detailed discretization (3I subcatchments) of the 400 acre 
tributary area was used. The area, which is mainly comprised of runways 
and aprons, is about 90 percent impervious. Snow on these areas is 
ploughed. However, it was assunoed that the melt from the windrows could 
be modelled, since this would flow over the adjacent impervious area 
and into the sewer systems. The degree-hour melt equation was used with a 
melt coefficient of O.OO'tJ and a base temperature of 28 F. (These 
resulted from several calibration runs.) With the expection of the surface 
roughness coefficient in Manning's equation, which was adjusted to 0.026, 
the other characteristics of the area were assigned the SWMM default values. 
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The simulated and measured hydrographs for each of the snowmelt 
events are shown on Figures 67 and 68. The agreement is reasonable 
considering the simplifying assumptions involved. 

8. A. 3 Verification of the SWMM Snowmelt Quality model 

The data used for the verification of the quality model were 
collected at the Brucewood site during the winter of 197^-75. The results 
of the simulations are summarized below. 

(a) March 12, 1975 (Figure 69) " The simulated pol lutographs 
for suspended solids and BOD were of the same general 
shape and order of magnitude as the recorded curves. The 
peak suspended solids concentration was underestimated, 
while the computed BOO was consistently slightly high. The 
initial chloride concentrations recorded were considerably 
higher than those computed, but subsequent values agreed 
fairly well. The high initial concentrations are pro- 
bably a catch basin effect since, follov-^ing road salting, 
there is likely to be some immediate snow and ice melt, 
causing an increase in the CI ion concentration in the 
catch basins. The salt application factor used for this 
simulation was 0.75 tons per mile, which was consistent with 
the available salting records for the area. It should be 
noted that the initial surface chlorides may be dissolved 
in the immediate melt from the roads. 

It was possible to reproduce the recorded lead concentra- 
tions reasonably well, with the exception of the peak 
values (Figure 70). However, this involved specifying 
a very high lead accumulation rate (0-38 lb/100 ft/day) 
to the model. This is probably not a reasonable estimate 
when compared with other pollutant accumulation rates. 
The apparent discrepancy may possibly be due to the localized 
accumulation of lead in the immediate curb area. The 
actual contributing area for lead compounds may be fairly 
small and may contribute only during the first part of an 



18't 



0«:00 

o«n 



oe:oo 



10-00 

» rr— 



IZ'.OO 



14:00 



16:00 



22:00 



E4:00 



02:00 



2S»F 



90"F 



m 



3 



ol 

06:00 






BASE TEMPERATURE 



TORONTO 
INTERNATIONAL AIRPORT 

DRAINAGE AREA I 
SNOWMELT MARCH l5. 1975 
FIG. 67 



SNOWMELT USING EQUATION (2) 



BASE TEMPERATURE 




: ae'F 


MELT RATE 




- 00043 7 hr. 


RIPE CONDITION FOR 


MELT 




TOTAL DRAINAGE AREA 


= 398.5 Ac. 


IMPERVIOUSNESS 




= 90% 


SNO* COVERED IMP. 




•- 80 % 


SNOW COVERED PER. 




= 100% 


INITIAL W.E. ON iMP 




- 04 in. 


INITIAL A.E. ON PER. 




- 0.5 in. 


RECORDED 


VOLUME = 14 3. 100 cu ft 




PREDICTED VOLUME ( SWMM ) 

- 1 7 I, 807 cu. ft. 
DETENTION STORAGE = 0062" 
MANNING'S n ^ 026 



\ 



\. 






08; 00 



10:00 



12:00 I4;00 

CLOCK TIME ( hrs. ) 



1600 



18:00 



20:00 



22:00 



24:00 



02:00 






zo«o 



22:00 



24:00 




02 -OO 



it- 



Q 



0600 



oeoo 



10:00 



TORONTO 
INTERNATIONAL AIRPORT 

DRAINAGE AREA I 
SNOWMELT MARCH 16, 1975 
FIG. 68 



SNOWMELT USING EQUATION 

WIND SPEED 

BASE TEMPERATURE 

MELT RATE 

RIPE CONDITION FOR MELT 

TOTAL DRAINAGE AREA 

IMPERVIOUSNESS 

SNOW COVERED IMP 

SNOW COVERED PER, 

INITIAL WE. ON IMP 

INITIAL WE, ON PER 



f2) 
: S mph 

-' 0.0043 7 hr 

: 398.5 Ac. 
^ 90 % 

- 80 % 
= 100% 
= 016 in. 

- 025 in. 



RECORDED VOLUME = 104,300 cu ft 
PREDICTED VOLUME ( SWMM ) 




= 101,830 cu.fl. 
^DETENTION STORAGE = O062 \ 
^ MANNING'S n = 0.026 / 



12:00 14:00 

CLOCK TIME ( hri > 



16:00 



18:00 



2000 



2200 



24:00 



02 00 



a 
d 



t-Ofc" 



0.8 014 



0.6 



I 0,4 



0.2 



E 120 



g 8o^ 



o 

z 

Si 
3 



40 



20 



- iSf 



10 




1300 



^ 1100 



u 900 

o 

S 

3 

5 700 



500 



300 



I.IZ 



Z£ 



SNOW COVERED IMP. 



i^ SNOW COVERED PER. 
AMOUNT OF SNOW 
INITIAL FW 

F W H C 

IMP. DETENTION STORAGE 

Na OF DRY DAYS 

CLEANING FREQUENCY 

C.B. BO.D. 

CB. VOLUME 

SALT APPLICATION 

TEMPERATURE 

WIND SPEED 

(ANDERSION'5 EQUATION ) 



itoo 



: 25 % 

: 85% 

= 04 



HO 



400 



_ 




500 



600 



TOO 



RAIN a SNOWMELT SIMULATION 
BRUCEWOOD CATCHMENT 

STORM OF MARCH 12, 1975 

FIG. 69 



RECORDED FLOW 



COMPUTED FLOW 



RECORDED SS 

COMPUTED SS 
( Its =0 ) 




RECORDED BOO. 

COMPUTED B.O.D. 
{ Itt = O) 



RECORDED CHLORIDE 
COMPUTED CHLORIDC 




100 



200 



300 400 

TIME ( min.) 



500 



600 



700 

187 



C9 



1.2 



U 



1.0 



0.9 



CONCENTRATION 
o o 






a 06 

< 

-I 






0.5 


/ 

/ 


A n 


0,4 




^^ 


0.3 




V I 


02 


r 




Q.I 


1 












100 



200 



300 



RAIN a SNOWMELT SIMULATION 

BRUCEWOOD CATCHMENT 

STORM OF MARCH 12, 1975 

(LEAD SIMULATION) 
FIG. 70 



RECORDED LEAD 
COMPUTED LEAD 




400 

TIME { mln ) 



500 



600 



700 



800 



event. Runoff from more remote areas of the subcatchment 
may not contain significant quantites of lead. Consequently, 
the total amounts of lead on the catchment compared with 
those of other pollutants may be fairly small, but still 
contribute significantly to runoff contamination. The 
model may be used to compute peal< concentrations of lead, 
but should not usually be applied for determination of 
total loads of this constituent. It is possible that 
additional washoff equations may have to be incorporated 
in the model when a greater data base is available. 

(b) For the remaining simulations (March 19, 1975 and February 
17, 1975), suspended solids were reproduced fairly well. 
BOD was overestimated in both cases; the recorded values 
may have been suppressed by the high chloride concentrations 
measured in the runoff. Chloride levels were simulated 
fairly well with calibration being limited to small variations 
in the salt application factor. 

8. 5 Conclusions 

(a) A Snowmelt Quantity model has been selected and integrated 
with the SWMM RUNOFF block. Preliminary applications based 
on a limited amount of data indicate that the model can 
adequately simulate snowmelt and rain-on-snow events. 

(b) A Snowmelt Quality model, based on the SWMM pollutant 
washoff equations, has been integrated with the SWMM 
RUNOFF block. There are Insufficient data for compari- 
sons between model simulations and measurements, and specific 
conclusions cannot yet be established. 

(c) More synoptic sets of data describing snow distribution, 
snowpack quality, rainfall, snowmelt and runoff quality 
should be collected. Further verification of the SWMM 
snowmelt routines should be carried out. 
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9. EQUIVALENT CATCHMENTS (LUMPING) 

9. 1 Genera) 

The costs of setting up and running a SWMM simulation are 
largely determined by the level of discretization used for a particular 
catchment. The number of subcatchments into which the catchment is divided 
should depend on the information the user requires from the simulation. 
If flow and quality information describing conditions in each of the 
major conduits of a system is of interest, then it is necessary to employ 
a fine discretization, defining the subcatchments tributary to each of 
the conduits modelled. Conversely, if the object is merely to simulate 
the outflow from the entire catchment, it may be desirable to reduce 
the time and effort required for simulation through use of a few 
larger subcatchments and conduits (employing a coarse discretization). 
The purpose of this part of the study was to investigate the level of 
detail necessary to adequately represent an urban watershed, and to 
illustrate the effects of reducing the number of subcatchments on 
the accuracy of runoff simulation. 

A methodology was defined to achieve a representative equivalent 
catchment from theoretical considerations. Verification of the procedures 
involved a series of applications on both hypothetical and real areas. 
The results of the simplified methodology were compared with those of 
detailed simulations. The term "lumped" is used frequently in this chapter 
to describe simulations based on equivalent subcatchments. A lumped model 
may be defined as one which does not take into account the spatial 
variation of the parameters that determine hydrologic processes (for 
instance, when a single catchment is used to describe the entire runoff 
producing area) . 

9. 2 Discussion of Previous Work 

A review of previous work relating to the application of SWMM 
revealed that the concept of lumped simulation has been mentioned by 
several authors. 

Metcalf Eddy et al ['tj] obtained a reasonably good simulation 
with a five subcatchment discretization of the Northwood area, compared 
to a more detailed 12 subcatchment discretization. It was concluded that 
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the overall accuracy of the simulation was reduced with simplification of 
the input data. Keeps and Mein [)l] observed an increase in peak flow of 
about 20^, accompanied by an advance in both the rising and falling limbs 
of the hydrograph, as the number of subcatchments was reduced. This was 
attributed to the elimination of the conduit storage available in a 
finely discretized system and it was concluded the degree of discretiza- 
tion has a significant influence on the simulated hydrograph. Shubinsl<i 
and Roesner [5I] discussed the simplified application of SWMM to large 
areas and concluded that this was feasible, although no comparative simula- 
tions were conducted. Jewell et al [llO] observed an increase in peak 
flow and time to peak as the level of discretization at the Greenfield 
test basin was decreased. It was concluded that the model was not extremely 
sensitive to variations in subcatchment size or to the location of 
gutters and pipes. 

This review of previous work indicated that no systematic study 
of alternate methods of aggregation had been conducted, and that the use 
of longer time intervals in the rainfall hyetograph had not been considered. 

9.3 Basic Concepts 

9.3.1 Overland flow 

Overland flow computations are made in the RUNOFF block of the 
SWMM. A provision is made in this block for minor flow routing computations 
in gutters and small pipes. The storage available in these is relatively 
small and, consequently, their effects were not included in the lumped 
applications (except in the smallest study area). The discussion of 
routing in the TRANSPORT block appears in the next section. 

The overland flow hydrograph for a subcatchment is determined 
in SWMM over a series of time steps during each of which the outflow 
rate is averaged. This rate is computed from Manning's equation, with the 
average ground slope of the subcatchment taken as the slope of the hydraulic 
grade. The excess depth of water (following subtractions for infiltration 
and detention storage) is assumed constant over the flow plane for a given 
time step and is used as the hydraulic radius for the flow plane. The 
following equations are used: 
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V = • De • S ( ) 

Qw = V-W-De (2) 

where: V = veloci ty 

n = Manning's roughness 

De = excess water depth overflow plane 

W = width of the flow plane 

S = slope of the flow plane. 

The equation of continuity is solved for the depth of water remaining over 
the flow plane at the end of each time step. 

De{2) = De(l) - -^ • At | (3) 

where: A = subcatchment area 1 

Qw = outflow rate during the time step 

At = duration of time step 

De{l) = excess intermediate depth after accounting for infiltration 

De(2) = excess water depth resulting from rainfall, infiltration 
and outflow during the time step. 

The parameters controlling the overland flow from pervious or impervious 
areas are: 

- ground slope; 

- Manning's roughness; 

- retention depth; 

- inf i I trat ion; and 

- width of overland flow. 

Of these parameters, all except "width" are assumed to be 
Spatially constant over the area of the subcatchment. Consequently, when 
aggregating subcatchments , it is necessary to assign these spatially 
constant parameters values representative of the aggregated area as a 
whole. The "width" parameter, however, is defined as that dimension 
of a subcatchment across which overland flow travels when draining to 
the inlet manhole or gutter draining that area (See Figure 71)- It is 
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SCHEMATIC SWMM SURFACE 
DRAINAGE ARRANGEMENT 

FIG. 71 




qL= RATE OF OVERLAND FLOW/ UNIT WIDTH 
W = 2L = TOTAL WIDTH OF OVERLAND FLOW 
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important to remember that this is not necessarily the physical width 
of the subcatchment itself. The "width", therefore, is a measure of the 
ease with which surface runoff may discharge from the subcatchment, and 
is a function of the drainage density and number of inlets. From 
equation (2) it may be seen that the rate of outflow in any time interval 
is directly proportional to the "width". A reduction in "width" will 
result in a reduced surface runoff and, hence, cause a greater amount of 
excess water to be built up on the flow plane in response to a given rainfall 
input. Conversely, an increase in "width" will allow water to drain more 
rapidly from the flow plane, and less surface storage will result. This 
effect can be used to attenuate the overland flow hydrograph and to com- 
pensate for conduit storage which is neglected in the lumping process. 
Figure 72 shows the effect of variations in the "width" on the overland 
flow hydrograph for an 11 acre test area. If the correct width of ]^]h 
feet is reduced, the hydrograph is delayed and the peak reduced due to 
increased artificial storage. Conversely, if the width is increased 
the hydrograph is advanced and the peak flow increased. 

The above discussion leads to the conclusion that, if it is desired 
to combine several smaller subcatchments into a single "lumped" catchment 
with the same surface runoff hydrograph, it would be logical to use 
average catchment values for ground slope, roughness, retention depth, 
and infiltration, but to use the sum of the individual subcatchment 
"widths" as the width of the "lumped" catchment. 

9.3.2 Routing 

in order to maintain an accurate simulation of flow through 
a sewer system it Is necessary to specify an equivalent transport system, 
which results in the same attenuation and delay of the inlet hydrographs 
as the actual system. A certain amount of conduit or channel aggregation 
is generally required in most detailed simulations because of computational 
limitations of the SWMM. The conduit dimension and slope are the most 
important properties involved in the SWMM routing computations. The 
general procedure recommended in the EPA SWMM User's Manual is to replace 
several of the smaller conduits with an equivalent element having an average 
size, slope, length, and roughness. It may be anticipated that as 
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fewer subcatchments are employed, more smalt conduits must be neglected, 
with a consequent reduction of conduit Storage in the routing computations. 
Such a reduction will result in higher peak flows and a shorter time 
to peal< in the outlet hydrograph. Two methods may be employed to introduce 
additional storage into the simulation and compensate for the neglected 
smaller conduits. The length of the aggregated conduit may be increased 
to add more conduit storage to the system, or the overland flow "width" 
of the equivalent or "lumped" catchment can be reduced in order to introduce 
the desired attenuation and time delay into the computations. It is 
important also to consider the relative effects of both the overland flow 
routing and sewer system routing. In many small urban catchments, the 
sewer routing effects are less significant than overland flow effects 
and, hence, the outflow hydrograph will not be very sensitive to changes 
in the sewer system. 

9. ^ Procedures and Alternative Methods for Simplified Simulation 

A lumped simulation may involve a single equivalent catchment 
or, in the case of a number of distinctly different drainage subareas, 
may utilize several aggregated areas. This latter case may be necessary 
where upstream controls, storage or diversions necessitate special consi- 
deration and the use of both overland flow and sewer routing computations. 
However, in many cases it would be possible to use the RUNOFF block alone 
to compute an outlet hydrograph. By reducing the lumped catchment overland 
flow width, an extra measure of routing effect can be introduced to make 
up for the sewer routing effects in the TRANSPORT block. In order to 
assess the utility of these concepts, a series of procedures were devised 
for performing lumped simulations. Detailed simulations for different 
areas were available against which the results of lumped simulations 
could be compared. Lumped simulations were prepared according the 
following guidelines. 

S.k.\ Simplified simulation with RUNOFF and TRANSPORT 

(a) Following detailed simulations, the individual subcatch- 
ments v/ere aggregated to provide different levels of discre- 
tization (i.e. detailed, five and one catchments). Para- 
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meters describing surface characteristics {infiltration, 
detention depth, ground slope, Manning's 'n') were 
computed for each equivalent catchment based upon 
subcatchment areas. 

(b) The overland flow width for each aggregated catchment was 
assumed to be the sum of the widths of the individual consti- 
tuent subcatchments . 

(c) The equivalent transport system was established by computa- 
tion of the weighted (by length) average dimensions and 
slopes of the detailed constituent elements. Where a 
single aggregate catchment was modelled, a corresponding 
single equivalent transport element, representative of 

the real trunk system was employed. The length of this 
element was approximated to the length of the actual 
major trunk conduit. 

9.2.^ Simplified simulation with RUNOFF block only 

(a) A single equivalent catchment was aggregated from the 
detailed discretization as before. 

(b) Several equivalent widths were used for testing purposes. 
These ranged from the sum of the subcatchment widths to 
the simple dimension of the lumped catchment, transverse 
to the overall sewer flow. 

(c) No transport element was used. 

These procedures are schematically illustrated in Figure 73- Case I 

indicates a small area discretized into four subcatchments. Overland 

flow is assumed to occur perpendicular to the drainage conduits and the 

corresponding "widths" are shown as broken lines. The "width" would 

be determined for each subcatchment according to the User's Manual. The 

sum of the individual overland flow hydrographs is represented in (A), 

the total surface runoff hydrograph having a peak of Q . Subsequent routing 

of this hydrograph through the TRANSPORT system results in an outflow 

hydrograph having a peak, Q^yj - These have been given peak values of 

60 cfs and 35 cfs, respect ivel /, for purposes of illustration. 



197 



SCHEMATIC ILLUSTRATION OF 

METHODS OF LUMPING 
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Case 2 illustrates the lumping of the individual areas to 
form a single equivalent catchment. In this example, the equivalent surface 
parameters are computed using an area weighted average and the width is 
assumed to be the sum of the four constituent widths, in order to preserve 
the overland flow hydrograph, and is shown as the broken line. The original 
transport elements have been replaced by an equivalent conduit, the length 
of which is approximately that of the major trunk in the finely discretized 
system. The results of the lumped simulation are similar to those of the 
detailed simulation (Case 1). 

Case 3 illustrates a lumped simulation using only the RUNOFF 
block. The overland flow width is reduced to account for the conduit 
storage neglected by not using the TRANSPORT routines. The reduced width 
is indicated by the broken line on the diagram. The diminished width 
causes the overland flow hydrograph to be retarded and attenuated. Selec- 
tion of an appropriate width results in a RUNOFF block surface hydrograph 
which is very similar to that of the TRANSPORT block in Cases 1 and 2. 

The simplifying procedures outlined above were applied on several 
test areas, and are discussed in the following sections. The lumped 
catchment parameters were based upon detailed simulations. In practice, 
the lumped parameters would be developed directly from basic maps and 
plans to represent the entire catchment, and the lumped "width" would 
be an additional calibration parameter. 

9.5 Application of Lumping Concepts 

The simplifying procedures discussed above were applied on 
hypothetical and real test areas. The test areas (Table 16) ranged 
in size from 6 to 2330 acres. Detailed simulations were performed 
for each area and, subsequently, lumped simulations were performed and 
the results compared to those from the detailed simulations. 

9.5,1 Quanti ty effects 

For each of the test areas a number of storm hyetographs were 
selected, covering a range of intensities and durations, in order to 
assess the validity of simulation over a range of hydrologic occurrences. 
A series of synthetic storms were applied over the hypothetical test 
areas, following the work of Brandstetter [7]. Table 17 summarizes 
the storm patterns used in each area. 

w. 



TABLE 16. TEST AREAS FOR STUDIES OF SIMPLIFIED PROCEDURES 



Area 


Level of 
Discret izat ion 


Number of 
Subcatchments 


Numbe 
Rout i 


r of Transport 
ng Elements 


Hypothetical 
t6'i3 acres) 


detailed simulation 
coarse simulation 
single catchment 
single catchment 


35 

5 

i 




38 

30 
1 



West Toronto 
(2330 acres) 


detailed simulation 
single catchment 
single catchment 


45 

1 
1 




73 
I 



Bannatyne 
(5'»2 acres) 


detailed simulation 
coarse simulation 
single catchment 


3 
T 




^2 
3 



Sma) I Wi nnipeg 
Area (6 acres) 


deta i led simulat ion 
single catchment 


6 

1 




6 





TABLE 



STORM PATTERNS FOR TEST APPLICATIONS 



Area 


Pe= 
in 


k Intensity 
inches per hour 


Hypothetica 1 




3.67 

1 .92 

.98 

1.00 

.50 


West Toronto 




2.16 
2.28 


Bannatyne 




.33 
1 .20 


Small Winnipeg 
Area 




5.65 



Durat ion 
in hours 

1 (triangular shape) 

2 (triangular shape) 
k (triangular shape) 
2 (rectangular shape) 
k (rectangular shape) 

2.1 (September 23, 1973) 

.h (October 2, 1973) 

1-0 (May 10, 1973) 

3.3 (June 19, 1971) 

1.5 (September 5, 1971) 

1.0 (5 yr Design Storm) 
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9-5.2 Hypothetical test area;. 

The hypothetical test cirea was designed to have uniform proper- 
ties in subca tchments in order to facilitate the process of aggregation. 
Three levels of discretization were employed. These are indicated in 
Fiugre 7^- Each diagram depicts the complete subcatchment and transport 
system used in the lumped simulations. Five synthetic storms vjere applied 
to the three discretizations of the hypothetical area. Both RUNOFF 
and TRANSPORT blocks were used in these simulations. At both levels of 
lumped discretization the total overland flow width was 79000 ft, the 
same as in the detailed case. The equivalent transport elements used in 
the simplified simulations were determined from an aggregation of the network 
used in the detailed system. 

The results of the individual simulations for the five rainfall 
patterns are shown in Figures 75 to 79. The hydrographs computed by the 
lumped model are very similar to those of the detailed simulations. A 
trend of increasing peak flow with decreasing number of equivalent sub- 
catchments was observed. The greatest increase in peak flow occurred for 
the intense one-hour storm, as indicated in Table l8, 

TABLE 18, RESULTS OF SIMULATIONS FOR FIVE RAINFALL PATTERNS 



Peak Flow Total Runoff 
Ra i nstorm Subcatchment s in (cfs) Volume (ft 3) 

1 hr 37 llOit 3.00x10^ 



5 1163 3.00x10^ 

i 1176 3.00xlo6 



2 hr 37 633 2.75x10^ 

5 629 2. 7^4x106 



1 



6^+7 2.7'+xlO^ 



h hr 37 |S| 2.36x10^ 

5 306 2.36x10^ 



f 310 2.36x10^ 



2 hr rect. 37 ^03 2.62x10^ 



5 '+0'+ 2.62x10^ 

1 J4O9 2.62x10^ 



k hr rect. 37 I63 2.28x10^ 

,5 IM- 2.2bxl0° 



t Mt 2.28x10^ 
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The increase In peak flow is less significant in the case of 
storms of longer duration and lower intensity. This is because, as the 
storm duration increases, steady state conditions occur in which the 
outflow rate is equal to the applied rainfall minus losses. In the four- 
hour storm simulations, a steady state of runoff from the lumped catchment 
was observed. The equivalent transport system provided a fairly good 
approximation of the detailed system in view of the relatively small 
increases in peak flow. These could be eliminated by an appropriate 
increase in length of the equivalent transport system. Such a measure 
could be employed to remedy the minor advance of the entire outflow hydrograph 
in time as fewer subcatchments are employed. This advance, which is 
also attributable to a loss of conduit storage, was limited to 10 minutes 
in the worst case. This would probably not be a significant factor in 
most practical modelling applications. ' 

Additional simplified simulations using only the RUNOFF block 
were conducted in order to investigate the complete elimination of 
the TRANSPORT block. These simulations all used the two-hour triangular 
storm. Results were again compared with those of the detailed discretization 
involving 37 subcatchments. No TRANSPORT elements were used in these 
simulations. The overland flow "width" was reduced to induce more surface 
storage to compensate for the exclusion of all conduit storage. Several 
simulations were necessary to determine the optimum "width". The results 
of these simulations are summarized in Table 19, and the complete hydrographs 
are shown in Figure 80. | 

TABLE 19. SIMPLIFIED SIMULATION USING RUNOFF BLOCK ONLY 
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79,000 




633--- 


2.7^x10^ 


.569x10^ 
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Detailed simulation using RUNOFF and TRANSPORT blocks 
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These results indicate that, in this case, a reduction in 
overland flow width from 79,000 ft to 56,000 ft is necessary to 
compensate for the conduit storage neglected through elimination of 
the TRANSPORT block. There is no direct relationship between conduit 
storage and surface storage and, consequently, it is not possible to 
determine beforehand the exact value of overland flow "width" required 
for an accurate simulation. The peak flow obtained in these simulations 
is relatively insensitive to the reductions in "width" over much of the 
range. A 1^% reduction in "width" to 20,000 ft results in a peak flow 
reduction of only \S%- The total length of trunk conduit in the hypo- 
thetical area was 20,000 ft. The overland "width" which produced the 
best simulation was 2.b times the length of main drainage conduit in the 
real system. | 

These results also illustrate the variation of surface storage 
with width. As the width is reduced, more water is left on the overland 
flow plane at the end of each simulation with a consequent reduction in 
total surface runoff. if the simulation were allowed to run on for a much 
longer time period, most of this volume would be recovered. 

9.5.3 Real test areas J 

The real test areas provided a more realistic set of parameters 
with which to test the lumping procedures. 

9.5.3.1 West Toronto . Three storms were selected for the simplified 
simulations for this area. Initially, a detailed simulation employing 45 
subcatchmen ts and 73 conduits (in TRANSPORT) was conducted for each 
event. Imperviousness and ground slope were measured from maps and plans 
of each subarea, while default values were assumed for infiltration, 
detention and roughness. Equivalent parameters were then computed to 
allow the entire 2330 acre area to be modelled as a single catchment. 
The overland flow width was taken as the sum of ^5 individual widths and 
the equivalent transport element was aggregated from the main trunk 
sewers in the detailed system, which consisted of five major trunks joining 
approximately 1500 feet from the outlet. The largest of these was of III" 
diameter leading to a 96" conduit which approximately bisects the entire 
catchment, and is some 9,000 feet long to the junction with a 75" and a 
51" pipe. The dimensions of the equivalent transport element selected 
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to represent this large trunk sewer, were determined from weighted average 

sizes and slopes of the real trunk. As a first approximation, the actual 

length of 9,000 feet was selected for the equivalent element. The 

results of the simplified simulations are shown in Table 20 and in Figures 
81 to 83. 

TABLE 20. WEST TORONTO SIMPLIFIED SIMULATION 
(RUNOFF fi TRANSPORT) 

Total 
Total Total Surface 
No. of Peak Runoff Infil- Storage 
Sub- Flow Volume tration at end of 
Storm catchments in cf s (f t^) (f t^) Storm (ft^) 

Sept. 23, 1973 '♦S 366 I'l. 59x10^ U.'iSxIO^ 2.19x10^ 

Sept. 23, 1973 1 369 14. '('4x10^ I'4.'i5xl0^ 2.32x10^ 

Oct. 2, 1973 45 A21 8.57x10^ 9-37x10^ 2.25x10^ 

Oct. 2, 1973 1 455 8.33x10^ 9.37x10^ 2.'«7xlO^ 

May 10. 1973 45 710 27.16x10^ 25.'*xlO^ 2.32x10^ 

May 10. 1973 I 9^7 26-78x10^ 25.'*xlO^ 2.6itxl0^ 

The results indicate a generally good agreement between detailed 
and lumped simulations. The total runoff volume computed for the single 
catchment was within 3% of the detailed result. Peak flow rates showed 
a minor increase for the simplified simulations for the storm of 
September 23, and October 2. Initial simulations (using an equivalent 
transport element of 9000 feet) produced a slightly greater increase in 
peak flow and, consequently, the length of the element was increased by 
10%. The accuracy of the simulation of May 10, 1973 appears to depart 
substantially from that obtained in the first two simulations. Figure 83 
shows the results of the simulations for this storm. In the detailed 
simulation many of the upstream pipes were surcharged. This resulted 
in a truncation of the peak flow as indicated (the SWMM simply stores 
flows in excess of pipe capacities). Since only the major trunk line was 
represented in the simplified simulation, there was no surcharging and 
the outflow was greater. Nonetheless, similar values of total runoff 
volume were computed at both levels of discretization. 
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As observed in the hypothetical area simulations, a slight 
advance of the outflow hydrograph occurred. This is attributable to 
conduit storage which is neglected in the lumped simulations. 

Additional simulations were performed using the RUNOFF block 
alone. A reduction of the overland flow width was used to compensate 
for the loss of conduit storage. Several simulations were done for 
the September 23 storm to find optimum value of equivalent width, and 
the results are shown in Table 21 and Figure 8^4. 

TABLE 21. WEST TORONTO SIMPLIFIED SIMULATION 
(RUNOFF ONLY) 

Total 
Overland Total Total Surface 

No. of Flow Width Peai<. Runoff Infil- Storage 
1973 Sub- All Subc's Flow Volume tration at end of 
Storm catchments (ft) in cfs (f t3) (f t3) Storm (ft3) 

Sept. 23 'tS 200,000 366* 14.59x10^ 14.59x10^ 2.19x10^ 

Sept. 23 1 130,000 k2k l^t. 57x10^ lA. 45x10^ 2.23x10^ 

Sept. 23 1 '40,000 319 13-99x10^ ]k.k5^]0^ 2.66x10^ 

Sept. 23 1 60,000 373 l^t. 29x10^ l^t.^SxlO^ 2.4^4x10^ 

•■ Detailed simulation using RUNOFF and TRANSPORT blocl<s 

These results emphasize the sensitivity of the results of the 
single equivalent catchments simulations to variations in overland flow 
width. A reduction in "width" to 130,000 feet was insufficient to 
account for the neglected conduit storage and the lumped peal<. flow was 
too high. Two further simulations were required to achieve a close 
agreement. A "width" of 60,000 feet, or roughly one-third of the 
detailed simulation overland flow "widths" produced acceptable results. 
This was 1.7 times the length of the main trunk conduit in the actual 
system. The total runoff volume decreased with decreasing width, which 
is to be expected since additional surface storage has been introduced. 
This is confirmed by the fact that the total surface storage remaining 
at the end of the simulation increased as the width was reduced. The 
total infiltration was about ]% less for the lumped simulations, 
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indicating that the aggregated catchment parameters closely represented 
those used in the detailed discretization. 

9. 5- 3-2 Winnipeg areas . Simplified modelling procedures were applied 
to two test areas in Winnipeg: Bannatyne (5't2 acres), and a small 
subsection of a neighbouring drainage area (6 acres). Simulations were 
performed to investigate the effect of lumping on smaller catchments in 
which conduit routing effects were relatively minor. Detailed simulations 
were conducted for both areas, and lumped simulations using both the 
RUNOFF and TRANSPOf^T blocks were subsequently performed. Equivalent 
parameters were computed for the single catchment used to represent 
each test area. The TRANSPORT block was used for flow routing for the 
Bannatyne district, while the GUTTER routine of the RUNOFF block was 
employed for this function in the small, six-acre area. The single 
equivalent transport element was based on the length and size of the main 
drainage conduit in each area. The results of the lumped simulations 
are shown in Table 22, and on Figures 85, 86 and 87. 

TABLE 22. WINNIPEG SIMPLIFIED SIMULATION 
{RUNCFF & TRANSPORT) 
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These results indicate a good agreement between the simplified 
and detailed simulations in total runoff volue computed. As experienced 
in previous simulations, there was an increase in peak flow in the lumped 
simulations, primarily due to neglected conduit storage. 
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9.5.3-3 Qual i ty effects . In the SWMM, storm water quality is computed 
from an exponential decay surface washoff equation and simple mass balances 
to establish concentrations within the sewer system. The amount of pollutant 
washoff is a function of the rate of storm water runoff, and the initial 
amount of accumulated surface pollutants. Assuming the correct rate of 
storm water runoff is maintained for the equivalent catchment, the 
correct pollutant washoff rate may be maintained by providing the 
weighted average pollutant loading rate and the aggregated gutter length 
in the lumped catchment simulation. This would result in the same total 
amount of accumulated dust and dirt at the start of the storm as in the 
detailed simulation. 

Sensitivity analyses discussed in Chapter 5 have shown that 
pollutant decay within the sewer system is not significant. Consequently, 
outlet pol lutographs computed for both detailed and simplified simulations 
should be similar, even though the transport system maybe reduced or 
neglected in the coarse discretization. 

Simulations of the hypothetical area, including quality computations, 
were conducted for two storms. Aggregated quality parameters were the 
dust and dirt loading rates and gutter lengths. These dust and dirt 
loading rates were computed using a weighted average rate, based upon 
the area in each of the five land use categories. All five categories 
were modelled in the detailed simulation. An immediate problem was 
encountered in that these loading rates, which are a function of land use, 
are fixed internally in the SWMM and only a single land use may be assigned 
to each subcatchment . Consequently, there was no means to supply the 
appropriate loading rate for the equivalent catchment. The SWMM was 
modified to permit direct card input of the desired loading rates for 
each land use. The aggregate catchment was arbitrarily assigned a single 
land use, and the average loading rate specified. The total gutter length 
was maintained in the lumped simulation, ensuring the same initial 
surface pollutant load as in the detailed simulation. Surface 
runoff quality is also affected by the volume of water stored in catch 
basins. The same number and volume of catch basins was specified in the 
lumped catchment as in the constituent subcatchments . Lumped simulations 
were performed for two synthetic storm patterns. Table 23 summarizes 
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the results of these simulations. Default values were used to describe 
runoff quality parameters in the detailed simulations. 

TABLE 23- SIMPLIFIED QUALITY SIMULATION - HYPOTHETICAL AREA 



Storm 

I 
I 

2 
2 



No. 


of 


Sub- 


- 


catchments 




37 




T 




37 




1 



Total 
BOD 

Computed 
in lb 

5.565 
5,590 



Total 

SS Computed 

in lb 

62, lit i. 
62,l't2 
85.122 
85,190 



These results show a close agreement for total BOD and SS 
washoff at the different levels of discretization. A small increase in 
peak concentrations, and a slight advance of the pol lutographs was 
observed in the simplified simulations. This corresponded to fluctuations 
in the outflow hydrograph noted in the previous sections. 

The foregoing discussion shows that lumping procedures may 
be applied to the surface runoff quality computations without incurring 
significant discrepancies in the resulting poi lutographs . However, 
other factors also influence the storm water quality computations. Dry 
weather flow flushing of sewer deposits are modelled in the TRANSPORT 
block. The deposition and scour effects which determine the first 
flush of pollutants depend upon dry weather flow velocities, conduit 
slopes, and conduit sizes. Adjustment of the equivalent transport 
element of a single subcatchment simulation in order to duplicate 
the detailed simulation first flush effect appears to be unfeasible, 
since such action would affect the outflow hydrograph. The amount 
of solids deposited in conduits at the start of a storm is sensitive to 
conduit slope. If only the RUNOFF block is utilized, these dry weather 
flow effects cannot be modelled. Therefore, it is felt that it is 
not possible to model dry weather flow quality effects in a lumped simula- 
tion. 
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9.6 Temporal Effects of Aggregation 

The successful use of the SWMM in a spatially aggregated mode 
questions the necessity for short interval rainfall data and short computa- 
tional time steps. tt is normal practice to match the input rainfall 
time step with the computational time step. The lumped SWMM no longer 
performs detailed routing in the sewer system and, consequently, longer 
time intervals might be employed for the input rainfall hyetograph if the 
appropriate RUNOFF and TRANSPORT routines prove stable. The potential 
benefits to be derived from use of a longer time step are: more readily 
available rainfall data (in hourly intervals), and the possibility of 
using the SWMM for continuous simulation based on hourly rainfall data. 

Simulations of the hypothetical test area were performed for 
time steps ranging from five minutes to one hour. The area was modelled 
using a single equivalent catchment including a single transport element. 
The GUTTER routing routine in the RUNOFF block was not used. Table 2'* 
summarizes the results of this series of simulations, and the resulting 
hydrographs are shown on Figure 88. 

TABLE m. HYPOTHETICAL AREA EFFECT OF INCREASING TIME INTERVAL 











Time 


Peak 






Sub- 






Step 


Flow 


Vo 1 ume 


Storm 


catchnen 


■ts 


(min) 


(cfs) 


(ft3) 


U hr triangular 




37 




5 


306 


23.6x10^ 


VII 




1 




5 


310 


23.6x10^ 


tl 




I 




15 


30** 


23.6x10^ 


•r 




1 




30 


292 


23.5x10^ 


II 




1 




60 


278 


23.1x10^ 



The results of these simulations indicate that larger time steps 
may be used to model the storm with a duration longer than the basin time 
of concentration, without a significant loss of accuracy. Both RUNOFF 
and TRANSPORT blocks were used in these simulations, and no computational 
stability due to the use of 60-minute time steps was observed. In the 
event of instability in the TRANSPORT block, the RUNOFF block could 
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be used in isolation for the simplified simulation. Table 2U indicates 
that the total runoff volume computed for the 5"minute and 60-minute 
time intervals is similar. Tiie peak outflow rate displays a decrease 
of approximately 10?; with the larger time step. This result is 
expected since use of larger time invervals attenuates the peaks in the 
hyetograph, giving lower overall intensities. 

In order to further investigate the effect of longer time 
steps over a range of hydrologic conditions, simulations were conducted 
using the two-hour and one-hour triangular synthetic storms. These 
durations were, respectively, equal to and less than the time of concen- 
tration of the hypothetical area. It might be expected that the shorter 
the storm duration, the greater would be the error involved in using a 60- 
minute time step. Table 25 shows the results of these simulations. 

TABLE 25. SIMULATION RESULTS USING HOURLY TIME STEPS 









Time 


Peak 


Volume of 




No 


. of 


Step 


Flow 


Runoff 


Storm 


Sul 


^catchments 


(min) 


(cfs) 


(ft3) 


1 hr 






5 


1176 


2.99x10'' 


1 hr 






60 


438 


2.90x10^ 


2 hr 






5 


mj 


2.7'+xlO^ 


2 hr 






60 


ill 4 


2.33x10^ 


k hr 






5 


310 


2.36x10^ 


k hr 






60 


278 


2.31x10^ 



The results demonstrate that for short duration storms, significant 
reductions in peak flow are caused by using hourly time steps. As the 
duration of the applied storm increases, the relative accuracy of the 
longer time interval simulations improves. Thus, for the one-hour storm, 
the peak flow was reduced by 62%, whereas for the four-hour storm a reduc- 
tion of only 105; was observed. The total runoff volume was reasonably 
close in all simulations. The maximum deviation in runoff volumes com- 
puted for short and long interval simulations was approximately 15^. 
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These results show that for long duration or uniform rainfall 
events, use of an hourly time step produces good results. For shorter 
and more sporadic events, however, use of hourly steps will result in too 
much dampening in the rainfall-runoff process and can result in significant 
inaccuracy. This would be especially true for smaller catchments. 
Simpler models, such as STORM, are better suited to such a purpose. It 
is suggested that, in general, the maximum time steps employed in a 

continuous SWMM application should be limited to 20-30 minutes. 

I 
9.7 Conci usions 

(a) An examination of the equations used in SWMM indicated 
the summation of the overland flow hydrographs generated 

for a number of individual subcatchments would fae similar 
to that generated for an equivalent catchment based on 
the aggregated properties of those individual subcatch- 
ments . 

(b) The principal effects of routing the overland flow through 
a detailed sewer system in the TRANSPORT block are an 
attenuation in the hydrograph peak and a retardation of 
the hydrograph due Lo conduit storage. These phenomena 
may be approximated in two ways: 

(i) extension of the equivalent conduit element used 
in RUNOFF + TRANSPORT simulations; and, 

(ii) reduction of the width for overland flow used 
in RUNOFF only simulations. 

(c) Application of the simplified procedures on both hypothe- 
tical and real catchments substantiated the intuitive and 
theoretical basis for "lumping" and indicated that accuracy 
of simulation can be maintained by careful aggregation of 
subcatchment parameters. 

(d) Practical lumped simulations may be conducted for a given 
catchment by using average values of the runoff parameters 
for the entire catchment. An initial estimate of the 
width of overland flow for the lumped catchment can be 
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taken as twice the length of main trunk sewer in the area 
being modelled. This initial estimate of width is an 
approximation based on the values of 2.8 and 1.7 reported 
in Sections 9.5.2 and 9.5-3- Independent lumping studies 
[ill] have recommended twice the longest dimension of the 
lumped catchment as an initial estimate of "width". Both 
of these approximations may be applicable in any one case, 
depending upon the extent of the drainage system and the 
relation between overland flow and channel routing effects. 
The "width" should be considered as an additional calibra- 
tion parameter, and used to calibrate the tumped model with 
measurements, where possible. 

(e) The lumped SWMM is a planning tool rather than a design 
or analysis tool. It should not be applied for design 
since it does not directly consider the effect of sewer 
surcharging or backwater, which may influence storm 
outflows. 

(f) Simulations using an increased time interval for 
hyetograph input and computational methods appear to be 
practical in cases where the time step used exceeds the 
time of concentration of the basin, and where the computa- 
tional stability is preserved. For short intense rainfall 
distributions peak flows will be inaccurate due to averaging 
of the rainfall hyetograph over a longer time interval, 
especially for small watersheds. 



10. A GENERALIZED SWMM QUALITY MODEL 

10. 1 General 

As indicated in Chapter 5, the current state-of-the-art of 
storm water quality modelling rarely justifies complex discretization 
of the catchment or the consideration of in-system storage. It was 
shown in Chapter 9 that is it possible to combine several small subcatch- 
ments into a single equivalent catchment without loss of accuracy in 
storm water quantity and quality simulation. A "lumped" version of 
SWMM, allowing aggregation of subcatchments and variable pollutant 
accumulation control parameters, also shows promise as an intermediate 
tool for multi-event simulation. 

As a next step toward simplified quality modelling, a new computer 
model was developed using the basic concepts of the quality algorithm 
for the SWMH, in a less sophisticated formulation. An attempt was made 
to overcome several disadvantages of the SWMM quality routines as compared, 
for example, with those used in STORM. (This model is described in 
Chapter 11.) In SWMM, the surface pollutant loading rates are built into 
the program, whereas in STORM different loadings rates can be read into 
the model for various land uses. STORM also performs quality computa- 
tions based on a specified input hydrograph, if required. This permits 
a comparison of the pol lutographs obtained from measured hydrographs or 
based on hydrographs from any runoff simulation model. 

The Generalized Quality model represents a practical alternative 
to the use of SWMM with a coarse discretization. The Generalized model 
is easier to calibrate than SWMM because all parameters may be readily 
adjusted. The new model is potentially more accurate for the modelling 
of pollutant washoff than STORM [66] because of the shorter time intervals 
used for washoff computations. The model was tested using data from the 
Brucewood area. 

10.2 Discussion of Basic Concepts 

The pollutant washoff rate is computed from the exponential 
decay equation which is used in both SWMM and STORM. Two different methods 
for the computation of suspended solids are included in the Generalized 
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Quality model. The first (ISS=0) uses the exponential washoff equation 
and allows for a variable availability factor. The second {ISS=1) is 
based on an empirical equation. (Both equations are given in Figure 

35.) 

In the Generalized Quality model, the amounts of surface 
pollutants available from each of the five possible land use areas are 
established from a daily dust and dirt loading rate for each land use 
classification. The proportions of each pollutant in the accumulated 
dust and dirt are specified by the user. The daily loading rate is 
multiplied by the number of dry days and the total surface pollutant 
load is the sum of the contributions from all the land uses. This total 
amount can be modified according to the frequency and efficiency of street 
sweeping as in the SWMM. This approach is equivalent to estimating a 
weighted average pollutant loading for a single equivalent catchment. 

In the case that no local data describing dust and dirt 
loading rates and pollutant composition of dust and dirt are available, 
the values given in the APWA survey may be used as a first estimate [61]. 

10.3 Description of the Generalized Quality Model 

The Generalized Quality model ccmputes surface runoff quality 
fn basically the same manner as the SWMM. Some of the important aspects 
of the new model are summarized below: 

a) aggregated single catchment model; 
b} five land uses possible; 

c) separate input hydrograph forms basis for quality 
computat ions ; 

d) user supplied dust and dirt composition and loading rates 
for each land use; 

e) no quality routing or pollutant decay computations; 

f) no erosion or deposition of sediments computed; 

g) two methods for suspended solids computation; 
h) catch basin contributions modelled; 

i) quality calculations may be based on either flow from 
the total area, or only on flow from the impervious 
area; 
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j) ten pollutants may be simulated (BOD, COD, suspended solids, 
settleable solids, coliforms, N, PO, , CI, lead, oil and 
grease) ; 

k) no default values supplied; and, 

1) pel lutographs , mass curves and surface load statistics are 
available for each pollutant. 

The principal advantage of the Generalized Quality model is 
the reduced cost of computations and reduced data preparation time 
achieved by aggregating the properties of the entire study area and 
modelling i t as a single catchment. The model ts extremely flexible 
and may be used for the study of many aspects of storm water pollution. 

Five different (and uses are allowed in the single aggregated 
catchment : 

single family residentials 
multiple family residential; 
comme r c i a 1 ; 
industrial ; and 
open space. 

Dust and dirt loading rates and compositions are required for each land 
use and the total accumulated amounts are computed according to the 
total area in each category. Generally, only regional average 
pollutant loading rates will be available. However, the user can readily 
adjust these for calibration purposes. The model does not supply default 
values for parameters such as street sweeping frequency and efficiency, 
number of catch basins, number of dry days preceding the storm event 
and pollutant availability and washoff coefficient. Consequently, 
known values or calibrated estimates may be readily supplied and the 
model can be quickly adjusted for local conditions. 

The overland flow hydrograph has to be supplied to the Generalized 
Quality model. This means that this hydrograph may be generated using 
the SWMM or any other urban runoff model, or measured flows can be used. 

Two flow options have been included in the model. In the first, 
the washoff intensity is computed assuming the hydrograph is generated 
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only from impervious areas. This approach is used in STORM for quality 
calculations and would normally be used for low intensity storms, where 
infiltration on the pervious areas predominates and results in little or 
no runoff from these areas. The second method, which would normally 
be used for higher intensity storms or in areas with low infiltration 
rates, assumes the hydrograph to have been generated over the entire area. 

The model does not consider pipe flow or surface erosion prin- 
cipally in order to minimize time required for data preparation and 
computational costs. 

1 . 't Verification of the Generalized Quality Model 

Preliminary calibration and verification of the Generalized 
Quality model was carried out using data from the Brucewood and Barring- 
ton study areas. The detailed SWMM simulations for Brucewood are discussed 
in Chapter 5- The Barrington catchment in East York is served by a 
separate sewer system draining an area of about 56 acres. The principal 
land use is single family residential. The results presented in this 
section demonstrate the flexibility of the model and its application 
in assessing the sensitivity to various parameters in order to establish 
a well-calibrated model. The main parameters used for calibration were: 

a) the availability factor for suspended solids washoff; 

b) the removing coefficient CC; and, 

c) the exponent b. 

The measured dust and dirt accumulation rates and the avai- 
lable values for catch basin BOD were used for the Brucewood simulations. 
Otherwise the SWMM default values were used. 

The recorded storms were of relatively low intensity. In some 
cases the flows generated from previous detailed SWMM simulations were used 
as input to the model in order to compare the results based on these with 
those based on simulations using the measured flows. 

Figures 89 to 93 show the results of the simulations of the 
Brucewood flows. In addition. Figure 3^ , shows a simulation for the 
Barrington catchment. 
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(a) Figure 89 

The empirical equation (lSS=l) was used for suspended 
solids computations. The pol 1 utographs were generated 
based on both the measured and the SWMM hydrographs. 
There was only a marginal difference between the simula- 
tions resulting from use of these alternative flows. 
However, a significant improvement was evidently caused 
by a reduction in the removing coefficient, CC , from 
0.9 to 0.25. 

(b) Figure 90 

The best fit to the measured pol 1 utographs was obtained 
using the empirical equation for suspended solids compu- 
tations (lSS=l) v/ith a removing coefficient, CC, of 0.25. 
The equivalent SWMM results are shown for comparison. The 
accuracy of the Generalized model simulation was not 
significantly altered by use of the measured flows, as 
opposed to the SWMM hydrographs. A removing coefficient 
of 0.9 resulted in large overest imat ion of the peak 
concentrations, v/hiie use of the exponential equation for 
suspended solids (!SS=0) failed to reproduce the shape of 
the measured pol 1 utographs . This was probably due the 
low initial surface loads combined with low rainfall 
i ntens i t ies . 

(c) Figure 91 

The best simulations of the recorded pol lutographs for 
both BOD and suspended solids were obtained using a remov- 
ing coefficient of 0.5 with the empirical equation (lSS=l) 
for suspended solids. It was found that the use of the 
recorded flows instead of the SWMM hydrograph did not 
significantly affect the accuracy of the simulation. 
The equivalent SWMM results are shown for comparison. 
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(d) Figure 92 

In this simulation, the measured dust and dirt accumulation 
rate, during the dry period prior to the event, of 0.5 
lb/day/100 feet of curb was supplied to the model. In 
addition, the suspended solids composition of the dust 

and dirt was set at 830 mg of suspended solids per gm 
of dust and dirt, while the proportion of BOD remained 
the same as the SWMM default value. Good results for 
suspended solids and reasonable accuracy for BOD were 
obtained using a removing coefficient GC = 0.5. 

Once again, use of the exponential suspended solids 
equation lead to a distortion of the pol lutograph , 
emphasizing that for low intensity storms, even with 
an availability factor of 1.0 for suspended solids, 
this equation is not applicable. The results of the 
equivalent SWMM simulation are shown for comparison. 

(e) Figure 93 ' 

In this case, the measured dust and dirt accumulation 
rate prior to the storm was 0.6 lb/day/100 feet of 
curb (c.f. 0.7 in SWMM) with the proportion of suspended 
solids being 570 mg/gm. Computations were based on 
the flows resulting from a detailed SWMM simulation. 
The best results were computed using the empirical equa- 
tion (ISS=1) with a removing coefficient CC = 0,5. 
However, the suspended solids pol lutograph is not very 
wel 1 represented . l 

(f) Figure Sk ' 

These plots summarize the results of simulations for a 
single storm event recorded at the Barrington catchment. 
The low intensity rainfall and low runoff peak are typical 
of the events recorded for this area. The figure indicates 
the results of computations based on the empirical equation 
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(1SS=1} with CC = 0.5, since the low initial surface 
pollutant loads lead to a significant underestimation 
when using the other option {1SS=0). The peak suspended 
solids concentrations computed based on both the SWMM 
hydrograph and the recorded flows were higher than those 
measured, while both methods underestimated BOO levels- 



10.5 Conclusions 



(a) A generalized surface runoff quality model has been 
developed based on the mathematical relationships and 
computational options currently employed in the SWMM 
and STORM models. 

(b) The Generalized Quality model Is more flexible than 

the SWMM quality model as a 1 1 parameters may be adjusted 
without program modification. This facilitates calibra- 
tion and leads to slightly better simulations of surface 
runoff quality than those obtained using SWMM. 

(c) The model can operate with any input hydrograph. This 
work has shown there is little difference between the 
results of the Generalized Quality model based on the 
recorded flows or SWMM hydrograph. 

(d) The shape of the pol lutograph Is generally better estimated 
using the empirical equation (tSS=l), although it appears 
that superior results may be obtained by relating the value 

of the removing coefficient, CC, to runoff intensity. A value 
of CC of 0.25 to 0.5 gave better simulations than those 
based on the SWMM default value of 0.9 for the low 
Intensity storms modelled. 
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ri.: CONTINUOUS SIMULATION 
11.1 General 

Single event simulation using the 5WMM or other detailed nrodels 
provides information on the characteristics of the storm water runoff 
from a catchment for a particular rainfall event. The initial conditions 
must be supplied to the model. These include the number of dry days 
preceding the event, the infiltration potential, the available depression 
storage and capacity of artific'al storage facilities. Often these 
values must be assumed because of insufficient data, and these assumptions 
can often significantly influence the simulated result. The antecedent 
dry period controls the pollutant accumulation {as indicated in Chapter 
5) while the initial infiltration and storage conditions may alter the 
simulated peak flows. From the standpoint of pollution control, statistical 
data describing a large number of events may be more important than single 
low frequency design events. Therefore, it is often necessary to investi- 
gate the long term runoff history and the associated pollutant accumulation 
and depletion. This can be accomplished using a simulation model capable 
of processing a long term continuous precipitation record. From the 
simulated series of events, it is possible to define a critical or "worst" 
overflow event in terms of total pollutant load, total volume, peak flow 
or as some combination of these. 

Several continuous simulation models were reviewed for possible 
interface with SWMM. The STORM model was selected. This is a simple 
continuous simulation model capable of performing both quantity and 
quality computations and assessing the influence of different storage and 
treatment capacities on the number and magnitude of the overflows. This 
model was then tested against both the SWMM model, and against overflow 
records from two test watersheds. Previous studies reported in the 
literature did not present a comparison of the two models. An important 
concept developed in this study was the use of the continuous simulation 
model for screening a long period of runoff events to isolate critical 
events for subsequent SWMM simulations. Information produced in this 
screening process may be used to set up initial conditions for detailed 
modelling. A methodology for the interfacing of continuous and single 
event simulation models is discussed at the end of this chapter. 
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11.2 Review of Available Models 

Continuous simulation models appear to have been first developed 
for rural watersheds, and only recently modified for urban conditions and 
for quality simulation. Several available urban models were reviewed 
for suitability under the terms of this study. Some of these were 
examined in more detail and will be briefly described below. For a 
more detailed discussion of all long term simulation models see Brand- 
stetter [7]. 

11.2.1 Storage, overflow and treatment model (STORM) 

The STORM model was developed jointly by the U.S. Corps of 
Engineers and Water Resources Engineers Inc. [112]. The model computes 
storm water runoff from a single catchment in hourly time steps based 
on the record of a single rain gauge. The rainfall depth in excess of 
the depression storage is transformed to direct runoff through the use 
of a specified runoff coefficient at each time step. Runoff from both 
pervious and impervious areas of the catchment is simulated. Snowmelt 
computations based upon the "degree-day" method may also be performed. 
The water balance between storms is determined via the recovery of 
depression storage based upon specified potential evapotranspi rat i on rates. 

The model performs no routing computations, and all direct 
runoff computed for each time step is assumed to drain from the catchment 
in that time step. Various combinations of storage and treatment 
capacities may be modelled and the effect of these on storm water 
overflows investigated. Quality computations may be performed in each 
time step, based upon the pollutant washoff from different land uses. 
Five common pollutants can be simulated. The quality computations are 
essentially the same as those performed in the SWMM. Dry weather 
flow is not considered. Between storms, the amounts of available 
surface pollutants are modified according to the number of dry days 
for accumulation and the number of street sweepings. 

On the basis of the simulations performed in this study, 
the STORM model is relatively inexpensive to operate and the computer 
cost for processing an entire year of rainfall data is less than $10.00. 



2i»2 



The computer program is nonproprietary and available from the U.S. Army 
Corps of Engineers, Hydrologic Engineering Centre. 

11.2.2 Chicago flow simulation program (F5P) 

The FSP model [113] simulates runoff from a multiple subcatch- 
ment watershed using precipitation data from one or more gauges at any 
time interval. Runoff from both pervious and impervious areas is simula- 
ted. Dry weather flow and snowmelt are accounted for in computing the 
total outflow. The water balance between storms is also considered through 
the use of infiltration and evapotranspi rat ion computations. 

Routing of combined runoff through circular pipes and trape- 
zoidal channels is performed, and reservoir storage at the outlet may 
be modelled. No water quality computations may be performed with the 
model. The computer program is available from the authors, 

11.2.3 Hydrocomp simulation program (HSP) 

The HSP model, Hydrocomp Inc. [lOI], simulates runoff from 3 
mul t i -subcatchment watershed using precipitation data from one or more 
gauges. Short time intervals are usually used, and runoff from both 
pervious and impervious areas is simulated. The runoff computations 
are more sophisticated than those used in most continuous simulation 
models. Infiltration, interception, evapotranspi rat ion , and depression 
storage are computed and the m<x!el may calibrated through adjustment of 
the coefficients for each of these processes. 

Snow accumulation and melt are computed using the Corps of 
Engineers' energy equations. Dry weather flow is also computed and 
combined with the surface runoff. 

Complete routing computations are performed for sewers and 
channel networks using the Kinematic Wave Theory, accounting for 
upstream and downstream controls. Diversion and storage structures may 
be modelled at specific points in the network. Quality computations are 
performed for surface runoff, dry weather flow and constituent routing 
through the network. Up to 17 constituents may be modelled. The 
pollutant accumulation and water balance between storm events is also 
considered. The computer program is proprietary to Hydrocomp International 
Inc. 
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l|.2.'t Massachusetts Institute of Technology urban watershed model (MIT) 

The MIT model [114] simulates storm water runoff for any time 
interval from a mul t i -subcatchment watershed using precipitation data from 
a single rain gauge. Runoff from both pervious and impervious areas is 
simulated. Evapotranspi rat ion and infiltration are modelled during 
each event, but it appears that no water balance between storms is applied. 
Snowmelt is not modelled. Dry weather flow is computed and added to 
the runoff. 

The Kinematic Wave Theory is used to perform routing computa- 
tions in an open channel network, and information on stage and velocity 
is generated. No diversion control or storage facilities are modelled. 
No guality computations may be performed with the model, which is 
proprietary to Resource Analysis Inc. of Boston. 

11.2.5 University of Massachusetts combined sewer control simulation 
model (SCH) 

The SCM program [115] simulates runoff from a mul t i -subcatchment 
watershed using data from a single precipitation gauge. Runoff from 
impervious areas only is computed. Dry weather flow is computed and 
added to the st rm water runoff, but snowmelt is not considered. The 
water balance between storms is accounted for by varying the loss and 
depression storage functions. Hourly time steps are used in the 
computat ions . 

Flow routing is performed using the Dynamic Wave Theory which 
permits consideration of upstream or downsfeam controls, but no storage 
or diversion functions are modelled. No quality computations are per- 
formed in the model. The program is available from the authors. 

11.3 Selection of a Continuous Model 

The selection of a continuous simulation model was based mainly 
on the ability of the model to perform both quantity and quality computa- 
tions, and to simulate the effects of different storage and treatment 
capacities on storm water overflows. Computer time and data requirements 
were also considered. The model selected, after review of all criteria 
shown on the summary Table 26, was STORM. Although at present, this 
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model does not account for dry we;ather flow, a new version is being 
developed which does, and which uses new runoff computations. The 
quality computations in STORM are essentially the same as those of SWMM, 
This is a desirable feature since it facilitates the use of STORM for 
screening a series of events for subsequent SWMM simulations. The 
STORM mode) is well documented, is freely available, and is inexpensive 
to run. 

TABLE 26. LONG TERM SIMULATION MODELS 
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n.4 Overview of the Storm I'^odel 

ll.'t.l Runoff quantity 

In STORM, runoff is computed hourly, based upon the average 
watershed runoff coefficient, the rainfall within the hour and depression 
storage, according to the following formula: 

R = C(P - f) (1) 

where: R = urban area runoff in inches per hour 

C = composite runoff coefficient dependent on 

urban land use 
P = rainfall plus snowme It in inches per hour over the 

urban area 
f = available urban depression storage in inches per 

hour . 
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The runoff generated in each hour is assumed to drain from the water- 
shed within that hour, and may then be modified by any treatment or 
storage option specified. 

n.'j.2 Runoff qua! i ty 

Runoff quality is also computed in hourly time steps. The 
rate of removal of a pollutant from the watershed within each time 
step is assumed to be exponentially related to the amount remaining after 
the preceding step. For each of five pollutants {suspended solids, 
settleable solids, BOD, nitrogen, PO, ) , the relationship is: 

Mp = A P^^j (1-e "^^^^)/At m 

where: M = the amount of pollutant washed off in this time 
P 

step, At 

■Jk] = the availability coefficient 

P, . = the amount of pollutant on the watershed at the start 

of this time step 

R, = the runoff rate from impervious areas 

E = the urban washoff decay coefficient, 
u 

The user may supply the various coefficients or rely on the default values 
in the program. Reference should be made to the User's Manual for a 
more detailed description of these [112]. The amount of pollutant accumu- 
lation on the watershed is governed by the number of dry days, the 
total length of curb and gutter, and the dust and dirt accumulation rate 
on the watershed. The various pollutants are expressed as fractions of 
the dust and dirt (as described in Chapter 5)- The maximum permissible 
amount of pollutant is limited to that accumulated in 90 dry days. Non- 
urban default values are discussed in Volume II. 

11.^.3 Storage and treatment ! 

The computation of the volumes of runoff stored, treated and 
overflowed are based upon the computed runoff and the specified storage 
volume and treatment rate. No storage is used until the specified treatment 
rate is exceeded. If, subsequently, the specified storage capacity is 
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filled, excess runoff overflows, to the receiving water body. When the 
runoff rate is less than the treatment rate, the excess treatment 
capacity is used to empty storage. The pollutant load in the overflow 
is determined from a mass balance of the amount of a constituent in the 
runoff and the amounts in storage and lost to treatment. Storage and 
treatment are treated simply as volumetric functions. 

The program provides a description of each event and a statistical 
summary of all events encountered in the input record. An event is 
considered to begin when the treatment rate is exceeded and end when 
the storage unit is emptied, or when the runoff falls below the 
treatment rate. 

1 1 . ^ . ^4 Other computations and features 

At present, STORM does not simulate the dry weather flow, 
which for small storms may be a significant contributor to outflow 
quantity and quality. The "first-flush" phenomenon is primarily a result 
of the flushing of the accumulated solids, deposited in the sewer pipes 
from the DWF and from the flushing of solids and BOD in catch basins. 
Consequently, it may be expected that STORM will not produce a distinct 
"first-flush" effect. The U.S. Army Corps of Engineers is currently pre- 
paring a new version of STORM which includes dry weather flow computations 
and unit hydrograph techniques [ll8]. 

The shape of the watershed is not considered by STORM, nor is 
the time of concentration taken into account. It is assumed that all 
runoff flows out of the catchment during the time step in which it is 
generated. For larger watersheds, with a concentration time greater than 
one hour, the computed hydrograph will generally occur earlier than 
the observed one. The reverse would be true for smaller watersheds with 
a concentration time less than one hour. This is not usually of great 
concern in many studies, however. 

Snowmelt is computed using the degree-day method, according 
to the formula: 

MELT = COEF x (T - T.^) ClI 
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where: MELT = snowmelt in inches over the basin 

COEF = degree-day coefficient, ranging from 



0.05 to 0.15 inches per degree-day 
T = average daily air temperature F 
T = temperature at which snow begins to melt. 



Snowmelt is computed only for those days when the average daily temperature 
is above the temperature at which snow begins to melt; otherwise the 
precipitation is added to the snowpack. The computed snowmelt is 
distributed uniformly throughout the period, 9 am to 5 pm. 

The Universal Soil Loss Equation [116, 117] is used to compute 
soil erosion. This emperical equation was developed for agricultural 
land. Many empirical coefficients may require calibration. However, 
the method is an approximative technique which may be used to define range 
or limits of watershed sediment yield. Use of this routine is optional. 

The runoff model is calibrated by varying the following runoff 
parameters: 

(a) runoff coefficients for pervious and impervious areas; 

(b) evapotranspi rat ion rates; 

(c) depression storage; and, 

(d) rainfall reduction factor used to relate point pre- 
cipitation to basin average rainfall. 

Because the model considers the entire watershed as a single computational 
catchment, it is quite sensitive to each of these parameters and rapid 
calibration for period totals is generally possible. 

11.5 Application of Storm to Test Watersheds 

STORM was applied to two real test areas on which the SWMM had 
been previously applied for other parts of this study. The objective 
was to verify the STORM results over the longer period against measured 
overflows and to make a comparison between the detailed SWMM simulation 
and the STORM results for selected storms. 
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11.5. 



West Toronto test area 



The West Toronto urban watershed was selected for STORM 
simulation because overflow measurements were available for a relatively 
long period and because detailed SWMM simulations had been conducted. 
The period of April 1, 1973 to October 31, 1973 was simulated by processing 
the hourly rainfall record from the Old Weston Road gauge, situated roughly 
in the middle of the watershed. Quality simulation was not carried out 
since no measurements were available. 

The program was calibrated by mal<ing slight adjustments to the 
depression storage, the runoff coefficient, and the rainfall reduction 
factor relating basin average rainfall to gauge rainfall. The program 
was then run for the entire seven-month period. The results are summarized 
in Table 27- The measurements have been corrected to account for intercepted 
and diverted flows. No storage or treatment was considered in the 
STORM analysis. 

This table indicates a reasonably good agreement between simula- 
tions and measurements. The apparent extra duration of computed over- 
flows could be attributed to the fact that STORM totals are computed 
in hourly increments, whereas measured durations can be accumulated in 
partial hours or minutes. Thus, a one-half hour rainfall may produce a measured 
one-half hour overflow, but STORM would assign this a one hour duration. 

TABLE 27. COMPARISON OF MEASURED AND COMPUTED 
OVERFLOWS FOR THE WEST TORONTO AREA 



Measured 



Computed wi th STORM 



1973 



Volume 
I nches 



No. of 
Events 



Duration 
Hours 



Volume 
I nches 



No. of 
Events 



Durat ion 
Hours 



April 


.55 


7 


17 


May 


1. 17 


11 


28 


June 


.77 


9 


17 


July 


.7'. 


6 


20 


August 


I.3A 


5 


15 


September 


.76 


7 


19 


October 


1.65 


8 


65 



Total 



6.98 



53 



18) 



.69 


7 


^+6 


1.52 


10 


'.8 


.82 


10 


19 


-69 


7 


16 


-99 


k 


11 


.69 


6 


20 


1.95 


10 


67 



7.36 



S^ 



111 
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A comparison between the measured hydrograph, the SWMM computed 
hydrograph (produced under other sections of this study) and the STORM 
computed runoff was made for five selected events, covering a range of 
rainfall intensities and storm durations. The comparisons were based 
on the total runoff hydrograph, that is the measured overflow hydrograph 
adjusted for intercepted flow at the outlet of the West Toronto area 
(see Table 28) . 

TABLE 28. COMPARISON OF MEASURED AND COMPUTED 
RUNOFF FOR THE WEST TORONTO AREA 









Measured 




SWMM 






STORM 




% 


Volume 


"p 


Vol 


1 ume 


% 


Volume 


1973 






CFS 


1 nches 


CFS 


1 nches 


CFS 


Inches 


May 10 






893 


.26 


7^0 




.27 


707 


.'46 


June 22 






703 


.13 


673 




.15 


37'« 


.15 


August 1 


1 




1064 


.69 


816 




.81 


928 


.89 


September 


23 


383 


.19 


367 




.18 


202 


.06 


October 


2 




509 


.12 


^457 




.11 


368 


.18 



Total 



1.39 



1.52 



1.7i« 



The SWMM results generally compare better to the measurements 
than do those of STORM. On September 23, October 2, June 22, and May 10, 
the peak flows computed by SWMM were closer to the measured peaks than 
those of STORM. On August 1, the heavy surcharging computed by SWMM 
resulted in a significant reduction in peak flow, since the SWMM simply 
stores flow in excess of pipe capacity. The total volumes computed using 
SWMM for the five events compares reasonably with the measured volumes, 
but those produced by STORM were higher. This comparison suggests that 
STORM compares fairly well with SWMM in terms of total volume. Figure 
95 shows a comparison between the measured, SWMM and STORM hydrographs 
for these events. 

11.5.2 Bannatyne test area (Winnipeg) 

STORM was applied to this test area for the four-month period, 
June to September 1971- Measurements of both quantity and quality of 
overflow were available. However, it has been noted in Chapter 5 that 
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DWF deposits control the pollutant loads and so the use of STORM for quality 
simulation is 1 lively to be somewhat inaccurate. The area has also been 
modelled using the SWMM, for selected storms during 1971- Table 29 
compares the total observed and computed runoff for each of the four 
months. In this case, a good agreement between the total number and 
volume of overflows was obtained. The computed durations were less than 
those observed. I 

TABLE 29. COMPARISON OF MEASURED AND COMPUTED 
OVERFLOWS FOR THE BANNATYNE AREA 



Measured 



Month 


Volume 




Duration 


(1971) 


Inches 


Number 


Ho 


urs 


June 


M 


t 




m 


July 


1 .17 


12 




fi 


Aug. 


.27 


2 




IT 


Sept. 


.26 


1 




k 



Computed by STORM 



Total 



2.17 



li 



Vol ume 




Durat ion 


1 nches 


Number 


Hours 


M 


i 


20 


m 


ti 


34 


.22 


f 


9 


.21* 

1.92 


J 


_1 
66 



The results of the detailed simulations using SWMM are shown 
in Table 30 along with the corresponding STORM results. Both models 
produced a reasonably good simulation of the runoff volume for six 
events. The total runoff volume computed by STORM agreed with the 
measured total somewhat better than that computed using SWMM. Figure 
96 compares the hydrographs computed by STORM to the recorded flows for the 
three major events in the study period,. i 

Direct comparisons of total pollutant loads were not possible 
because STORM does not consider dry weather flow. However, a comparison 
of the total pollutant emissions computed by SWMM and STORM for surface 
runoff only initially showed large discrepancies. This was due to the 
fact that street sweeping is applied in STORM only when the number of 
dry days exceeds the cleaning frequency. Thus, when using the reported 
frequency of k2 days, no cleaning occurs unless there is a dry period 
in excess of kZ days and, consequently, pollutants build up to maximum 
values on the watershed in a short time. 
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TABLE 30. COMPARISON OF MEASURED AND COMPUTED 
HYDROGRAPHS FOR THE BANNATYNE AREA 





Mea 


sured 




SWMM 




STORM 


Month 


Op 


Volume 


Op 


Volume 


Op 


Vo 1 ume 


(1971) 


CFS 


1 nches 


CFS 


Inches 


CFS 


Inches 


June 19 


kj 


.15 


31 


.15 


37 


.18 


July 3 


60 


.12 


k6 


M 


27 


.05 


July 15 


15 


.02 


8 
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17 


.05 


July 17 


32 


.10 


22 


.10 


27 


.11 


July 28 


35 


.10 


30 


.08 


27 


-09 


Sept. 5 


100 


.26 


95 


.19 




.2^4 



Total 



.80 



.63 



72 



The simulation was repeated with a short cleaning interval, which 
permitted the accumulated surface load to be cleaned more realistically. 
A better agreement between SWMM and STORM surface runoff pollutant computa- 
tions was evident in the simulation with the shorter cleaning interval. 
The amount of pollutants allowed to accumulate was greatly reduced by 
sweeping, resulting in smaller total discharge. 

Table 31 shows the comparison between total pollutant loads 
in the runoff using both the SWMM and STORM. The SWMM simulations were 
based on the actual number of dry days. 

TABLE 31. COMPARISON OF SS DISCHARGES COMPUTED BY 
STORM AND SWMM FOR THE BANNATYNE AREA 



Total 



Suspended Solids in Pounds 



STORM 



168385 



U'lll 



SWMM 



Month 


w 


Day Sweep! 


ng 


13 


Day Sweep i 


ing 




(1971) 


Interval 




In 


terval 




Runoff Only 


July 19 




3816't 






2253 




2710 


July 15 




ejsk 






850 




658 


July 17 




I50it3 






2738 




2112 


July 28 




1A600 






3336 




2931 


Sept. 5 




9382i« 






li23k 




15195 



23606 
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! I . 6 Continuous Simulation Using Lumped SWMM 

The encouraging results obtained using the simplified "lumped" 
SWMM are discussed in Chapter 9. It has been demonstrated that the SWMM 
may be operated in a simplified manner with respect to both area and time 
{larger time steps). The conclusion of these results was that the 
application of the SWMM could be extended from the simulation of the 
single-event to multi-event simulation. Because of reduced data prepara- 
tion, and computer storage and computer time required, many simulations 
of preselected events can be accomplished at a cost comparable with 
that of a detailed single-event simulation. Initial watershed conditions 
are still required for each simulation. 

A further extension of this aspect of the SWMM application will 
be possible when the model is modified to account for the water balance 
between storms. This may be done by re-establishing the soil infiltration 
capacity between storms and also by accounting for the loss of depression 
storage. This will permit the correct antecedent moisture conditions 
to be modelled continuously and will result in appropriate initial 
conditions for each storm in the; continuous record. In order to model 
quality continuously, it is necessary to allow the pollutants to accumu- 
late during dry periods. These modifications would permit the SWMM to 
be operated in a continuous fashion with a reasonable degree of accuracy. 

The advantages of using the SWMM for continuous simulation 
depend upon the level of accuracy required in the results. The SWMM 
provides more detailed storage/treatment options and information 
than does STORM. The SWMM also allows for cost estimation, which 
is an important feature in evaluating alternatives. However, STORM 
provides the basic data which is necessary in preliminary studies of 
ungauged urban systems. Normally, the primary concern of municipalities 
when evaluating their drainage systems is the overflow magnitude and 
frequency. 

Current research at the University of Florida [111] is 
directed toward making these and other modifications to enable SWMM 
to be extended for continuous simulation. Initial results are encouraging, 
and further modification and testing of the continuous simulation SWMM 
should be continued. 
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11.7 



Interfacing Continuous and Single-Event Stmuiatlon Models 



The proposed methodology consists of three stages: data 
preparation, planning stage and design/analysis stage. The data prepara- 
tion stage is based on the Data Analysis Program discussed in Chapter 12. 
The planning stage is based on the STORM model and partly on a lumped and 
modified version of the SWMM (or the Generalized Quality model described 
in Chapter 10, for surface runoff applications). For the design/analysis 
stage, the SWMM or the WRE model is used. The methodology is schematically 
outlined in Figure 97. The description of individual components and of 
their application is presented in the following sections. 

11.7.1 Preparation stage ! 

The application of most urban runoff models usually requires 
3 large volume of data and the cost of preparation of these data can 
amount to a significant portion of the total study costs. Therefore, it 
is desirable to simplify and computerize this part of a study to the 
maximum possible extent. CI imatological data are typically the most 
voluminous data involved in urban runoff studies. Their processing has 
been computerized by means of the Data Analysis Program (DAP) 
described in Chapter 12. The DAP serves as an interface between the 
existing data ban!<;s and both the planning and design stage models. For 
the planning stage, hourly precipitation and temperature data are 
required. These are available from the Atmospheric Environment Service 
(AES). 

The output from the DAP model consists of the punched cards 
required for STORM and SWMM models, and of event summaries. These 
summaries list the clock time of the onset and end of each storm, its 
duration, the total depth of rainfall, the peak intensity and the 
antecedent dry period. The event summaries are useful for a rapid 
review of the precipitation data and, eventually, for the identification 
of critical rainfall/runoff events, where this task is not done using 
the STORM model . 

The remaining physical and operatibnal input data describing 
the runoff model parameters are processed manually. In order to assess 
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the importance of individual parameters In the SWMM model, sensitivity 
analyses and level of catchment discretization were investigated (Chapters 
k and 9). 

11.7-2 Planning stage 

In the planning stage, various land use alternatives, drainage 
systems and the resultant pollutional impacts on the receiving waters 
are evaluated. Typically, at this stage only limited information 
describing the watershed is available and, consequently, a detailed runoff 
simulation would not be feasible or appropriate. At the same time, it is 
important to establish the probability of occurrence of runoff events 
of specific magnitudes. This can be achieved by continuous simulation 
of the urban runoff over a long period using STORM. 

The SWMM model would next be applied in a simplified (lumped) 
manner for a limited number of critical precipitation/runoff events, the 
frequencies and antecedent conditions of which would have been 
determined previously by STORM. If the planning only involves surface 
runoff simulation, the Generalized Quality model may be used in lieu 
of the lumped SWMM for quality simulation. 

At the planning stage, the SWMM can be applied in various degree 
of detail depending on the purpose of the study. Work described in 
Chapter 9 indicates that the SWMM can be applied as a spatially lumped 
model without a significant sacrifice in the accuracy of simulation. 
Reduction of the number of elements in the RUNOFF and TRANSPORT 
blocks to a minimum results in a considerable saving in data reduction 
time and in computer time requirements. Application of the lumped SWMM 
model for the Bannatyne and Toronto West areas showed that it was possible 
to represent both areas as single aggregated catchments, in spite of their 
relatively large areas. 

The long term effectiveness of various alternative storage 
volumes and treatment capacities is broadly assessed using STORM. 
A more detailed analysis is possible using the STORAGE block of the 
SWMM, which can simulate both in-line and off-line storage. The SWMM 
also allows the simulation of various treatment processes. 



Consequently, in a comprehensive planning study the role of 
STORM should usually be limited to that of screening for identification 
of critical events, which are subsequently simulated using SWMM. 

Once the hydrographs and pol lutographs have been simulated 
by the SWMM, different runoff (or overflow) control alternatives may 
be studied. The SWMM can consider seven levels of treatment and 
estimates and, very approximately, the costs of implementing these control 
alternatives. In some instances, a preliminary analysis of the impact 
of critical discharges to the receiving waters would also be carried 
out in the planning phase. 

At the completion of the planning stage, the user has a 
good indication of the nature of the runoff or overflow problem in 
the study area and also has a feeling for the effectiveness of various 
runoff/overflow control measures for a number of critical rainfall/ 
runoff events. The information obtained during this phase of a study 
is at a planning level, where the relative effects and magnitudes are 
more important than the absolute values required for subsequent design 
purposes . 

11.7-3 Design analysis stage 

In this stage, the design of a drainage system and pollution 
control facilities is carried out, as well as a detailed study of 

probable impacts in the receiving waters. Consequently, it is desirable 
to produce accurate hydrographs and pol lutographs for selected events 
using a calibrated, detailed sinulation model. At this level, SWMM 
would be used at a high level of discretization. 

Sewer surcharging would not have been considered at the 
planning stage, as sewers are considered as an open channel network at 
that stage. Critical events for pollution abatement are not usually the 
low frequency storms used in design. However, surcharging, becomes very 
important when analyzing methods to reduce flooding in an existing sewer 
system of insufficient capacity, or when evaluating the response of 
a drainage system to a storm of a frequency lower than the design frequency. 
Under such circumstances, the SWMM is not applicable and other more 
sophisticated models are required. For simulation of backwater effects 
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which might occur in the sewer system, it becomes necessary to use a 
model with dynamic wave routing, such as the WRE model or Dorsch HVM. 

Once calibrated hydrographs and pol lutographs have been 
obtained from the detailed model, the storage and treatprient design can 
be finalized, following an evaluation of the environmental impact of the 
effluent or overflow on receiving waters. 

11.7.'* Case study 

The hypothetical case study used data from the Bannatyne combined 
sewer district for which overflow records for four months were available. 
An abbreviated output from the STORM model is shown in Table 32. A 
visual inspection of the input indicated two storms of particular 
interest - Storm No. II (July 7, 1971) from the quantity point of view, 
and Storm No. 2't (Sept. 5, 1971} from the quality point of view. The 
latter storm not only produced a significant runoff, but also had the 
longest antecedent dry period, 16 days. Such a dry spell would allow 
a considerable accumulation of pollutants prior to the storm. Storm No. 
2k was modelled using the calibrated SWMM. The results of this simula- 
tion are shown in Figure 98. In the same figure, the reductions in the 
SS pollutograph resulting from the control of runoff through storage, 
swirl concentration, and microstrain ing (treatment level four in the SWMM) 
are also shown with the corresponding costs. Finally, it was assumed 
that the catchment discharges into a lake and a receiving water body 
simulation was conducted. Suspended solids were simulated at a point near 
the drainage outlet and the results are shown in Figure 99- For an 
antecedent dry period of 16 days with no effluent treatment, the SS 
concentrations were as high as 6IO mg/1. This value was reduced to '«60 
mg/t for an eight day dry period and to 6O mg/1 for only one dry day. 
In comparison, the level four treatment (mi cros Era i ner) reduced the 
maximum suspended solids concentration in the first case to 210 mg/1. 
This raises some interesting questions regarding other possible 
pollution controls, e.g. interception of the dry weather flow (say three 
DWF) with and without the treatment of the overflow, or flushing of 
sewers with diversion of the flushings to the v/aste treatment plant. 
Such control alternatives could also be studied at the design analysis 
stage. 
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TABLE 31. BANNATYNE DISTRICT QUANTITY ANALYSIS 
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11,8 Surnmary and Conclusions 

(a) A review of available continuous simulation models has 
been carried out, and the STORM model selected as 
being the most suitable model for application under the 
terms of this study. 

(b) STORM has been tested against measured data from two 
watersheds, and against the results of detailed simulations 
by SWMM. The model was found to give reasonably good 
estimates of the overall magnitudes and frequencies of 
overflows. In some instances, close simulation of measured 
hydrographs was also achieved. However, this is not a 
principal objective of the application of STORM. 

(c) For surface runoff, STORM may be calibrated to produce 
similar total pollutant loads to those computed by the 
SWMM over a given series of events. The state-of-the-art 
of quality modeliing is such that emphasis should be placed 
on the overall pollutant discharge rather tlian on instan- 
taneous concentration values. In this respect, STORM may 
be used to provide a useful summary of runoff quality over 
a long period. 

(d) The present 1acl< of a dry weather flow routine in STORM 
precludes any accurate characterization of combined 
sewer pollutant discharges. In this case, a lumped 
version of the SWMM may be used as an intermediate 
planning model following the identification of critical 
events by STORM. It has been noted that a version of 
STORM incorporating dry weather flow will be available in 

1976. 

(e) A methodology for interfacing STORM and SWMM has been 
described and shown to be applicable to a simple 
case study. 
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12. DEVELOPMENT OF A METEOGROLOG tCAL DATA ANALYSIS PROCESS tNG PROGRAM 

12.1 General 

The meteorological data requirements of both the SWMM and STORM 
nodels may be quite extensive, depending upon the number of events simulated 
and the degree of modelling detail. The STORM model usually uses several 
years of hourly precipitation and temperature as a data base, while 
the SWMM requires short inverval (often five minutes) rainfall data for 
each event simulated. Since the SWMM may also be used in a larger 
time step "lumped" fashion for the simulation of many events, organization 
of large amounts of meteorological data will often be necessary. The 
large quantity of these data creates opportunities for both random and 
systematic errors, which may lead to poor simulation and greater debugging 
efforts. 

The terms of reference for this study required that subroutines 
be written for the creation and management of a data base for modelling 
purposes, using data from existing data banks. The most important source 
of meteorological data for real storm simulation is the Atmospheric 
Environment Service (AES), Environment Canada [119]. This chapter 
describes the specific routines which form the Data Analysis Program. 
This is intended to provide an interface between the urban runoff model, 
requiring data in a specific format, and the data bank, in which informa- 
tion is stored in specific formats. 

12.2 Meteorological Data Processing in Canada 

C 1 imatolog i ca I data have been machine processed in Canada since 
about 1950. Prior to this, all recorded data were reduced by clerical 
methods. The basic data input medium is nov/ the punched card, which 
is created from recorded data forms. Punched cards are computer processed 
and thfc d-ita is recorded on rapid access disc or storage devices, while 
t ho cord', -in: retain'.-d as a secondary source fjf information. AES holds 
V i r t u-:! M / .J I I c I imatfj log i ca 1 dat<i collected in Canada. 
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12.2.1 Maintaining the data bank 

The four principal activities [t20J in the monitoring the 
data bank are: 

a) generation of punched cards; 

b) qual i ty control ; 

c) processing; and, 

d) storage. 

Both manual and automatic methods are used to produce punched 
cards containing the information recorded on data sheets. The majority 
of the 350,000 new cards produced per month are punched manually, although 
direct automatic generation of cards from experimental gauges is being 
tested at a number of stations. Quality control programs [121] have 
been developed to check the data stream for accuracy, completeness and 
random errors. Transcription errors, discrepancies and missing data 
are identified and the pertinent records completed manually. At present, 
the homogeneity of the data, for a single gauge or between gauges, is 
not assessed. Consequently, the user must perform mass analyses and 
similar checks to establish the suitability of a particular record or 
group of records. The five operations carried out by the AES for the 
identification of erroneous data are: 

a) The data are scanned and suspect records are Isolated 
for further action. 

b) These data are examined in order to justify their 
accuracy, if possible. 

c) Where accuracy cannot be justified, an attempt is made to 
determine the incorrect procedure used by the observer 
and to correct the data accordingly. 

d) When the remedial measures outlined in the two previous 
steps fail, the recorded value is rejected and is not 
included in further processing. 

e) At this stage, estimated data are entered in lieu of 
rejected or missing records depending on the type of data 
and the necessity for completeness. 
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The edited record Is stored on disc storage devices and eventually 
transferred to the archives, usually in the form of punched cards, although 
the actual discs are being used increasingly for ultimate storage. Data 
is recalled from the archives when required for the preparation of 
reports, tabulations and distribution files to satisfy requests from 
individuals or organizations. In spite of these quality control operations, 
some errors do occur in AES supplied data-. 



12.2.2 



Data formats 



Twenty-one different data formats are presently used by the AES. 
These are referred to as "card types" because of the cards onto which 
all data are originally punched. All the "card types" are described in 
Volume M of this report. Card types 3 and ^ contain the hourly preci- 
pitation and temperature data required for urban runoff modelling 
[122, 123]. Shorter interval precipitation data have not been abstracted 
and digitized into the bank at present. These rainfall records must 
be obtained from photocopies of the original tipping bucket records for 
the desired storms. All data can be obtained from Environment Canada, 
AES. 



12.3 



The Data Analysis Program 



The Data Analysis Programs have been developed to accommodate 
meteorological data from the AES data bank, and to create and manage the 
user data base for the SWMM and STORM models. The three basic functions 
involved are: I 

a) reading the data from the AES bank and screening for 
fau I ty data ; 

b) analysis of the data for consistency and suitability, 
and combination of the records of several gauges; and, 

c) creation of a card data base in the appropriate format 
for the models. I 

Figure 100 shows the functional flow chart for the Data Analysis Program. 
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12.3.1 Reading and screening 

The quality control procedures employed by AES have to be 
supplemented by a data screening process before the creation of the 
user data base for modelling. This function is indicated as block 
"A" in Figure 100. Overpunches are used by AES to denote special 
conditions, such as freezing precipitation. When these are encountered, 
appropriate messages are printed and the datum is replaced with a zero. 
In this manner unintelligible or missing data, which would be rejected 
in subsequent modelling, are removed. In addition, the user may choose 
to replace the rejected data with more realistic values. Frequently, 
gaps occur in the precipitation record due to gauge or recorder malfunc- 
tion. In the absence of card type 3, hourly rainfall data, the Data 
Analysis Program abstracts the six-hour precipitation totals from card 
type k and distributes these equally throughout the missing period. 
(This technique is described in more detail in Volume III of this 
report.) The resulting precipitation record may be subsequently modified 
using techniques described in the following section. These more precise 
techniques should be used when the period of missing data is significantly 
long, such as those resulting from winter shutdown of certain gauges. 
Raw data are screened to identify occasions on which temperature and 
precipitation values fall outside user defined limits. Simple plots may 
be requested to facilitate visual inspection. These methods allow rapid 
detection of any random errors present in the data. 

12.3-2 Analysis and processing of the data 

Following the preliminary screening, it may be desirable to 
combine the records of several gauges to produce a record more representa- 
tive of some intermediate location. This function of the Data Analysis 
Program is shown in Figure 100, block "8". Combination of records may 
be accomplished on the basis of the "normal ratio" method, the "Theissen 
Polygon" method, or by using some weighting factor based upon the distance 
of the gauges from the study area. The consistency of long periods of 
data or the individual gauges should be investigated before combination. 
Single or double mass analysis techniques are available to facilitate these 
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investigations and to determine the significance of variations in the 
relative catch of a gauge over the period of record. A summary analysis 
of the edited and combined precipitation record is generated by the 
model. This is useful for identification of critical events and for 
esbablishing initial conditions for one-event simulation. 

The capabilities provided in the Data Analysis Program may be 
surmiarized as follows: 

(a) Combination of Gauge Records - A variety of techniques may 
be used to combine different gauge records for the purposes of estimating 
missing data or constructing a weighted average precipitation record for 
a specific location. The basic computations are performed using the 
general formula: 



p = 

n 
i = 1 



EX. 
I 



where: P = the combined precipitation 

P. ....P = the individual gauge precipitations 

X X = weighting factors. 

in 

The weighting factors (X.) are user supplied and may be computed by the 
technique appropriate to the desired purpose. Both hourly and short 
interval records may be processed in this manner. 

Missing precipitation data may be estimated from the records 
of three or more neighbouring rain gauges using the "normal ratio" 
method. The long term precipitation totals at the particular stations, 
including the station for which the record is to be generated, are 
computed manually and the ratio of the normal at each station to the 
normal at the subject station is used to determine the relative 
contribution from each record. The weighting factors applied to each 
record are computed from: 

X =Ne ^^ 

Ni ^ ' 



where: Ne = the long-term normal at the gauge for which the estimate 
is required 

Ni = the average long-term normal of the surrounding gauges. 

I 
The average rainfall record for a specific area may be estimated from 

the records of neighbouring gauges. The general combination formula in 

the Data Analysis Program may be used to compute the average precipitation 

record by either the "Theissen Polygon", or "arithmetic average" method. 

For each case the appropriate gauge weighting factors are computed 

separately. 

I 
(i ) Theissen Polygon Weighting 

Xi = Aj_ 
At 

where: Ai = the area of the polygon associated with 

each gauge. 

At = the total catchment area. 

(ii) Arithmetic Weighting 

Xi = 1 (i.e., all gauges have the same weight i ng) 
(iii) Distance Weighting 

Xi = Li 

where: Li = the distance from each gauge to the centre 
of the study area. 

(b) Mass Analysis - Changes in gauge location, exposure, 
instrumentation or manual recording procedures may result in variations 
in the measured precipitation catch from a particular gauge. Such 
changes are often difficult to detect and not analyzed in published 
records. Double mass analysis may be used to check the consistency 
of a gauge record by comparing its accumulated precipitation record 
with the concurrent accumulated average values for several surrounding 
stations. 
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A change in the relative catch of a station may be observed 
when the accumulated records are plotted. The Data Analysis Program has 
been formulated such that a double mass analysis may be performed for any 
selected input station record and a plot of the resulting curve produced 
for visual inspection. 

(c) Event Summaries - It is often useful to preprocess the 
precipitation record before modelling in order to define the antecedent 
conditions for specific events. The Data Analysis Program allows the 
user to supply a storm definition in terms of the minimum amount of 
hourly rainfall necessary to signal the start of a storm event and the 
number of dry hours signifying the end of that event. Thus, the 
program can identify the storms within the given period of record and 
print out summary statistics for each. In some cases, this initial 
analysis is sufficient for identification of critical periods for 
subsequent detailed modelling. The event summaries consist of 
information on: 

- storm index number; 

- start time, end time and duration; 

- antecedent rainfall; and 

- peak rainfall intensity. 

12. 't Conclusions 

(a) A Data Analysis Program suitable for the construction 
and management of a user data base suitable for both 
single event and continuous simulation has been developed. 
The program is operated independently of both the SWMM 
and STORM models to provide the necessary means for 
organizing the large volumes of meteorological data 
required for running these models. 

(b) The program is simple to operate and satisfies an obvious 
practical need. The program can process short interval 
data as well as hourly records and, with simple output 
format changes, may be used to generate precipitation data 
for any urban runoff model. 
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1|, GENERAL CONCLUSIONS AND RECOMMENDATIONS 
tt* I Conclusions | 

(1) An analysis of the available measurements of storm water 
runoff events and combined sewer overflows on urban 
watersheds in Canada indicated that no one set of data 

is perfect for a general validation of all SWMM routines. 
The precipitation and flow measurements from three 
watersheds, of 48, S^2 and 2330 acres were found to be 
acceptable for comparison with the SWMM quantity simula- 
tions for those areas. Pol lutographs were compared for 
both a combined and a storm sewer system. Additional 
information from past studies, and independent studies 
conducted in parallel to the present worl<, was also 
considered. i 

(2) The comparison of measured flows with those simulated 
by the SWMM over a wide range of watersheds of varying 
size and characteristics confirmed that the model (with 
the appropriate level of discretization) predicts peak 
flows, time to peak and total volumes sufficiently accura- 
tely for planning applications and for design work {when 
surcharge is not required). 

(3) The sophisticated WRE TRANSPORT routine, which is being 
incorporated in the SWMM, and the Dorsch HVM mode! give superior 
simulations to those resulting from the use of the initial 

SWMM TRANSPORT routine for surcharged sewer systems. 
The results of simulations for free surface flow obtained 
from all three models are considered equivalent. The 
WRE and HVM routing models were compared for si mi liar 
systems using identical input flows; the results of the 
two models are considered to be equivalent from a practical 
viewpoint. The inlet hydrographs generated by a detailed 
HVM simulation and a less detailed SWMM simulation were 
very s imi lar . 
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('+) The analysis of different urban runoff quality models 

indicates that the SWMM incorporates the most comprehen- 
sive quality routines of those reviewed. However, the 
accuracy possible in quality modelling is not comparable 
with that obtained in quantity modelling, due to the 
greater complexity of the phenomena involved. 

(5) The SWMM quality model, when calibrated, gives estimates 

of the order of magnitude of peak pollutant concentrations, 
and of the total pollutant load released to the receiving 
water, during a storm water discharge event. 

(6) Methods resulting in significant simplifications of both 
the quantity and quality simulation procedures have been 
developed for the computation of catchment outfall 
hydrographs and pollutographs. It has been show that, 
if these methodologies are used, large urban watersheds 
can be modelled in an aggregated manner with a minimum 
of input data, without a significant loss of accuracy 

as compared to equivalent detailed simulations. 

{7) A Generalized Quality model based on the equations used 
in the SWMM has been developed and documented. This 
model requires less data preparation and computer process- 
ing time than the SWMM and may be more rapidly calibrated 
than the corresponding SWMM routine. This model may be 
used in conjunction with flow simulation models or with 
recorded flows. 

{8) Snowmelt quantity and quality submodels have been integrated 
with the SWMM as user options. These models have been 
tested against preliminary field data. 

(9) All of the other SWMM routines have been tested and 

debugged and, subsequently, assessed and/or modified. 
In particular, the storage/treatment cost routine was 
adapted to include cost information based on Canadian 
data . 
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(10) Comparison of the results of STORM with those of detailed 
SWMM simulations and with measurements emphasize that 

the very simple STORM model is an excellent tool for 
predicting the quantity and number of individual overflows 
from a particular district. 

(11) A new model for the processing and analysis of meteoro- 
logical data for use with STORM and SWMM has been 
developed and documented. 

(12) The interfacing of the models investigated has been 
discussed, and a package of models capable of simulating 
all aspects of storm water flow and quality, from the 
planning to the design phase, has been presented. Parti- 
cular emphasis has been placed on the applicability of 
the package to Canadian conditions, as well as on the 
level of sophistication appropriate to the different 
stages of a study and on the selection of the "right 
model for the job". 



13.2 Recommendations 



(1) SWMM flow simulation has proven to be accurate in a 

large number of applications for areas of greatly differing 
size, land use and sewer system configuration, and over 
a range of meteorological conditions. The incorporation 
of the WRE TRANSPORT routine will extend the capabilities 
of the original SWMM to include the analysis of surcharged 
systems. Consequently, it Is considered that the use of 
the SWMM can be recommended for the simulation of flows, 
where required, in most drainage planning and design studies, 
although particular advantages offered by other models 
in certain specific situations should not be overlooked. 

(2) The use of STORM is recommended in planning studies for 
simulating the volume and number of overflows over the 
period of the long-term meteorological record. STORM 
should be used for the generation of a statistical 
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anaTysis of the effects of storage on reducfng overflow 
volumes and for an approximate guide to the probable 
associated reduction in the total pollutant loads dis- 
charged. STORM is particularly useful for the screening 
and comparison of alternatives and for the isolation of 
critical events for subsequent SWMM simulations. 

(3) At present, the modelling of runoff quality cannot be 

expected to be as precise a science as the simulation of 
runoff quantity. In some specific situations different 
from those examined in the report, a limited number of 
flow measurements may be required for quantity simulations. 
However, it is considered that some sampling of the quality 
of runoff and overflows will usually be required to 
supplement the final modelling activities in a study 
directed towards the design of pollution abatement 
facilities. Past experience with the interpretation of 
the results of sampling programs, and subsequent efforts 
for calibration, have emphasized that the collection of 
samples should be very carefully planned in view of the 
high costs involved. Preliminary modelling is often useful 
in determining the specific sampling locations required. 

{k} Simulation of storm water quality will be useful in a 

variety of studies, such as the assessment of the require- 
ment for upgraded treatment plants, the need for new plants, 
the benefits of sewer separation or new interceptors and 
the control of runoff from new developments. 

The costs of treatment and, to a lesser extent, storage 
facilities, are partly a function of the design flows. 
Although the modelling of quality i s by no means exact, 
accurate flow simulation is very useful for the design of 
any pollution abatement facility. Even if the calibra- 
tion of the parameters controlling quality is limited, 
quality simulation can still be of value for illustration 
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13.3 



of the severity of a problem and comparison of alternative 
pollution control schemes. Consequently, it is recommended 
that quantity and quality modelling form an integral part 
of all studies for pollution abatement. 

(5) Any wel i -va I idated computer model which has proven 

suitable for specific location conditions will be a valuable 
tool and an aid to making decisions. The particular advan- 
tages of the SWMM are its wide availability, the excellent 
training programs offered and, most importantly, the 
potential for the study of both quantity and quality 
throughout the urban system and in the receiving waters. 

Implementation of modelling techniques will be most 
successful if the advantages and limitations of different 
models are understood. The interfacing of models outlined 

in this report, such as the complementary use of STORM 
and SWMM is considered an example of the hierarchical use 
of models for different stages of a project. 

Further Work 



It is considered that additional research is needed in the following 
areas : 

a) additional efforts and research for the improvement 
and refinement of existing models; 

b) the further accumulation of data to facilitate these 
improvements; and, 

c) research into, and documentation describing, the range 
of pollution abatement and drainage design alternatives. 

The majority of past research has been concentrated in the first category; 
that is, the tools have been studied to a greater extent than the problemSi: 
However, there are definite needs for further work in the first two 
categories, the improvement of models and the extension of their capabilities. 
These include: 
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- development in the SWMM of a routine to account for the 
water balance between storm events to permit use of simpli- 
fied versions of SWMM for continuous simulation; 

- refinement of the WRE model to account for a wider variety 
of sewer system structures and conduit shapes; and, 

- formulation and testing of new relationships for the 
simulation of quality in the SWMM RUNOFF and SNOWMELT 
routines. 

It is considered that significant improvements in the approach to 
quality modelling will only result from the acquisition of substantial 
amounts of new sampling data from a large number of areas of diverse 
characteristics. 

The extension of the uses of the models investigated in this 
report for a number of additional studies is a most promising area for 
further research. Some examples include: 

- the use of the SWMM for continuous simulation (this 
work is being pioneered by the University of Florida); 

- the application of the SWMM for rural basins, or for 
areas with a low level of urbanization (this work would 
eventually lead to uniform methodologies for all fiooding 
studies and would assist in the more rigorous analysis of 
the effects of urbanization); and, 

- the use of the SWMM for the prediction of inflow into sanitary 
systems and, in general, for the more accurate simulation of 
the hydraulic loading on treatment plants. 

Future research into the solutions to specific problems will 
be more productive from the practitioner's view point than efforts directed 
towards the improvement of the available tool s . For instance, the 
existing versions of the models currently available may be used without 
additional changes for: 

- the examination of new drainage policies and optimum 
methods for flood control; and, 

- the assessment of established design criteria, such as the 
design storm frequency and the design storm itself. 
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In addition, existing models may be used in the quantification of problems 
of broad environmental concern, such as: 

- the relative Importance of different sources of pollution, 
nonpoint sources, rural drainage, urban runoff etc.; and, 

- the most suitable locations for, and patterns of, new 
developments at a municipal and regional scale. 

The widespread acceptance and application of urban runoff models 
within the profession will eventually be better established by the 
savings resulting from modelling techniques, and the achievement of 
well-defined policies, rather than by the exaltation of the theoretical 
merits of particular models. 

13 .'+ implementation on Computer 

The computer tape which resulted from this study contains the 
Storm Water Management Model System in a version which will run under IBM 
Fortran IV Version H, and is presently with the Water Resources Branch, 
Ontario Ministry of the Environment. ' 

All systems' development work was done on one machine and under 
one compiler, which was most appropriate in view of the fact that consid- 
erable effort was required to debug the various programs before effective 
testing could begin. To have done this on several different computer 
systems would have been an expensive task. 

The conversion of a system as large as this for use on other 
computers will require a very considerable effort for each new machine 
(the term "compatible Fortran" is a relative one). 

If the Scientific Authority wants to improve SWMM as new studies 
may require, the problem of maintaining different versions on different 
computers could be very considerable. Problems of portability, compati- 
bility, and maintenance are now being affected by two technological 
advances, namely in in terminal operations, and in communications 
capabilities. By Implementing the system on one central computer and 
permitting people to access it via telecommunications, the problems des- 
cribed are almost wholly obviated. 

The system as it now exists is in a form which will be readily 
available to most engineers or organizations who wish to use it, since 
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it will run on the most common large commercially available computer, 
the IBM 370 series. For example, Datacrown, the Toronto data centre used 
for most of the computer work on this study, has or will have within the 
next few months access to its Toronto computer via terminals in all of 
the Provincial capitals, giving access to the Toronto system with no 
long-distance telephone charges. 

In view of the above, it is our recommendation that the SWMM 
not be implemented in different Fortran versions at this time. 
Instead, it should be made available in a method similar to that described 
above in the version prepared in the study, and consideration should be 
given in the future, if it seems appropriate, to conversion to other 
machines. This would permit close control over the model in the early 
stages of use across the country and permit dynamic modifications to 
be handled with least effort- 
Should it occur that pressure is exerted from commercial 
data processing houses with competitive large scale computer systems for 
the use of the system on their machines, it may be appropriate to make 
the present version available to them at no cost on the understanding 
that they be responsible for conversion to their own system. They might 
then provide a service to users, charging for use of their system but 
not applying any surcharge for use of the program. Alternatively, some 
royalty system might be devised dependent on the manner in which the system 
is made ava i lab 1 e . 

This would ensure that the conversion is made by people most 
familiar with the specific compjter system involved. It would ensure 
that the decision to incur the cost of conversion would be based on 
market potential, and finally it would permit the conversion to be 
done at no cost to the government. 

The relative merits of batch processing, with normal turn 
around in the Toronto area, and interactive time sharing considering the 
various phases of urban watershed modelling, were also considered. 

Some definition of terms is necessary, since usage throughout 
the industry varies somewhat. For the purposes of this study, batch 
processing is defined as a process in which a run is submitted to the 
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computer either through a high-speed device such as a card-reader or 
through a low speed device such as a keyboard terminal, and the job is 
then run and the output returned to the user either immediately or at some 
later time, with no communication between the user and the program 
from the time of submission to the receipt of results. Interactive 
computing is defined as a process which may be initiated in either of 
the ways described above for batch processing, but in which the program 
will halt at certain predetermined points to permit the user to make 
certain decisions before the run is continued. One further term should be 
noted, i.e. conversational computing. This is similar to interactive 
computing, but specific prompting questions are written into the program 
so that a pseudo-conversation can be carried on between the program 
and the user in an effort to prompt or guide the user in the making of 
the run time decisions. The decision as to whether batch or interactive 
computing should be used depends on two issues: the most appropriate 
method of use for the engineer, i.e. whether or not it is important 
for him to be able to interact or converse with the program in terms of 
the design function which he is attempting to carry out; and, the amount 
of data which is required to set up the various runs of the system. 

Regarding use of SWMM, the engineering user does not seem to 
have a great need for interaction, seeming, in our testing, to operate 
perfectly adequately in a batch mode. The amount of data required to set 
up the runs of the SWMM is very large and unsuitable for entry 
via a keyboard terminal. With the system as it now exists, it is 
possible to access the program in an interactive mode via keyboard terminals 
using existing control software owned, for example, by the commerical 
data centre, Datacrown, used for the development work. This could permit 
an engineer to set up a major run in a batch mode and then, having completed 
a first run, modify certain input parameters via a keyboard terminal, 
retain the rest of the data unchanged, and initiate another run. 

To permit conversational computing would require major changes 
to the SWMM Itself and, although such a process is most helpful for 
the neophyte user, it is not considered that It would be worth the effort to 
implement such a system at this time. 
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Accordingly, it is our recommendation that batch processing 
continue to be the prime mode of operation for the system. 
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